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Some Factors Affecting the Surface Area of Hydrated 
Portland Cement as Determined by Water-Vapor 
and Nitrogen Adsorption 


L. A. Tomes,' C. M. Hunt, and R. L. Blaine 


Brunauer-Emmett-Teller (BET) surface areas of hydrated portland cement have been 


calculated from water-vapor and nitrogen-adsorption data. 


The adsorption and desorption 


of water vapor caused measurable decreases in specific surface areas by both water-vapor 


and nitrogen adsorption. Changes 


were also produced by 


wetting and drying. The 


adsorption and desorption of nitrogen at the boiling point of liquid nitrogen did not produce 


similar effects 


The specific surface and nonevaporable-water content of a hydrated cement are also 


somewhat dependent on the initial drying treatment. 


The effect of measuring water-vapor 


surface and nonevaporable water with specimens dried for different periods of time is also 


considered, 
1. Introduction 


Measurements of water vapor ;1, 2, 3a, 3b, 4] ? and 
nitrogen sorption [4, 5], as well as other evidence 
(6, 7), have demonstrated that the hydration products 
of portland cement are largely colloidal in nature. 
Powers and Brownyard [3a, 3b] have carried out 
water-vapor adsorption measurements on a number 
of hydrated cements and found that 7,,, the parameter 
in the Brunauer-Emmett-Teller (BET) equation [8], 
from which a value for specific surface is calculated, 
was proportional to w,, the nonevaporable-water con- 
tent of the hydrated cement. They found the 
ratio of 7,,/w, to be relatively constant and not 
greatly affected by the age or water-cement ratio 
of the specimen. The constant proportionality be- 
tween v,, and w, indicates that surface computed 
from water-vapor adsorption data is proportional 
to the amount of hydration products present. This 
evidence suggests that kind of structural 
similarity exists between the hydration products 
formed during the first few days of hydration, and 
those formed during the later stages of hydration 
{3a}. 

Although the constant ¢,,/w, ratio for specimens 
made from a given cement provides an important 
clue in determining the submicroscopic morphology 
of hydrated portland cement, both ¢7,, [9] and w, 
(3a, 9, 10] are somewhat dependent upon the ar- 
bitrary conditions under which they are determined. 
For example, it was noted by Powers and Brownyard 
j3a] that successive sorption isotherms by the same 
specimen of hardened cement paste were not re- 
producible, a fact that suggests strong interaction 
between water vapor and one or more of the hy- 
dration products of cement. Thus, whereas specific 
surface is a function of the amount of hydration 
products present, nonreversible sorption behavior 
Suggests that it also depends on the history of a 
specimen after hydration has stopped. The present 


some 


TT 
Reseurea Associat repres 
National Bureau of Standards 
Figures in brackets indicate the 


enting the i te “hloride Institut st sot re 


it the end of th 


paper 


HS 1TST 


paper considers the nonreversible sorption behavior 
of hydrated cement and gives further examples that 
illustrate the effect of wetting or exposure to water 
vapor upon subsequent adsorption of water vapor by 
hardened cement paste. Nitrogen-adsorption meas- 
urements have also been made as an independent 
method of detecting small changes in the colloidal 
structure of hydrated cement. 

Another factor to be considered in the determina- 
tion of v, and w, is the difficulty of determining 
how much of the water in a hardened cement is 
chemically combined [8a, 9, 10, 11]. However, by 
carefully controlling the conditions of drying, it is 
possible to obtain a reproducible but arbitrary 
dry weight, which affords a baseline from which 
to determine v,, and w,. Powers and Brownyard 
(3a, 3b] have determined these quantities after drying 
hardened cement paste in a vacuum system con- 
taining magnesium perchlorate dihydrate and tetra- 
hydrate, whereas Copeland and Hayes [9] have 
described a procedure for vacuum-drying cement 
pastes, which uses a dry-ice trap in the system to 
establish the water-vapor concentration. According 
to Copeland and Bragg [12], nonevaporable water 
in specimens dried by this latter procedure prob- 
ably corresponds rather closely to chemically com- 
bined water. However, hardened cement pastes 
and hydrated calcium silicates reach constant 
weight almost asymptotically when dried in vacuum 
at room temperature. It is common practice to 
discontinue drying while the specimens are still 
losing weight at a very small but detectable rate 
(9, 13]. Therefore, this paper also considers the 
manner in which ¢,, and w, vary when they are 
determined on specimens of the same paste that 
have been dried for different periods of time. 


2. Methods and Materials 


Water-Vapor and Nitrogen Adsorption 
Measurements 


2.1. 


Two types of water-vapor sorption apparatus 
were used in these experiments, a vacuum apparatus 
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similar to the MecBain-Bakr Balance [14], and a 
gas-stream apparatus. The essential 
the vacuum apparatus, including one of SIX 
quartz springs, is illustrated in figure | The 
apparatus could accommodate specimens weighing 
0.4 to 0.8 g. Spring extensions were measured with 
an external cathetometer. The gas-stream appara- 
tus, illustrated in figure 2, was similar to the one 
described by Powers and Brownyard [3a]; however, 
it used nitrogen instead of air to transport the water 


its 


vapor. Five-gram specimens were used with this 
apparatus, although it could handle much larger 
specimens. Five to seven days were allowed for 


sorption equilibrium with the gas-stream apparatus, 











whereas 24 hr were allowed with the vacuum ap- 
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Fiat RE 2 Sorption frain for er posing specimens lo controlled 


water-vapor concentration 
4, Bubble counter containing mineral o B, 125-m! Erlenmever flask 
taining specimens; ©, glass wool trap; |), sulfuric acid solutions for establishit 


water vapor pressurt 


features of 


paratus. Exposure of hardened cement pastes 


ach sic 
longer periods of Sione at these relative humiditig ach 


t 
molds a 


has indicated that very little additional Water| 04 
was taken up, and a few specimens actually Jog | the ae 
weight, a behavior previously noted by Powos|! for a 
and Brownyard [3a] Both the vacuum and ras, hydra 
stream apparatus were operated in a cabinet main. | Y = 
tained at 25.0 O.1°C. por | 
Specimens were vacuum dried in an apparaty| ae 


similar to the one Phas ‘ribed by Copeland and Haye! 


experi 
(9), in which a trap maintained at the te mperatuy “The 
of dry ice was placed between the — cimens and thy nt 
ceme 
pump. The system maintained a vacuum of aboy “" ble 
2 w of mercury in the final aan of drying, andj)" ” 
was operated in a room maintained at 21° C. Thy 
vapor pressure of water at 78° © is 0.56 uw. Whe 
the vacuum apparatus illustrated in figure 1 was 
used, final drying and water-vapor  adsorptio; 
measurements were carried out In the same ap- 
paratus. 
Values of ¢,, were calculated by the BET equation 
[RQ] ; CaVU 
) | hh | ALO: 
: it ; P/ Po 50s 
?’ Po P) v,,¢ r,,¢ MeO 
NasO 
: . . mae : KO 
where pis the pressure of water vapor in equilibrium Loss ( 


vapor pressure of pur 
related to the heat of ad- 
amount of water adsorbed 
Because the adsorption of} 26°" 
and 7,, have} The valu 


with the specimen, jp» is the 
water, C’ is a constant 
sorption, and ¢ is the 
per gram of specimen. 
water was measured gravimetrically, 7 


been expressed in grams of water per gram off 
specimen. To compute surface areas, an area] 
coverage of 10.6 A® was assumed for the water] 
molecule. 3.1 
Nitrogen adsorption measurements were made | Tal 
with an apparatus similar to that described by la 
Emmett [16]. Specimens were outgassed overnight we © 
at 100° to 110° © before measuring their nitrogen sever’ 
adsorption. Surface areas were calculated, as- ine ; 
suming each molecule to cover an area of 16.2 —o 
» aay 
2.2. Nonevaporable Water Determination e : 3 
ere 
The values for nonevaporable water appearing in , (TY \ 
figure 6 were determined by igniting the specimens dried 
to constant weight at 1,050° C after they were comb 
vacuum dried in a system containing a trap at the |! thi 
temperature of dry ice. This is essentially the pro- | Sorp! 
cedure of Copeland and Hayes [9], except that drying | Wco 
time has been purposely made variable. The non- weigl 
evaporable water differs from the loss on ignition by dried 
an amount equal to the ignition loss of the original woul 
cement. As Powers and Browny ard have pointed inde 
out [3a], there is some uncertainty in this correction, lt 
which means that the absolute values of nonevap- | &dso 
orable water are known with less certainty than the | Wate 
relative values the « 
| anot 
2.3. Preparation of Specimens ~~ : 
hum 
Neat cement pastes, having a water-cement ratio | dry; 
of 0.5, were prepared in a mechanical mixer and cast | thar 
in the form of square wafers 1 em thick and 5 em on | ¢ape 
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each side. The specimens were removed from the PABLE 2 ——— of wate 7 adsorbed by cement paste during 
is after 24 hr and cured in closed quart jars in several cycles of drying and adsorption 
molds ¢ ‘enw 
' the presence of high concentrations of water vapor 
‘ , : : rh 
for periods indicated in the experiments. To stop » Grams of water adsorbed per gram of 
i( . . : l : “thei dry paste at relative vapor pressure 
hydration, the specimens were vacuum ¢ ried for Specimen Cyek pip 
ydavs. They were broken into smaller pieces and 
stored several months in sealed glass tubes before 0.0703 | 0.1105 | 0.2245 | 0.3300 
sampling. In this way, storage changes were not 
| taking place at a significant rate at the start of the 4 aa a ah ' 4 
xpermments. | 3 0178 0201 0248 . 0283 
experim a : 1. Lumps of paste 5to / { 0169 0190 0236 0266 
The chemical composition ot the clinker and the 8S mm in diametet \ 5 0154 O183 0225 0256 
eement from which the pastes were prepared is given | > aa. aaa — ie 
n table 1. 8 0165 0192 0228 . 0257 
f | 0229 (2h2 0325 0364 
TABLE | Composition of ceme nker * = 0196 0231 0273 0308 
2. Pulverized paste | 3 0173 0198 0241 0271 
less than | mm in } t O169 0193 0229 0260 
liametet 5 0146 o172 0210 0238 
( ompositior 6 O157 O76 0210 0233 
Percentage by weight 7 0140 O158 0193 0220 
Or. by we ht S O1S52 O17 0209 0230 
( ‘ Percentage 
Che specimens were obtained from paste that was cured in moist air for 1 
week. Specimen 1 weighed 0.6636 g and specimen 2 weighed 0.7181 g in their 
Cad 65. 13 $ CaO-siO oi) dryest state (cycle 8 
310 22. 62 $ CaO0-AbO ; ' va)/Wwas; w =Weight of conditioned specimen, wa=dry Weight of specimen 
(bO 1.41 § CaO. AhO>y FeO 1th for respective cycle, and was=dry weight of specimen for cycle 8 
FeO 4 2 { Osi) 2 
sO 0. 
MgO 1.71 ” ” 
VaO 0.04 : 
K2O 2s 
Insoluble i 14 
90 4 
nker W ! ground Ww gypsum a i to a 1.8-percent SO - | Y 
rhe ig , of the blended cement at 1,050° C was 1.69 percent 
Com ilculated from equations four Bogue, The Chemistry of = \ 
Portland Cement, p. 548 (Reinhold Publis! ( I ,» New York, N. Y., 1947 x . 
ies represent potential compositiot um crystallization 80 
o a 
: A 
3. Results 
o 70 ‘a 
; . 3 
3.1. Sorption and Desorption of Water Vapor rs , ey, 
mn , . 4 x S 
lable 2 gives the amounts of water adsorbed by ” — 
' , iel bhi . 60 - i 
wo specimens of paste which were su jected to L A o 
several cycles of drying and exposure to water vapor © Pe 
The experiment was carried out in the apparatus 
illustrated in figure 1 The specimens were dried 05 3 3 4 . “"' i? 9 0 
0 days from the stored state before initial exposure LES OF DRYING AND ADSORPTION 
| ‘ i > ° . ‘le ») | ° co 7 ’ ’ ‘] : ' ° ; 
owatel vapor. In cye le S$ 2 through r the spe cimens Figure 3 S pecifu surface by water-va por adsor ption obtained 
were redried iis nearly as possible to their evele | by repeated measurement of the same specimens of hardened 


dry weight before exposure to water vapor. When 
dried in’ this wav, permanent loss of chemically 
combined water would appear in table 2 as a decrease 
nthe amount of water adsorbed, because all of the 
sorption initially available would not be 
incovered by redrving specimens to their initial dry 
weight. In cycle S, however. the specimens were 
(ried 41 days, so that any initially available sites 
would be uncovered, the specimens had 
indergone change during the experiment. 

It may be seen from table 2 that the specimens 
adsorbed less water with each repeated exposure to 
water vapor for the first 4 or 5 eveles. Thereafter, 
the effects were less marked This Is illustrated in 
another way in figure 3, where calculated surface 
areas have been plotted for each cycle of drying and 
humidification. In the eighth cycle, after the 41-day 
drying treatment, slightly more water was adsorbed 
than in the preceding cycle, but the original sorptive 
capacity of the specimens was never regained. The 
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sites 


unless 


cement paste 


dry weight of the specimens did not decrease by more 
than 2 to 3 mg per gram of paste between the seventh 
and eighth eyeles. There may be uncertainties in 
this estimate, due to gain of small amounts of carbon 
dioxide or other atmospheric contaminants during 
the course of the experiment, but the change in dry 
weight could be much larger than 2 to 3 mg and still 
be too small to account for change in nonevaporable 
water proportional to the change in computed surface. 
Therefore, the loss in sorptive capacity probably 
involves small changes in the colloidal structure of 
the paste. 

The first specimen in table 2 consisted of lumps of 
hardened paste 5 to 8 mm in diameter, whereas the 
second specimen was crushed in an atmosphere of 
nitrogen until most of it was reduced to a powder 
much finer than the 1-mm upper limit given in the 
table. Comparison of the results obtained with 


) 








these two specimens shows them to be slightly differ- 
ent. However, it evident that the amount of 
water adsorbed and the change in sorptive capacity 
during the course of the experiment was not ex- 
tremely dependent on the state of subdivision of 
the specimens. 

Surface areas determined by nitrogen adsorption 
also showed the effect of exposure to water vapor 
In table 3, surface values are given for specimens 
that were measured after initial drying, and speci- 
mens that were exposed to water vapor at relative 
humidities up to 24 percent and redried. The re- 
dried specimens had smaller surfaces, affording 
independent evidence of change in the colloidal 
structure of the paste. There is no evidence that 
the adsorption and desorption of nitrogen at low 
temperatures brought about changes similar to those 
produced with water vapor. 


Is 


TABLE 3 Effect of adsor ption and de sorption of wate 


Papo 


on the nitrogen surface €s of cement pastes 


Measured after 
initial drying 


Measured after adsorp 
tion and desorption 


Curing time of water vapor 


Specimen * Surface Specimen * Surtact 
m-g nd 
1 day 3A 19. 6 3C 11.4 
1 day 3B 18. 0 3D 
1 week 4A 19.9 1C 15.9 
1 week iB 18. 5 iD 16, 1 
1 month. 5A 21.8 5C 17.5 
1 month 5B 21.8 5D 18. 0 


ind D specimens were taken from the same sample of hardened 
uniformity as were samples in the 


» The A, B,C, 
paste, which was not blended to produce 
ensuing experiments 


3.2. Wetting and Drying 


Because the adsorption and desorption of water 
vapor resulted in changes in hydrated cement, the 
effect of wetting and drying was also investigated. 
Table 4 gives the amount of water vapor adsorbed 
by three specimens of paste, one of which was 
measured after the initial drying, and the others 
after rewetting and redrying. The experiment was 
performed with specimens, and adsorption 
measurements were made with the apparatus illus- 
trated in figure 2. Rewetting and redrying were 
accomplished by covering the specimens with water 
for 1 hr at 21° +2°C and carefully reevaporating the 
water in vacuum. The entire operation was _ per- 
formed in the adsorption flasks, so that any ex- 
tracted material would be redeposited with the 
specimens. The results in table 4 have been ex- 
pressed in grams of water adsorbed per gram of 
ignited paste, to permit closer comparison of different 
specimens. The paste from which these specimens 
were obtained was crushed and carefully blended 
before sampling. Analysis of data from 35 hardened 
pastes has shown that the amount of water adsorbed 
by different specimens from the same blended paste 
can be determined with a standard deviation of 0.3 
mg per gram of ignited paste. This reproducibility 
indicates that the results in table 4 are highly 
significant 


5-¢ 


TABLE } Effect of wetting and drying on the subsequent 
adsor ption of water va por by a dried ce¢ ment paste { 


3 


Number ignited paste at relative 


Vapor pres. 


Specimen * of times sure (p/ pr 
wet and 
dri« d 
0.0980 0.1462 0.1946 0.238 
6A 0 0. 0275 0. 0304 0. 0327 y 
6B 1 0265 0293 0316 . = } 
ac ; . 02455 0276 (297 0314 | 


* These specimens were obtained from paste that was cured in 
] week 


Moist air fe 

(‘omparison of table 4 with table 2 shows thg 
wetting and drying produced less change in water. 
sorptive capacity than was produced by adsorption 
and desorption of water vapor. This difference sy. 
gests that the adsorption and desorption process 
may have occurred in a relative bumidity range wher 
shrinkage effects predominated, whereas wetting an¢ 
drying produced strong swelling and shrinking during 
each evele. . 

The nitrogen surface was larger after a single re. 
wetting, but smaller after the process was repeated 
two more times. The results are given in table 5 
It is believed that these changes do not represent a 
spurious effect, because they could be produced with 
other pastes. The change in water-vapor surfae 
after a single rewetting was quite small, as shown in 
table 5. This suggests that a small shrinkage of th 


and 


TABLE 5 Comparison of water vapor nitrogen surface 
areas oft ce ment paste Ss showing the effec fs of wetting and 
drying 

Water vapor Nitrogen 
Number of times wet 
ind dried 
Specimen * Surface Spe Surtace 
me) ¢ mig 
0 HA 97.7 6D 12.7 
l 6B 04, 5 Hk 23. 2 
; ( 87.3 Ht 10 
* These specimens were obtained from paste that was cured in moist air for 

1 week 

Specific surface based on ignited weight of past 


regions accessible to water may increase the size or 
number of channels accessible to nitrogen, whereas 
further shrinkage reduces accessibility to both water 
vapor and nitrogen. Another possible explanation 
of the changes in nitrogen surface is that most of the 
water-soluble materials would probably be dissolved 
from the pores of the paste during the first rewetting 
whereas in subsequent wetting and drying the water- 
soluble materials have already been removed and the 
shrinkage effect would predominate. It 
reasonable that water-vapor sorption would be rela- 
tively unaffected by the presence or absence of water- 
soluble materials in the pores of the paste 

It was anticipated that wetting the specimens 
would lead to further hydration, because they were 
covered with water for 1 hr before drying was started 
However. they redried to weights that were 2 parts 
per thousand less than the original dry weight, and 


secms 
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ue ni 


le the effect was very reproducible from specimen to 
It was doubtful if additional hydration 





specimen. . \ additiona 
eam ,, | took place during wetting, but if it did, it was com- 
tram 9 ( . ) 
°" pletely obscured by structural changes which facili- 
| tated drying. 
7 3.3. v, and w, as Functions of Drying Time 
Ht and w, were determined on specimens that were 


wi | dried for different periods of time. Due to interac- 

“ion between water vapor and hydrated cement as 
oist air! indicated in the preceding experiments, it was neces- 
sary to take a separate specimen for each drying time. 
The specimens were prepared from a single carefully 








S tha : 
wate blended sample ol pulverized hydrated cement paste, 
rption and each specimen weighed about 5 g. The results, 
10 . : : ; =. 
eh which include drying time, dry weight, and amount 
€ Sug. —— satay ' ' 
went of water adsorbed, are given in table 6. All weights 
"ERS ° . ‘ ‘ 
wl ‘| are expressed in grams per gram of ignited cement in 
le ; Pneen shell 
be re order to compare data from different specimens. It 
and 
+ may be seen from table 6 that the specimens which 
"ing -. 
> received longer drying treatments adsorbed more 
le re.) Water. Values of v,, and w, were determined for 
le "e. ° . ° . ° 
each specimen and plotted as functions of drying 
eated . . ; ° : . 4 ¥. 
ble | time in figure 4. It is evident that after a week of 
> oO . " _ 
pe drving, the rates of change of v,, and w, were small. 
‘nt a - jis 
wit! TABLI 6 Dy 7) weights and amounts oT water adsorbed by 
rfaci specimens dried for diffe rent periods of time 
mn in 
f th Dry Grams of water adsorbed per gram of ig 
Cumula weight, nited paste at relative vapor pressure 
tive dry rams per p/p 
irlace ne time ram of 
: nited 
ane " 
past 0. 0999 1478 1954 0. 2400 
Days 
2 6 1. 1373 0.014 172 0.0192 0.0214 
$5 1. 1234 O23 259 (284 0305 
7 1.4 1. 11S 0254 IRs 0304 0329 
7D 4 1.1171 254 284 0314 0336 
6.4 1.115 0264 04 0320 uS41 
go 1145 (27 {O04 O424 0349 
2. 1 1. 1126 1285 ] 1336 0359 
»s These nen ‘ I | pia ured t for 
des 2 days of predryi in hydrat 
15¢- 
oO! 
as . 
el w — 
7 V—v> 
n a nt? WV ry 
dq 
1e a 
“l us 
M4 2 
oO 
4 ° 
cs 
ie = 
is 0.05 
mo 
i- 
— | | 
a m 
ZA 
s / 
l — 
5 5 
: TOTAL DRYING TIME, DAYS 
Ficure 4 and u iw functior drying time 
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LIT : — + —— 
1.16 4 
7A | 

| 

5 a 


anal 


TOTAL SPECIMEN WEIGHT, g/g OF IGNITED PASTE 














Lia 
1.13 + 
1.12 7“ 
aT r | 
0 0.1 0.2 0.3 
P/Po 
Figure 5. Adsorption isotherms up to 24-percent relative 


dried to 


humidity of two specimens of the same cement paste 
different levels of dryness 


The dotted lines represent the dry weights of the specimens and the curves 
represent the specimen weights after exposure to water vapor. 


What is not immediately evident from the values 
in table 6, however, is the fact that the specimens 
did not resorb all of the water that was removed in 
the latter stages of drying. It can be readily shown 
that differences in the amount of water adsorbed by 
any two specimens at the same relative humidity are 
smaller than the differences in their respective dry 
weights. Failure to regain all of this weight may 
be illustrated in another way. Using data from 
table 6, the adsorption isotherms for specimens 7A 
and 7G up to 24 percent relative humidity. are 
plotted in figure 5. The amount of water adsorbed 
has been added to the dry weight, and the values so 
obtained are plotted as a function of relative vapor 
pressure. The horizontal dotted lines 7A’ and 7G’ 
represent the dry weights per gram of ignited paste. 
It is evident that the isotherm for specimen 7G, 
which was dried much longer than specimen 7A, did 
not retrace the 7A isotherm. The isotherms were 
completed in about 40 days, and there is no evidence 
from the rates at which they gained weight that the 
two specimens would have reached the same weight 
in a much longer period. 


4. Discussion 


The water molecule, by virtue of its small molecular 
size and polar nature, is able to penetrate into 
colloidal domains of many adsorbents to which 
nitrogen cannot penetrate. This is true of hydrated 
cement, as shown by the results of Gleysteen and 
Kalousek [4], as well as the values given in table 5. 








The polar and reactive nature of the water molecule 
raises some question as to the mechanism of sorption 
where certain adsorbents are involved [17,18,19,20]. 
This may also be true of hvdrated portland cement. 
However, Bogue [6] has briefly reviewed evidence 
that the surface area of hydrated cement, by water- 
vapor adsorption, is of the magnitude to be expected 
from electron microscope and low-angle X-ray 
scattering measurements. For the purposes of the 
present discussion, it is not essential that either the 
water vapor or nitrogen surface represent absolute 
surface values, to bring out the structurally significant 
fact that water vapor has access to sorption sites in 
hydrated portland cement that are inaccessible to 
nitrogen, or to investigate a few of the factors affect- 
ing the magnitude and reproducibility of surface 
values obtained. 

The present paper considers some characteristics 
of hydrated cement which are important in the 
design and interpretation of gas and vapor adsorp- 
tion experiments. They illustrate some of the 
problems involved in the determination of specific 
surface and nonevaporable water, and they also 
suggest that exposure to water or water vapor 
should not be one of the unintentional or uncontrolled 
variables in an experiment. 

The isotherms of Powers and Brownyard [3a 
mentioned in the early part of this paper, were 
measured at relative vapor pressures extending from 
dryness to slightly less than saturation pressure, and 
the specimens showed permanent gain in weight 
when redried. They attributed differences in suc- 
cessive isotherms to irreversible shrinkage at low 
relative humidity and additional hydration at rela- 
tive humidities greater than 80 percent. In the 
present experiments, the specimens were not exposed 
to relative humidities greater than 33 percent, ex- 
cept those that were rewet and redried.  Further- 
more, it was possible to redry all specimens to 
weights close to their original dry weights. In this 
connection, Powers [21] has presented data showing 
that the hydration of an unhydrated cement was 
virtually unmeasureable when the cements were 
stored at relative humidities below 30 percent. He 
also found that hydration was slight unless the 
cement was exposed to relative humidities greater 
than 70 percent. Thus, in the present experiments, 
changes in specific surface of the order of 10 to 30 
percent of the initial value were produced with little 
if any change in the degree of hydration. 

The changes in surface and p,, obtained here 
might be predicted from Powers and Brownyard’s 
results, or from the changes in structure that are 
believed to take place in other inorganic gels during 


and desorption of vapors [23]. However, this 
would not be predicted from the earlier isotherms of 
which indicated that once _ initial drying 
‘shrinkage had taken place, adsorption and desorption 
of water vapor were reversible [1, 2, 24]. However, 
some differences between Jesser’s experiments and 
the present ones might be pointed out. He measured 
comparatively large specimens of low water-cement 
ratio that were never subjected to vacuum drying. 


Jesser, 
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whereas the specimens in the present experimen gs a fun’ 
were broken into small pieces, or even pulverize | jinear re 
and dried in vacuum. Jesser’s experiments resemble | 
to some extent, conditions to which a CONCretp| ecause 
might be exposed in service, whereas the presen) drving t 
experiments and those of Powers and Brownyay} yhile t! 
employ rather conventional gas sorption technique | obtaine 
[14]. | all dries 

There are a number of mechanisms by which thp fgure 6 
accessibility of an adsorbent to nitrogen or water! 035 m 
vapor might be reduced, and a few examples might} during 
be cited. Rao found ‘‘drift”’ in the w ater-sorption | specific 
isotherms of titania [25] and ferric oxide [26] gels (‘opelal 
He attributed it to coalescence of smaller particles; there t 


the Umi" 


into larger particles, accompanied by a decrease jp | », corr 
total capillary volume. An example that may be! differet 


more clesely related to the colloid chemistry of} in figul 
hydrated portland cement is Taylor’s evidence [27|| the pre 
that calcium silicate hydrate I, sometimes referred 
to as tobermorite [28, 29], possesses an expanding 
layer structure similar to that found in some clay 


4 . . s( 
minerals. He also found that w ater-vapor isobars wr 
of this compound did not return to their starting age 

: the V 


point [30], and suggested that this was due to failure 
of the lattice to reexpand to its original dimensions 
when brought back to the original conditions of 


cemen 
nitrog 








temperature and humidity. Another process by preclu 
which the surface of an adsorbent might change is ically 
polymorphic transformation. Kalousek [31] noted | °° 
changes in the X-ray pattern of one of his caleium| "PT 
silicate hydrate preparations. Also Brunauer, Cope- — 
land, and Bragg [32] reported indications that re- caleut 
moval of lime frem their tobermorite preparation — 
resulted in formation of small amounts of afwillite. seve 
Whatever the mechanism by which the sorptive vk 
capacity of hydrated portland cement changes, the ; 
present experiments suggest changes in sorptiv em 
capacity, and the factors producing these changes | 
are fully as interesting as the surface values | °°? 
themselves. cease 
Another point should be mentioned in connection = 
with the initial drying of hydrated cement paste. | “??' 
If one takes the data shown in figure 4 and plots r, 
Tl 
~_ L. E 
ie ior 
™™ latin 
past 
F autl 
g Oe of 
ads 
O10} : 
>) 
3 
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FIGURE 6. v, as a function of nonevaporable water content of a 


cement paste dried la di Tere ni levels of dy yness 





"er ized | jinear relationship is obtained. 


Might} during the latter stages of drying 
rption} specific surface, 0.35 mg corresponds to 1.2 m?/g 
| gels Copeland and Haves [9] have pointed out that 








‘IMeng gs a function of w, as has been done in figure 6, a 


This line differs from 


plots of Powers and Brownyard [3b], 


the Um/Wn 


Meret} pecause all of the points have been obtained by 
Presen:| drying the same paste for different periods of time, 


while the plots of Powers and Brownvard were 


INigue{ obtained with pastes of different ages which were 


The slope of the line in 
0.35, which means that 7,, increased by 
every milligram of water removed 
Expressed as 


} all dried the same way. 
figure 6 1s 


035 mg tor 


there is no simple way of calculating a change in 
corresponding to a change in w, produced by a 
diferent drying treatment. The linear relationship 


iry of! in figure 6 suggests a possible empirical approach to 


the problem 
5. Summary 


Adsorption and desorption of water vapor and 
epeated wetting and drying resulted in decrease in 
the water-vapor sorptive capacity of a hardened 
cement paste. Changes were also produced in the 
nitrogen sorptive capacity. These results do not 
preclude the possibility of very small losses of chem- 
eally combined water or gains of small amounts of 
atmospheric contaminants; however, they probably 
represent shrinkage or other modification of the 
colloidal structure of cement paste. If surface areas 
calculated from sorption data are to be treated as 
functions of the amount of hydration products 
present, it is important to consider the history of a 
specimen after hydration has stopped 


‘he effect of varying the drying time before 
measuring the nonevaporable water content and 
surface area by water-vapor adsorption was also 
considered, When a hardened paste received suc- 


cessively longer drying treatments, a linear relation- 
ship between 7, and w, was obtained as the paste 
ipproached dry hess. 
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\{n expression for the field of an electric dipole located over a flat surface with a speci- 
fied surface impedance Z is derived from the formal integral solution by a modified saddle- 
point method. Using the value of Z appropriate for a homogeneous conducting ground, 
the general expressions reduce to those given by Norton. In this case the phase of Z lies 
between 0 and 45 When the phase exceeds 45°, as it may for a stratified ground, the radi- 
ated wave of the dipole becomes partially trapped to the interface. This effect is most 
pronounced for an inductive surface where the phase of Z is 90°, in which case the energy 
of the wave is confined within a small distance from the surface. Such inductive surfaces 
are a metallic plane with a thin dielectric film or a corrugated surface. 

This unifying treatment provides a link between the surface waves of Zenneck, Sommer- 
feld, Norton, and Goubau, and indicates that the phase angle of Z controls the extent to 
which these waves may exist for a dipole excitation. 


1. Introduction 


The propagetion of radio waves along the surface of the ground has been discussed from a 
theoretical standpoint for many vears. As long ago as 1907, Zenneck [1] ' showed that a wave, 
which was a solution of Maxwell's equations, traveled without change of pattern over a flat 
surface bounding two homogeneous media of different conductivity and dielectric constants. 
When the upper medium is air and the lower medium is a homogeneous dissipative ground, the 
wave is characterized by a phase velocity greater than that of light and a small attenuation in 
the direction along the interface. Furthermore, this Zenneck surface wave, as it has been 
called, is highly attenuated with height ebove the surface. In 1909 Sommerfeld [2] solved the 
problem of a vertical dipole over a homogeneous ground (half-space). In an attempt to explain 
the physical nature of his solution, he divided the expression for the field into a “space wave” 
and a “surface wave.’ Both parts, according to Sommerfeld, are necessary in order to satisfy 
Maxwell's equations and the appropriate boundary conditions. The surface-wave part varied 
inversely as the square root of the range, and it was identified as the radial counterpart of the 
plane Zenneck surface wave 

lor many vears it was believed that this Sommerfeld surface wave was the predominant 
component of the field radiated from a vertical antenna over a finitely conducting ground. It 
was not until 1935 that an error in sign in Sommerfeld’s 1909 paper ? was pointed out by Norton 
3], which also partly accounted for the unusual calculated field-strength curves of Rolf [4]. 
At about this time there was a series of papers by Weyl, Van der Pol, Niessen, Wise, and Norton, 
deriving more eccurate representations for the field of the dipole. A discussionof this later 
work has been given by Bouwkamp [5]. Norton [6, 7], in particular, has developed his formula 
for the field components to a stage where numerical results can readily be obtained. It is now 
generally accepted that the Sommerfeld surface wave (or the radial Zenneck wave) does not 
bear any similarity to the total field of a vertical dipole over a homogeneous conducting earth. 
In fact, the field excited by a dipole varies as 1/d, where d is the distance for low frequencies, and 
it varies as 1/d® for high frequencies. Norton has suggested that the field in air of a dipole 
over a homogeneous ground be expressed as a sum of three components: A direct ray (or pri- 
mary influence), a reflected ray which is to be modified by an appropriate Fresnel reflection 
coefficient, and a correction term. Norton has described the first and second components as 


the space wave; the third or correction term, the surface wave. This seems to be a logical step, 
Figures in brackets indicate the literature references at the end of this paper 
It has been suggested by B. M. Fannin (Final Report, Pt. III, Investigations of air to air and air to ground experimental data, Elec. Eng 
Dept., Cornell University, 1951) that Summerfeld did not actually make an “‘error in sign’”’ in his 1909 paper. Fannin’s contention is that Rolf and 
ther workers simply misinterpreted Sommerfeld's formula for the attenuation factor. This can hardly be so, as Sommerfeld himself in 1910 (The 


} 


propagation in wireless telegraphy, Jahrb. d. Drahtl. Tel. 4, 157, 1910) used the incorrect formula and its resulting series expansion 
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although, teken separately, the space and surface waves of Norton are not solutions of Max- 
well’s equations. On the other hand, his “space wave” is the contribution that would be derived 
on the basis of geometrical optics, and his “surface wave” is the correction from wave theory. 
This latter term will be called the “Norton surface wave” as distinct from the Zenneck and 
Sommerfeld surface waves, and the trapped surface waves discussed below. 

It was pointed out by Sommerfeld [8] nearly 60 vears ago that a straight cvlindricel con- 
ductor of finite conductivity can act as a guide for electromagnetic waves. Some 50 vears later 
Goubau [9, 10] demonstrated that such a evlindrical surface can be launched with reasoneble 
efficiency from a coaxial line whose outer surface is terminated in a conical horn. The improve- 
ment in the trensmission properties by coating the wire with a thin dielectric film has been dis- 
cussed in detail by Goubeu |10]. The plane counterpart of the Goubau-Sommerfeld cylindrical 
surface wave is obtained when a flat metallic surface is coated with a dielectric film. Atwood {1 1 
has discussed the nature of the surface waves that can exist in a structure of this type. When the 
film thickness is small compered to the wavelength, the single propageting mode has e phase 
velocity slower than that of light and is attenuated rapidly above the surfece in the air. Such 
surfaces heve been called inductive because the normal surfece impedance, looking into the 
surface, is almost purely imaginary for a low-loss dielectric on a well-conducting bese. A similar 
tvpe of surface wave can exist over a corrugated surface [12, 13, 14] which is inductive if the 
periodicity and depth of the corrugation are smell compared to the wavelength. 

The excitation of surface waves on dielectric-clad plane conducting surfaces has been 
discussed by Tai [15] for a line current source, and Brick [16] for a dipole source. Correspond- 
ing treatments for corrugated surfaces have been given by Cullen [17] for a line magnetic or 
infinite slot source, and Barlow and Fernando [18] for a vertical electric dipole source. In the 
case of dipole excitation, the field varies predominantly as the inverse square root of the dis- 
tance along both the dielectric-clad and the corrugated surfaces. Surface waves of this type, 
which have a phase velocity less than that of light, can be called “trapped surface waves” 
because they carry most of their energy within a small distance from the interface [19, 20]. 

It is the purpose of the present paper to study the general problem of a vertical dipole 
over a flat surface with a specified surface impedance, Z. A solution is obtained for the field 
in terms of the error function of complex argument. It is shown that, using an appropriate 
value for Z, the expressions for the Norton space and surface wave for a homogeneous conduct- 
ing ground are obtained as a special case [7]. By allowing Z to be purely imaginary and using 
asymptotic expansions for the error function, the field is shown to be of a trapped-surface- 
wave type [14, 19]. Another special case is for a two-layer conducting ground [21, 22] which 
can exhibit some of the features of trapped surface waves. Thus, a general connection between 
these various forms of seemingly unrelated surface waves is developed. 


2. Basic Formulation 


A vertical electric dipole of length ds and current / is located at a height / over a fiat 
ground plane which exhibits the property of surface impedance. Here the ratio of the tangen- 
tial electric and magnetic fields is specified on the surface of the ground plane and is denoted 
by Z. Choosing evlindrical coordinates, the dipole is located at (0,0,A) and the ground plane 
at z=0. As will be seen, the Hertz vector at (p,¢,2) need have only a 2 component, II, to 
satisfy the boundary condition; consequently, the fields are given by 


0"l] ; . Q? 
‘= dpa: E.=0, E.=(¢+—,)m 
at Rat 8-2 
tpw Op 


where uw=42r 1077 and B=27/free space wavelength. The resultant Hertz vector, Il, can now 


\ time factor exp (ivf) is suppressed throughout 
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be 











be represented as a sum of a primary and secondary part in the manner 


iI=TI,+Il,, (2) 
where 
- ie et, | 
i,=ce~ "1/9 with @ = end r,=|(z—h)*+ p?]’* (3) 
T We 


and II, is to be regarded as the influence of the ground plane. In the most elementary case 
of a perfectly conducting ground plane (i. e., Z=0), Il, would be given by the simple relation 
1,=ce~ 2/p, with re=|(2+h)?+ p*]*. (4) 

To find I, (or Il) for the more general case, it is necessary to express the primary field in 
terms of an integral representation [2] (sometimes called Sommerfeld’s integral), as follows: 


e~ *Pr ® Jo(Apje~ “2 hdd 


/ ( u 


where u ?— 8B and J)(Ap) is a Bessel function of order zero and argument dp Due to 
the presence of the branch point of u, the contour of the integral is indented upward slightly 
at\=6. Physically, this integral can be regarded as a means of representing the primary field 
as a spectrum of plane waves whose angles of incidence, 6, on the ground plane are related to 
the integration variable by \=8 sin 6. Consequently, the angles include both real and com- 
plex values, ranging from 0 to 2/2 along the real axis and from 2/2 to 2/2+7 along a line 
parallel to the imaginary axis 
The form of eq (5) suggests writing II, in the following form 


Jo(dAp)e~ "Et" PR(d)Add 


ul 


Il, (6) 


where (A) is now a reflection coefficient for plane wave incident on the ground plane at angle 


of incidence @=sin~'d/8. Consequently, 


, u—ipaA 
Rr a iBA’ i) 
where 
+ ie 
A ° 
p® 1207 


In the limiting case of a perfect conductor (Z=0) it can be seen that R(A)=1, and the integral 
is of the form of a Sommerfeld integral and TI, is again given by eq (4). 
The Hertzian potential II can then be written in the form 


pry pts 
uz‘ rh ee (8) 
ry V2 


where 2P is the effect of the imperfection of the ground plane which results from a finite value 


of the surface impedance 7 
3. Evaluation of the Integral P 
The integral 


iBAA é 
W- IBA)u 
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~ Jo(dp)dd (9) 











must now be evaluated. Before proceeding, certain conditions and restrictions must be clearly 
stated. At a later stage, these may be relaxed somewhat. For the moment 


l\Al<1, O< arg A 7 and Bp > >1. 


Since 
J (x) - | #13 (7)+H? (| 


HH} (a H”(—2z) 


and 


where 77)” and //” are Hankel functions of the first and second kind, respectively, it follows 
that 


2 (u+iBAyu 


BAL” he~*et 
P ad | H§? (Ap)dr LO 


Introducing the substitutions \=8 cos a and A=sin a, it is seen that 


p iBA (" " cos ae~ +h) sina HH? (Bp COSQa@) d 
a c 


{ a, Ll] 


x, a—a 
Sin > cos 


. ( 


Essentially P is now a spectrum of plane waves traveling away from the ground plane so that 


the vertical component of the propagation constant is 8 sin a. Again, complex angles are 
included in the spectrum. When |8p cos a! 1, the first term of the asymptotic expansion 


of the Hankel function can be employed. That is, 


s ? 
Hi Bp cos @)~ — Gee ag Ors, 12 
7TBp COS @ 


When this is inserted into eq (11), the equation for P becomes 


A re} 2 , COS a)2¢- BR cos (O—a 
P=—e nS (FYI do 13) 
2 2rp zs ; a-+-Ay, a a 
Sin = COS a 
where AR and @ are defined by o+-h RP sin @and p PR cos 86. In the usual case, the separation 


p between transmitter and receiver is large compared to their respective heights, h and 2, and 
consequently, @ is slightly less than 7. It can be seen that the exponential factor in the inte- 
grand of P is rapidly varying except for a region near a=@~z7. This is, of course, the saddle 
point of the integrand or the point of stationary phase. The important part of the integrand 
is in the region near the saddle point. In fact, this is the justification for employing the 
asymptotic expansion for the Hankel functions //)"; the argument 8p cos a@ is always large 
in the region near the saddle point if 8p itself is large and 2+A is less than p. 

The integral for P is now in a form where the saddle-point method of integration can 
be applied. The usual technique [23] is to deform the contour (which in the present case is 
along the negative imaginary axis, the real axis from 0 to x, and then along a line parallel to 
the positive imaginary axis) to a path of steepest descent. For example, if we let 


cos (a—Fé l ir 
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an 








and let 7 range from > to © through real values, an integral of the type 


ee 
G(cos aje~ BROS \a—® day (14) 
Nasttis 


where G@ (cos a) is slowly varying at a=9, is transformed to 


v4 Gi ‘ 
as\1 » 7m COS @) > 22 ” 
(27)2%e—t8R : e~ Rd, (15) 
» yl—tr?/2 


In this deformation of the contour, account must be made for the singularities of the integrand 
that are crossed. In the case of poles, 27i times the sum of the residues is added to the new 
integral. The final step in this classical procedure is to expand G(cos a)/y1—ir?/2 in a power 
series in + enabling the integration to be carried out term by term. The leading term of this 
resulting asymptotic expansion is 


(292)’2G(cos 0)e~ ®*®(BR)-* (16) 
and succeeding terms contain (6R)~”, (BR)~”, and so on. 


As is often the case, the contribution from the saddle point (which is the branch point in 
the original \ plane) cannot be separated from the pole(s) of the integrand. For example, in 
the present problem, 


1 2 cos 
(cos a) %¢ ipR cos(@—a) 


G (cos a (17) 
. ata a—a 
Sin > cos 


has a pole at a=2+ a which is near the saddle point a=@ since a is small and @ is near x. In 


other words, the integrand is not slowly varying near the saddle point on account of the factor 


cos (a—a,)/2. The integral is then of the form 


7 +i isk cos (4 a 
I da (18) 


after the slowly varying factors have been separated out. The necessary modification of the 
saddle-point method to treat integrals of this type was devised by Van der Waerden [24] and 
Clemmow [25]. On making the usual deformation of the contour via the substitution cos 
(A— a | —7r*, the integral now becomes 


] .-te/4.-@2./9 9—ay\ (** ePRrdr 
i af “¢ a y= COS ( ; ) 
a 


2 fey 0—a \ (19) 
tT*+12 cos ( = ) , 


It can readily be verified that no poles of the integrand are crossed in this deformation for the 
argument of a (or A) less than 45°. The integral to consider is now of the type 


dr (20) 
oT 


« ré 
e=12 | cos ( : — )]: 


Since | et dr (*) it follows that 
J 


where, for real r..2=—BR and 








Now it can be seen readily that 


d M2 
“ie “4=-() = 


and an integration with respect to £ from s to © leads to 
or, after a change of variable 


which holds for a// ¢ in the e plane with a cut along the negative real axis. This can be written 


A e“erley re, 25 
ve 


where erfe is the complement of the error integral given by 


2 
ET TCZ 6 é dz, 26 
V7. 
where the integration in the z plane is directed from 2) towards the right to Therefore 
I 12 wre IR o—Corfe(iw' : 27 


where 


w Ble 2B | os ( : — | , 


and finally, since R=r,, 


j iBAA 
P : ( . . Jo(Xp dx 28 
: Up t+i1BAa 
~) mp)’ erte rw ( 7 Mo. 
where 
oT h 
w=p\ i+ 29 
P| Ar, ) 
and 
iB; 
D : - J pe 


While the derivation was carried out for the argument of A lying in the range 0 to 7/4 or b 
in the range 0 to —x/2, the formula is actually valid for all values of 6 if the above definition 
of the error function complement is used. The justification for this step is based on the principle 
of analytical continuation. If, in the initial derivation, the argument of A was allowed to exceed 
7/4, we would find that a pole would be crossed in the deformation of the contour in the a 
plane. There would be a corresponding change in the integral A but, after a change of the sign 
of the variable 7, the expression vielded for J would be unchanged in form if we are consistent 
in the definition of the complex error function 


4. General Properties of the Solution 


The veneral behavior of the result will be considered now, For the moment the heights of 
the transmitter and receiver are assumed to be very small. That is (2+h)<<<p and, con- 
sequently, w~p and r,~r,~p. The Hertz potential then has the form 











where 


F(p)=1- i(rp)'*e-"erfel(ip (31) 


The vertical electric field is then given by 


, 1a Oli 
iD p (32) 
p Op Op 
and, since Bp > >1 and F\(:) is a slowly varying function, it follows that 
E. +- 2¢B*e- 1B p 'F(p). (33) 
For a highly conducting flat ground plane, [!p|<<1], the electric field denoted by E® is 
given by 
Ky = 2¢B*e~ ep. (34) 


This result can be obtained by elementary methods, being just the free space radiation field of 
an electric dipole of moment 2¢. The factor 2 arises from the imaging of the primary dipole in 
the perfectly reflecting surface. The ratio 


F(p)=E,/E (35) 
is by definition the complex attenuation factor and is a measure of the efficiency of trans- 
mission of waves along the interface for a vertical-electric-dipole source. 


The behavior of F(p) is verv dependent on the phase of p, and since 


p—|per——i8? y2__ ibe (ZY. 


2 2 \1200 (36) 


the surface impedance, 7, plays an important role. For example, in the case of a homogeneous 


ground of conductivity o« and dielectric constant e, 


| a - 
1200 n | l—x: ) (37) 


V k io °] 


The refractive index, N, is usually large compared to unity with a phase near, but always less 


with 


than —45°. Correspondingly, for a homogeneous half-space, the phase b lies between 0° and 

90°. This is precisely the case that has been treated in detail by Norton, who gives numerical 
results of /'(p) as a function of |p) for 6 varying from 0 to —180°.* In all eases for |p|>0, |F(p) 
was less than unity. In other words, the field /, was never greater than £> for the perfectly 
conducting ground plane. For purposes of numerical work, the error-function complement 
can be expanded in an ascending power series, as is well known. This leads to the following 


series formula for F'(p 


: 2p (2p 
I p | uUmp)e-? 2p P - E eee (38) 
1-3 1-3-5 
which is absolutely convergent 
The second term of this series, the quantity —7(mp)'*e~’, is rather interesting if it is consid- 


ered by itself. The vertical electric field, £3 , which would result by ignoring the other terms, is 


De 


Ds 2¢p- i(rp?*)e~? (39) 
p 
") Z 
2T ail 
Ce wb ’ Ae | . | ’. (40) 
\ Bp 
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This field has the characteristics of a Zenneck surface wave; it varies inversely with the distance 
along with an exponential attenuation, and it has a phase velocity greater than that of light. 
By retracing the steps in the derivation it is not difficult to show that /} is the residue at the 
pole ap in the complex a plane. Taken by itself, £} could easily exceed EY by a great amount 
if 6 was near —90° corresponding to the phase of Z (or A) being near zero. This corresponds 
to a small conductivity in the lower half space. 

The above argument would indicate that the Zenneck surface wave is excited by a vertical 
dipole source on a homogeneous half-space. Of course, when the other terms in the series for 
F(p) are considered, the total field does not have the attenuation characteristics of a Zenneck 
surface wave. For example, if the asymptotic expansion for the complement of the error fune- 
tion is employed, the following series for F'(p) results: 

| 1-3 1-3-5) 1-3-5-7 
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Fi Pp = 


which are also asymptotic and valid for large values of p and —2x%<b<.0. This would indicate 
that the field /, at large ranges varies inversely as the square of distance, whereas the Zenneck 
surface wave varies inversely as the square root of the distance. This viewpoint is confirmed 
by the work of Wise [26] and Rice [27], who develop rigorous expansions for the Hertz potential 
on a homogeneous earth. 

The absence of any features associated with a trapped surface wave on a bomogeneous 
earth is due to the restricted range of b (namely between 0 and — 7/2) for a homogeneous earth. 
It is of interest to note that Sommerfeld in his 1909 paper had an error in sign which was finally 
pointed out by Norton [3]. This error is equivalent to the interchanging of the 6 and —6. In 
other words, the complex conjugate of p is used in place of the proper quantity. Rolf [4] 
employed Sommerfeld’s formulas along with the incorrect sign. Rolf’s series expansions for 
F(p) would correspond to eq (38) with 6 replaced by —6. His calculations therefore correspond 
to values of /(p) for 6 from 0 to +90°. The appearance of his curves is indeed weird. As 
will be indicated below, the function F(p) can exceed unity in certain instances when 6 is posi- 
tive. This is indicated in Rolf’s curves. The irregularities and dips to zero, however, must 
be attributed to errors in computation. The series formula for /(p), although mathematically 
convergent, takes a great num ber of terms when p is greater than 5 or so. And, furthermore, 
the asymptotic expansion of Fp) is not accurate for the real part of p less than about 15. The 
anomalous features in Rolf’s curves occur just where the accurate computation of F(p) would 
require tens or hundreds of terms in the power series formula. This viewpoint is confirmed 
when more elaborate calculations are carried out [21]. 

The nature of the attenuation function /(p) when 6 is positive can be best seen from an 
asymptotic development. This is obtained by noting that 


wr 2 2 
erie ip é dr 2— ( dae 12 
\ We \ Te Pp 
The latter integral asymptotically vanishes for p-> if the real part of jp” is less than zero 
Therefore 
vr ie*! 1-3) 1-3-5 
erfelipyea— 2" (1414134195, ) 3 
\ 7p <P \<pr \<p) 
for 2x >b>0 and p >>1. Utilizing this result, it then follows that 
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in the asymptotic sense. The term —2i, xp e~” has all the characteristics of a surface wave 
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term almost vanishes compared to the term 1/2p. 








It is not present in the asymptotic development for /(p) when 6 is negative. Of course, for 


dissipative media such as a stratified conductor, the real part of p is large and the surface-wave 


The most interesting case is when the surface impedance is purely inductive. For example, 
if a flat perfectly conducting ground plane is coated with a thin layer (of thickness h) of pure 


dielectric (with refractive index N), then the surface impedance, Z, is given to a good approxima- 


Z~ ipwh | xs} (45) 


Another example is when a flat metal surface is corrugated in the form of periodic rectangular 
slots of width d and depth DV with a separation w between slot centers. The surface impedance 


for the magnetic vector parallel to the directions of the corrugations is then 


tion by [19] 


— , 
Z=1 ~~ No tan BD, (46) 
o- 
provided that there are at least five corrugations per wavelength along the surface [13]. For 
D< y/4, the surface impedance is purely inductive, and since 
iBp,- 7 
p z (Z/m)’, (47) 
p+ip and 6=7/2. The attenuation function then has the asymptotic development, 
7 un. , oo 157 
I Pp) — Je el pl ae rm te ee | (48) 
<\p tp) Sp 
where 
Bp xX? - , 
Pp -and Y¥=|Z 
2 % 
The vertical electric field at 0 is then given by 
FE’, =~ 2cB*e~ pF P). (49) 
The surface wave term now veries as 
l Pt 
e~i8( 1+55 po 
\ Bp 
which 7s not a Zenneck surface wave since the phase velocity, D, given by 
l l » 
’ (50) 


is less than the velocity of light, 1 \ Mo €- 


If the surface is purely capacitive such that Z iX’, the argument of the error function is 


ip’ ¢ iz ‘Bp 2 a nm)”, (51) 
so that erfe(ip’) can be expanded directly to lead to the asymptotic representation 
” 3 152 > 
cee (52) 


Fi p)=- 


which contains no surface-wave term. 
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5. Height-Gain Functions 


It is desirable now to return to the general formula for F(p) which is a (modified) saddle 
point approximation for the field of a vertical dipole over the impedive surface. For h=0, it is 


W=1 iy 1 pe “erfe(iw'), 53 
where 
eo 45, (14) 
For small heights and receiver it is of interest to develop Win a power series in z. For example, 
a , (dW 2 (fw \ : 
i Wi0)+2( 7 ) ta([5) +--: 54 
Now 
dW _p tetetteane dw _ 
ds - iy wpe ericviw? d- oo 
nnd 


dw ; : ] 
d: 2p ( - Ap Ap 


| | ik AW (0) =ikAF( p 


al ~} p l- tk A>). 56 


Therefore 


and, to a first order in 


The generation of higher-order terms using this method is cumbersome, and it is really better 


to return to the original integral for the Hertz potential II(z) given by 
a es Xo - 
I I Fre “dy 57 
u- IBA 


for h=0. Expanding the exponential, this is rewritten 
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Now since II satisfies the wave equation, 0°/02° can be replaced by 
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the contribution from the second-order term is 


22 1d 9d Qee- ise 2? i(Bz)? 2? F’(p) ...,F’(p) , 22 F” 
| 4 p | Ho o Fis) 4 5 oP) igge © Se 
2 p Op Op <p p- (Bp) 2p Fi —P) Fi —P) 2 Fi Pp) 
where the prime(s) indicate a derivative(s) with respect to p. It is not difficult to see that 
for 2x —b—0, F(p) is slowly varying compared to other factors and the term in curly brackets 


is much less than unity if Bp l and p* ~~. In the case of the inductive surface (i. e., 


In >h>0O), it is convenient to write 


F(p)=F (p)+ Fi (p), (62) 
where asymptotically 
| P)> 2iy mp ¢ : 05) 
is the surface wave term and 
; l 3 15 
F.(p) = —— — Cer (64) 
<p 4p Sp 
is the remaining nonsurface) contribution \For the case Ia< b- 0), F(p)=F,(p)| 
Now 
F'(p ] ] 
AP) +0 ( ) 
FE. P p p 
and 
Fp) : 
NE +0 ( ) 
} i p ) r 
f / 
but 
Fp) ipa 
F’.(p 2 
ana 
F* P) ( iB A” ) 
F’.(p) 2 
Therefore, the Hertz potential II,(2), corresponding to the surface-wave field, is given by 
; (7BAz)’ 
H(z) 10 I +iBAz+ — a *? | 65) 
~ ITO} 
whereas 
1,(2)~11,(0)[1+7BAz}. (66) 


In the case of a purely inductive surface (Z=iX), the exponential factor for the surface 


term is given by 
XYS8/n éewr (67) 


indicating a rapid decay of the field with height. 
By a slight extension of the above, the case of both transmitter and receiver raised yields 
Il. ~h)~TI.(0.0)¢ BACe+h) (68) 


and 


II, (2,h) ~ 11, (0,0) (1 + 7BAz) (1+ 7Bdh). (69) 

















Still another approach to the effect of raised antennas, which is particularly appropriate 


for large heights, is now considered. In general, 


where 


with F(w)= F,(w)+eF,(w), where e=0 for —24r<b<0, e=1 for 2x >b >0, and w= p(1+ C/A)?, 
CU=(2+h)/r.. Asymptotically, 





F (w) QiVarw ¢ - 72) 
and 
F(w) -- -- ra _ Lo — ee ee 73 | 
2w 4u- sw’ 
After some algebraic manipulations, 
Ii=II,+- ell 
TT, e 1 C—A e~ #2 l 1-3 1-3-5 ” lia = 
=f 4+.— 7 . — oe ne 8 i4 
c ry C+A fro pui+cya 2p°1+C/A)y’ 4p*(1+C/A) re 
and 
IT 24 1,- —_ _ 
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The first term in the expression for II, is, of course, the direct or primary influence. The 


second term is the secondary influence modified by a Fresnel reflection coefficient. The terms | 
containing inverse powers of p can be considered as corrections to reometrical optics. The 
surface-wave term is absent when 6 is negative. When 3 is positive, it is a finite contribution | 
to the total field I, and for a highly inductive surface (i. e., 5=90°) it is predominant. | 
The vertical electric field is simply given by 
B= 2S 81 —C?) 0 for Br, and 8B i. 

| 

and furthermore, 7,~r,—2hC. 
6. Conclusion 


It has been shown that the nature of the field radiated from a dipole over a flat surface 
depends markedly on the complex value of the surface impedance Z. For example, when the 
phase of Z is between 0 and 45°, the derived expressions for the field correspond to those of 
Norton and others for a dipole over a homogeneous flat ground. When the phase of Z exceeds 
45° the radiated field becomes trapped by the surface. This effect is most pronounced when 
Z is purely imaginary, corresponding to an inductive surface. A similar unifying treatment 
for a dipole over a curved spherical surface will be given in a sequel to this paper. 


The author thanks Kenneth A. Norton for his very helpful suggestions and Howaré 


Bussey for his constructive criticism. 
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1. Introduction 


Most of the commercially available permanent- 
| magnet materials of high coercivity are mechanically 
hard and brittle and must be cast or sintered into the 
required shape for use. A magnet alloy of approxi- 
mately 60 percent of copper, 20 percent of nickel, 
and 20 percent of iron (called Cunife or Cunife I by 
its producers in the United States) has an unusual 
combination of magnetic and mechanical properties. 
Although it has a coercivity of several hundred 
oersteds, Cunife is so ductile that it is possible to 
eold-reduce it considerably. 

Several investigators [1,2,3] ' have reported on the 
magnetic properties of Cunife wire and on the effects 


of cold-working and subsequent heat treatment. 
| Their results have shown that even if the material ts 
ecold-worked to the point at which the magnetic 
properties are adversely affected, the initial magnetic 
properties may be recovered or improved by a simple 
‘at treatment or baking. These reports, however, 
pertained only to the larger wire sizes. No data 
vere available on the magnetic properties of Cunife 
lrawn to wires a few thousandths of an inch in 
diameter 
The solution of a special magnetic problem re- 
quired evlindrical magnets approximately 0.50 in. 


long and 0.005 in. in diameter2 The ductility and 
high coercivity of Cunife made it appear suitable for 
this application, and the present investigation was 
undertaken in order to evaluate it 


2. Material 


The commercial material used in this investigation 
had a composition of 19.3 percent of nickel, 19.5 
percent of iron, and 60.4 percent of copper, and was 
supplied by the maker in the form of wire with a 
diameter of 0.025 in. and in the heat-treated condi- 
tion. This wire was cold-reduced from 0.025 to 0.020 
in. through tungsten carbide dies with five passes. 
Diamond dies were used for the further reduction 
from 0.020 to 0.005 in. with 21 passes. Potassium 
soap was used as the lubricant for all passes, and no 
heat treatment was given the wire at any time during 
the reduction from 0.025 to 0.005 in. At various 


—_—_——————____. 


it the end of this paper 


Figures in brackets indicate the literature references 
the Mine Fuze Branch 


This problem was related to a classified project of 
Diamond Ordnance Fuze Laboratory 


having a com position of 60 percent 
20 percent of nickel, and 20 percent of iron, was cold-drawn from wire, 25 mils in diameter 
Demagnetization curves were obtained for several sizes of the wire 
These wires were then given a heat treatment and again demagnet- 
Values of coercivity, 
product are compared for various wire diameters. 
they 
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of 


copper, 


retentivity, and maximum energy 
Although the magnetic properties were 


were substantially improved by the heat 


stages of reduction, samples were cut from the coil 
of wire for later measurement. 


3. Heat Treatment 


The magnetic properties of the cold-drawn Cunife 
wire may be markedly improved by properly heat- 
treating the wire. Since the maximum temperature 
at which the wire is heated affects the final magnetic 
properties, several heat-treating temperatures were 
tried. The resulting magnetic properties for samples 
of 5-mil wire when heat-treated for 1 hr at various 
oven temperatures are shown in figure 1. During 
this heat-treating process the wire was embedded in 
iron filings to prevent oxidation. In each case the 
temperature of a muffle-tvype oven was allowed to 
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Ficure 1. Influence of temperature of heat treating for 1 hr 


on coercivity and retentivity of 0.005-in.-diameter Cunife wire. 
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reach the heat-treating temperature before the wire 
and iron filings were placed inside. The wire and 
filings remained at the heat-treating temperature for 
| hr and then the heaters were turned off, the oven 
door opened, and the oven allowed to cool to 300° C 
This cooling usually required 20 to 30 min. The 
sample was then removed from the oven and cooled 
to room temperature. 

The results in figure 1 show a definite maximum 
value for the coercive force of Cunife when heat- 
treated about 620° C. This is in accord with the 
results obtained by Geisler [3] with the allov Cunico. 
He has shown an interesting curve for the change in 
coercive force of Cunico as a function of aging time 
when held at various aging temperatures. The 
maximum coercive force for the three isothermals 
occurs at different aging time; if the proper aging 
time is selected the coercive force will increase, go 
through a maximum, and decrease when plotted as 
a function of aging temperature. 

In order to determine the effect of oxidation on 
the magnetic properties, a specimen of 0.005-in 
wire was placed in an oven at 610° C in an atmosphere 
of air and was heated for 1 hr. The wire was then 
slowly cooled to room temperature. The surface of 
this wire was then gravish black. However, the 
magnetic properties were the same within 5 percent 
as those of wire of the same size, oxidation of which 
had been prevented by heating in iron filings 

It might also be expected that, since wire specimens 
of different sizes have had different amounts of 
cold-working, different heat treatment might be 
required to obtain the best magnetic properties of 
wire of each size. However, using 0.020-in. (36% 
cold -reduced) and 0.005-in. wire (96% cold-reduced 
as test samples, no such effects were found. The 
maximum heat-treating temperature of 610° C and 
| hr of heating time gave the optimum results for 
wire of both sizes 

As a result of these experiments, all wires were 
heat-treated in the same manner. The procedure 
was to heat the wire specimens, while 
embedded in iron filings, to a maximum temperature 
of 610° C, to hold them at this temperature for 1 hr, 
and then to allow them to cool slowly to room tem- 
perature. This produced wire that was bright and 
showed no noticeable signs of surface oxidation 


selected 


4. Apparatus and Measurements 


The normal induction and demagnetization curves 
of the Cunife wires were obtained by using a High-H 
permeameter [4] and a calibrated ballistic galvanom 
eter. The magnetizing force was measured with a 
flip coil, and the flux density in the wires was meas- 
ured with a coil surrounding the specimen. Air flux 
was compensated with an auxiliary coil. The forms 
for the B coil and compensator coil were of aluminum, 
each form being 0.75 in. long with an inside diameter 
. of approximately 0.1 in. and a wall thickness of 
0.010 in. Each form was wound with 15,370 
turns of AWG 44 varnish-coated copper wire. The 
Cunife specimens were cut 0.75 in. long, and a suffi- 
erent number of wires 


coil 


were used to have a cross- 
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sectional area of at least 0.00040 sq in 
netization curves were obtained from 
magnetizing force of 3,000 oersteds 

) 25-1n.-4 


5. Results 
The shi 


The dependency of the coercive force on the redue.|aved con 
tion of area of the wire is shown in figure 2. Sine This can 
the 0.025-in. wire was in the annealed condition |idea that 
initially, the percentage reduction in area is also the jepends 
percentage cold-working for the lower curve. Ag a] of ferrom 
result of the cold reduction of the wire from 0.025 to} yre of se¢ 
0.005 in., the coercive force decreased from 550 to 501 The re 
Subsequent heat treatment restored the | product, 
initial coercive force until the wire was reduced jy }iy figure 
area more than 85 percent, or to a diameter of approx. | decrease 
imately 10 mils. Further reduction in area caused | yptil a! 
a sharp decrease in the values of the coercive foree | peached. 
/7,; for the heat-treated 0.005-in. wire, /7. was only | ereases 
60 percent of the value for the 0.025-in. wire. 

This sharp decrease in coercive force ut the smallest 


7 he demag ver the 
au Maximun * ated, th 
re? 4 


i nereased 


oersteds. 


show a 


approac 


sizes of wires studied is of special interest.  Sinee } y rapid 
the major role in production of coercive force in The l 
Cunife is believed to be due to strain disregistry , js show 


between a precipitate and the matrix, it is quite ipper ¢ 








likely that for the smaller sizes of wire the severe | of the 4 
cold-working has so distorted the grains that some | Jater, t 
of the effects of distortion may persist even after 
the wire has been aged. As shown in figure 2. th 
coercive force of the cold-drawn wire continually 
decreases as the cold-working increases. However, | 
aging will more than restore the initial coercive force : 
because of the preferred orientations as a result of ¥ 
cold working [5, 6, 7]. If the cold-working exceeds 
the optimum amount for the small wires, then the 
initial coercive force of the larger wire size cannot be 
regained, indicating that the awnisotropy of the alloy 
was changed. - 
Figure 3 shows that the value of the residual in- 
duction, 2,, remained approximately constant until 
the reduction in area exceeded 80 percent, then de- 
creased sharply for further increases in cold reduction FIGURE 
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ifter the successive samples of wire were heat- 
IMup 


ated, the residual induction, B,, for each wire size, 
| opegsed to or exceeded the initial value for the 





925-in.-diameter wire. 

The shape of the residual induction curve for the 
"edue. | sved condition shows no dependence on wire size. 
Sing this can be explained on the basis of the accepted 
dito, dea that the value of residual induction for alloys 
SO the lepends primarily on composition and the amount 
As al of ferromagnetic phase, while strain and structure 


125 to| are of secondary Importance. 





toii} The relationship between the maximum energy 
d the! product. (B,f/Ze),, and the reduction in area is shown 
ed ip} figure 4. For the cold-worked condition, the 
prox-| Jecrease in the maximum energy product is small 
Nused | yntil a 5O-percent reduction in area of the wire ts 
fore | reached. Kor ereater cold-working, BH ’ de- 
only | creases quite rapidly. The heat-treated materials 


show a gradual increase until the reduction in area 
pproaches SU percent, and then the curve dips and 


illest 


Sine , rapid decrease 1D Bl occurs 
ey The information contained in figures 2, 3, and 4 
istry is shown in figure 5 in a different manner. The 
uite ipper curve is for the ratio (2B .),,as a function 
vere of the diameter of the Cunife wire Here. as well as 
Ome later, the subseript cw designates the cold-worked 
iftey 
the 
ally 
ver } 
res A 6 B 
of -“ 
eds | 
thi 
he 
i 
In- 
tl 8 
fe | 
n_ | Ficure 3. Retentivity as affected by cold reduction and hee 
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condition and Aé designates the heat-treated condi- 


tion. 


This curve shows that the value of B, is 


practically independent of the cold-working until 
the wire has been reduced to about 8 mils. 


This figure also contains the data for the ratios 
{7A 


tions of wire diameter. 


(F.) ns and (Balla) ml cu (Balla mnt as func- 


The points for both sets of 


data fall along the same curve and both indicate 
the improvement in magnetic properties that results 
from the proper heat treatment for the cold-worked 
wire 

The demagnetization curves for heat-treated Cunife 
and 0.020 in. 


wires of diameters 0.005. 0.012, 


shown in figure 6, 
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It is possible to improve the magnetic properties 
of 0.005-in.-diameter wire by starting the drawing 
process from a heat-treated wire having a smaller 
diameter than 0.025 in. This of course will require 
less cold-working in reducing the wire to 0.005 in 
For example, a 0.012-in.-diameter wire, in the heat 
treated condition, was cold-drawn 0.005 in 
Wire samples of diameter 0.010, 0.008, and 0.005 in 
were heat-treated and measurements made of thei 


to 


magnetic properties The results obtained for the 
coercive force, /7,, and residual induction, B, from 
wires drawn from 0.012 and 0.025 in. are given in 
the following table: 
Hl R 
Wire di- 
ameter Initial Initial Initial Initial 
wire di wire di wire di- wire d 
imeter, ameter, ameter aie 
0.012 in 0.025 in 0.012 in 0.025 
’ Oerste Oersteds KN ilogqausses A d 
0.010 625 580 16 1s 
OOS 610 190) L 6 18 
OOS 2) 360 L 6 a. 


Figcure ri shows the relative otte ets oft cold workine 
and heat treatment on the small sizes of Cunife wire 


when the starting diameters of the heat-treated 
wires are 0.025 and 0.012 in 
6. Conclusions 
Permanent magnets as small as 0.005 in. in diam- 


cold-drawn from commercially available 
Cunife wire If the cold-drawn wire is subjected to 
i simple heat treatment, the permanent-magnet 
properties of the material are substantially IM pros ed 
The magnetic properties for the smallest wire sizes 
resulting from the final heat treatment, are 
parently affected not only by the diameter of 
wire but also by the cold-working 


eter may be 


up 


thie 
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tra of the very complex mixtures of substances 
the temperature of liquid oxygen is described 


chemical compounds in such a mixture is greatly in excess of the 


number of masses reeorded, and it is accordingly impossible to identify all the compounds 


present However, the 
products from each source are of substantial 
a representative pattern of each and to proce 
the character and amount of pollutants in tet 


nd is made by treating the gases associated with 


s 


Introduction 


At an early stage in the study of air pollution in the 

Los (ngeles area, it was proposed by Martin Shep- 
erd 1] that condensible ASCs he collected from the 
umosphere by cooling and filtering at the tempera- 
ture of boiling oxvgen and subsequently subjected to 
analysis by spectrometer. During 1949-51, 
approximately 50 samples were taken in the vicinity 
of Los Angeles by this method, all of which gave very 

/ complicated spectra whose Interpretation was under- 
taken with only limited suecess 2. a Subsequently, 
similar samples were taken 1) several under 
various circumstances All of these have shown ra 
Dasic similarity, Which indicates that most of the 
polluting gases in each urban atmosphere are the 
same chemical in much the same _ pro- 
portions, and are of the general character of “products 


mass 


cities 


substances 

ite | of combustion 

| It is fairly obvious that in may expect 
large volumes of products of combustion to enter the 
air from three sources: (1) From fuel used for heating 
buildings and producing power, (2 from automobiles. 
and (3) from waste disposal. In most cities the first 
simportant only during cold weather and the third 
sintermittent or localized The 
ous and universal It would be desirable and may 
vet be possibli to find distinguishing characteristics 

but their principal 


ith \ city we 


second is continu- 


between these three sources, 


differences appeal 1o depend on the completeness of 


combustion of the fuel rather than on its type. 
whether coal, petroleum products, or wood and paper. 
When combustion is “complete,” only carbon 
dioxide, water vapor, and a little sulfur dioxide, 


mitre oxide, and nitrogen dioxide, are added to the 
ar, and oxvgen, in negligible amount, is removed 
irom it When combustion is “incomplete,” but as 
nearly complete as it can readily be made 
the conditions of burning, carbon monoxide is the 
chief component of the “products of incomplete 
combustion,” and hydrogen and various hydro- 
arbons are the most easily identified of the addi- 


“hevy Chase 15, 


fet 


Md 





significant sources of pollutants are much less numerous, 


under 


and if the 
ly constant composition it is possible to obtain 
ed thereafter to make a fairly good estimate of 
ms of their sources. The greatest simplification 
motor traffie as a unit 

tional products. As conditions are changed to make 
combustion less complete, products become more 
complex and more dependent on the character of the 
original fuel; but to the stage in which visible smoke 
appears, sufficient to attract attention and to induce 
remedial action, only moderate differences of compo- 
sition occur, 

The temperature at which samples can be “frozen 
out” without condensing oxygen from the air at the 
same time is too high to permit the collection of 
hydrogen or carbon monoxide, or methane unless it is 
present in explosive concentration. The amounts of 
other gases that escape condensation depend mainly 
on their at the boiling point of 
OxXVeen. 

Most of the air in the cold trap after the initial 
sampling is removed by partial evacuation at the 
boiling point of nitrogen, and treatment is_ then 
given the sample to eliminate some of the carbon 
dioxide and water. The sample is then warmed, 
usually only to room temperature, and the vapors 
are admitted to the evacuated inlet reservoir of the 
mass spectrometer. In this procedure another por- 
tion of the pollutant may be lost, i. e., that which 
condensed from the gas phase but did not re-evapo- 
rate at the time of transfer to the spectrometer. 
This loss includes a portion, but not all, of each 
high-boiling compound that occurred in the initial 
atmosphere to or nearly to the point of saturation. 
The losses that occur both through failure to condense 
in the first place and to re-evaporate later are dis- 
cussed in the next section, 

The mass-spectrometer a “spectrum” 
produced by the amplification of ion currents each 
of which is proportional to the total charge of ions 
having a single value of W/e, where M is the sum of 
the atomic weights of the atoms in the ion and é¢ Is 
the positive charge carried by it (the number of 
electrons removed in the ionization process). The 
charge is usually unity, and its measurable indication 
on the record is referred to as a “peak” designated 
by a “mass number,” which is the weight of the ion 
with unit charge that would produce the peak. In 
this paper, the mass number enclosed in parentheses 
will indicate the ion current represented by a particu- 


vapor pressures 


records 
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lar peak, usually expressed as a percentage of the 
total ion current of all peaks of interest. 

Chemical compounds produce ions either by the 
loss of an electron, leaving a “‘parent’’ peak that is 
only undetectably different in mass from the parent 
molecule, or by the loss of one or more atoms from 
the parent molecule, leaving a charged fragment. 
Each compound breaks up in a characteristic manner, 
producing a “pattern” by which a compound that 
occurs alone can almost always be identified with 
certainty. However, the ion currents with the same 
mass-to-charge ratios from a mixture of different 
compounds are additive. In theory at least, we 
could determine every compound present by the 
solution of an equal number of simultaneous linear 
equations if we had as many or more peaks on which 
to base the equations and patterns of all compounds 
involved. 

All of the spectra of urban pollutants so far 
obtained have contained 80 or more usable peaks; 
but in every case there has been evidence that these 
peaks represented many more than 80 compounds, 
and their complete identification was obviously 
impossible. Even when a sample was separated by 
multiple distillations into 50 fractions, the simplest 
fraction still contained traces of more compounds 
than can be identified with certainty. 

However, if we consider significant of 
pollutants, it is probable that their number is much 
less than that of the usable peaks, and if the products 
from each source are of substantially constant com- 
position, we may be able to obtain a characteristic 
pattern of each and to proceed thereafter to make a 
fairly good estimate of the character and amount of 
the pollutants in terms of their sources A procedure 
for making such estimates and some of the results 
obtained is the primary subject of this paper; but 
the main theme will be preceded by a discussion of 
sampling and laboratory procedure intended to make 
the analytical method understandable and essentially 
complete. 


sources 


2. Character of Sample Submitted to Mass 
Spectrometry 


As indicated in the introduction, the sumple that 
enters the mass spectrometer contains everything 
that was in the gas phase in the atmosphere except 
1) vapors below their concentration of saturation 
at the boiling point of oxygen, and (2) vapors 
condensed at that temperature but not re-evaporated 
at the time of introduction into the spectrometer 
The second group of losses is discussed first. 

The cold trap (Shepherd) that was generally used 
has a volume of about 300 ml. Typically about 100 
liters of air are passed through the trap. If the 
atmosphere is initially saturated with a given com- 
pound, its condensate, when warmed to the initial 
, temperature for introduction into the spectrometer, 
will be confined to about 0.3 percent of the volume 
occupied before sampling, and consequently only 
about 0.3 percent of the compound collected will 


immediately re-enter the Vapor phase. However. 


this vapor is introduced into an evacuated al 

relatively large inlet reservoir of the spectromete 
and as pressure in the saturator decreases during the| 
transfer, more of the condensate is evaporated. l 
the case most favorable for the detection of a hig} 
boiling compound, the freezeout tube contains op! 
as much total volatile material as will produce 

pressure satisfactory for analysis in the combing, 
volumes of tube and inlet reservoir. In this cag 
the compound should be detectable if the spectro. 
meter is only as sensitive as would be necessary 4 
show the presence of the compound in solid or liquii 
form in the reservoir itself. This is the Conditioy 
roughly represented on figure 1 by the notatio, 
“approximate limit of detection.”” Actually thy 


freezeout tube will usually contain enough material | 


for several spectrometer runs; the inlet  reseryoj 
will be filled from it through a passage small enoug! 
to prevent much mixing between the reservoirs, and 


the composition of the gas entering the inlet reservoir | 


will not change greatly during the transfer. 
tivity for the compound is then reduced by a facto, 
roughly equal to the number of times the inlet 
reservoir can be filled from the contents of th 
sampling tube. 

Generally a compound should be “detectable” if i 
occurs in the atmosphere to the extent of 5 percent of 
saturation and if, at the temperature of delivery inte 
the spectrometer (assumed to be 25° C), an atmos- 
phere saturated with it would contain more than | 
part per 10 million. This statement involves 
number of approximations, including the ease of 
forming ions from the compound, and implies that 
the substance contributes a measurable amount to 
one or more peaks that would serve for the detectior 
and presumably the identification of the compound 
if nothing else were present. It indicat 
any certainty that so small a quantity of the com- 
pound can be recognized in the complicated spectrun 
of the total condensate 


Sensi- 


does not 


In order to make use of this rough generalization 
in anticipating an indication of a compound in th 
spectrum, we need some knowledge of its 
approximate vapor pressure at 25° C. Relatively 
few , Vapor pressure include this 
temperature and those that do are scattered in the 
scientific literature. On the other hand, tables of 
boiling points for thousands of compounds occur in 
several commonly available textbooks The boiling 
point is closely enough related to the vapor pressure 
at another temperature for our purpose provided 
there is not a change of phase between the two 
temperatures. Accordingly, the vapor pressure al 
25° C was plotted with respect to the boiling point 
for enough substances to permit the construction of 
a satisfactory but entirely empirical curve, shown in 
two sections In figure | 


Hass 


observations of 


The compounds included in this curve were chosen 
almost at random, and many of the vapor pressures 
at 25° C were obtained by extrapolating a formula 
riven by the observer. often much beyond the limit 
of applicability specified by him. Little can be said 
for the plot except that the regularity shown is great 
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Picure | Rela nash n the vapor pressure at 15°C’ and the boiling point of com pounds listed in table 1 
enough to serve the purpose intended. Unfor- | because they are rare and because they would have 


tunately, it applies only to compounds that are 
liquid both at the temperatures at which their vapor 
pressures reach | atm and at 25°C. Table Lisa list 
of the included in making the chart. 
Elements, hydrocarbons, and organic halides, fall 
surprisingly close, considering the accuracy claimed 
by the observers of the original data, to a smooth 
curve that a straight line below about 0.05 atm. 
Vapor pressures are notably less than would be read 
from the curve in the ec: of water, the simpler 
aliphatic alcohols, and some other polar compounds. 

Of the 6,175 organic compounds listed in the “C- 
table” of International Critical Tables, not one that 
is liquid at 25° C has a boiling point high enough to 
indicate that failure of its condensate to vaporize in 
the sampling tube would seriously interfere with the 


Both 


substances 


Is 


isc 


determination of its vapor in the atmosphere 


to be widely dispersed to approach saturation in a 
significant portion of the atmosphere, it is improbable 
that high-boiling liquids not on the ICT list would 
occur as vapors in the atmosphere and then be lost 
in sampling. Unless all condensed liquids including 
water are completely vaporized in the sampler, 
solution would complicate the situation somewhat. 
The case is different with solids. No vapor pres- 
sure data were found for any solid that did not appear 
below the line of figure | by a factor in the range from 
1.5 (camphor) to 210° (red phosphorus). All but 
a few chemical substances may be expected to be in 
the solid form when collected at —183° C, and a great 
many of them will remain solid at room temperature. 
Some of these may readily escape detection although 
initially present in the atmosphere in greater quantity 
than more volatile substances that can be determined 
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Whena compound Is identified, the question arises 
the sample that 
the 
original sample total 
pressure ih the sampler and the best avatlable mito 
mation on the vapor compound at 
the working temperature, the sample in the cold trap 
appears to be the amount 
shown by the mass spectrometer may be low by as 
If saturation is 


whether the percentage found in 


same as that i 


enters the spectrometet! is the 
If. from a knowledge of the 


pressure ol thre 
saturated o1 nearly So, 


much as two orders of magnitude 
not indicated in the cold trap, no loss through failure 
to vaporize need be assumed. but there is still to be 
considered the amounts of highly volatile constit 
uents of the atmosphere that were not collected in 
the first place or were lost during the subse puent 
removal of air 
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The problem of loss of substance through failure | &™ ol 
to condense in the cold trap Was attacked in th samp 
same way as that of fathure of the condensate to water 
vaporize, that is, by constructing a plot relating sorbal 
eold U 


vapor IS8° CC to the boiling point — 
Nearly all vapor-pressure data that could be used The 


evel by extrapolation, relate to the liquid phase and 


pressure nt 


} remov 
warm 


(scar 


nearly all substances are solid ut the boiling point of 


OxVeen Liquid phase equations were extrapolated { 
and, as in figure 1, thev gave a very consistent curve } gen. 
which is a straight line within the limits of error of | ™™@° 
the data This line has been drawn in two sections | and tl 
in figure 2 Extrapolated data regarding the vapor | temp 
pressure of solids have been added With the PX. the W 
ception of Tike thane thes lie well below the plotted limed 
line, and solution and adsorption may reduce then (rapp 
furthes They give a rough guide for estimating highe 
the probable loss of material bv failure to condenss nf Ul 
taker 


Thus we should not collect) methan 


unless 10 


expect to 


percent is) present in the atmosphere 
The limit for 
ppm and for ethane 10 ppm 
expect to lose only about 0 y ppm ol acetvlene. 0] 


ppm of hydrogen sulfide, and less than a part pel TI 


4. 


the collection of ethvlene is about 30 


However. we should 


billion of carbon dioxide. The loss of enough of | redu 
anvthing that boils above —20° C to be measurable heig' 
with the mass spectrometer is improbable the 

It is concluded that the delivery to the spectrom- that 
eter of the vapors (1) of all substances that boil abov 3 eX 
about 20° ©, and are liquid at 25° C, and (2) of For 


substances that are solid at 25° C and boil between are 

20° and 300° C, is substantially complete. Out- | 
side this classification there may be significant loss are 
of material at both ends of the scale of volatility. 











{ 
3, Treatment of Sample Before Introduction 
| Into the Mass Spectrometer 


Samples of atmospheric condensates | have been 
ated in varlous Ways before introduction into the 
mass spectrometer. The data and 
norted in this paper will not differentiate between 
‘he yarious treatments because the differences that 
Joguited from them did not produce certainly sig- 
. } nificant changes in the interpretation of the samples 
| to be considered In all cases, the significant portion 
of the condensate must be separated from a much 
eater quantity of air and of condensed carbon 
dioxide and water These major constituents of the 
npolluted atmosphere are to be pri- 
marily as diluents of the sample rather than as 
j ‘oterfering substances. As diluents they carry most 
of the ion current at a few mass peaks and leave the 
remainder of the spectrum too feeble for use. Other- 
wise their interference its limited to a small number 
of peaks shared with pollutants that can be identified 


conclusiens re- 


considered 


most cases, by ions of other masses 
{ir is removed by evacuation while the cold trap 
s immersed in liquid nitrogen is im- 
possible to remove oxygen and nitrogen completely 
without losing low-boiling pollutants, the evacuation 
scarried only far enough to control dilution Car- 
bon dioxide and water have been reduced to workable 
concentrations in several ways. These include the 
ise of absorbing reagents ahead of the cold trap when 
sampling the atmosphere, freezing out a portion of 
water ahead of the cold trap, and the 
sorbants in the path of all of the sample between the 
old trap and the spectrometer. 
The treatment preferred is the following: After 
| removing the air from the cold trap, It Is allowed to 
warm to about 120° C while connected through an 
, Asearite tube to a receiver immersed in liquid nitro- 
ren. When of the dioxide has 
removed in this way, the Ascarite tube is bypassed 
and the transfer of the sample is continued at a higher 
temperature, but enough of 
he water behind as tce. Preferably the ice is sub- 
limed at least once to permit the escape of anything 
trapped in the original precipitate of frost. Later a 
higher-boiling fraction may be taken through a dry- 
ng tube in the same way that the first fraction was 
taken through Ascarite. 
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use of ab- 


most carbon been 


one low Lo leave most 


4. Numerical Arrangement of Analytical 
Data for Study 


The curves plotted by the mass spectrometer are 

| reduced to numerieal form by 
heights” in terms of seale divisions. The heights of 
the peaks considered useful are then ‘“normalized,”’ 
that is, they are added together and each peak height 
$s expressed indi dually as a percentage of the total. 
For this purpose, certain omissions and corrections 
are made. The volumes of air, carbon dioxide, and 
water left in the sample after preliminary treatment 
are each still greatly in excess of the total volume of 
“pollutants,” and their principal peaks are sometimes 


reading the “peak 


too high for the available scale. Even if these peaks 
are readable, they indicate nothing about the sample 
except variations in the success of the pretreatment. 
If they are included in the spectrum when it is 
normalized, two samples that contain identical pol- 
lutants will appear to be substantially different if 
different quantities of air, carbon dioxide, or water 
escaped removal. It is, therefore, convenient to 
include in the total for normalizing only the peaks 
that are subsequently to be used as significant of 
pollutants. For this reason Shepherd and his co- 
workers normalized the California samples with the 
omission of all peaks below (25), those between (27) 
and (37), and (44), (45), (46), (99), (101), and (102). 
To make subsequent observations comparable with- 
out recalculating the voluminous California data, 
the same omissions have been made in the total used 
in normalizing, but several peaks not included in 
that total are sometimes used in the interpretation 
of the results. These include especially (20), (22), 
(31), (35 (45), and (46). Actually, in most 
cases, it makes little difference whether this group 
of peaks is included in the normalizing total or not 
as together they seldom constitute more than a 
fraction of 1 percent. 

The omission of the principal peak of nitrogen and 
carbon monoxide (28), of oxygen (32), and of carbon 
dioxide (44) from the normalizing total and from 
other steps in the interpretation of the results, es- 
pecially the plotting of figures, is unquestionably 
desirable. Little has been made of (29) and 
30), which are large peaks but seem to be nearly 
constant in substantially all the spectra of urban 
atmospheres sO far studied. They are affected by 
isotopes related to (28) and not readily corrected for 
because (28) is an important peak of nitrogen, carbon 
monoxide, and dioxide. Peak (29) repre- 
sents ethyl ions which are possible fragments of a 
great number and variety of heavier compounds and 
have very little diagnostic value. Peak (30) has 
been given much attention, without practical result, 
because it contains the principal identifying peak 
of the nitrogen oxides and acids, and the parent 
peaks of ethane and formaldehyde. Figure 2 shows 
that nitric oxide, NO, like methane, is too volatile 
to be collected in the trap at any probable concen- 
tration but NO, most of which is dissociated to NO* 
in the spectrometer, may be indicated. The peak ts 
usually equal to about 9 percent of the normalizing 
total, and we have been unable to attach any definite 
significance to its variations. Peak (40) Is an im- 
portant one in the spectra of many potential pol- 
lutants and is the parent peak of argon, for which a 
correction must be made if it is to represent pollutant 
and not air. The correction is based on (32). A 
correction of (45) and (46) for isotopes in carbon 
dioxide is based on (22), a doubly charged ion that 
distinguishes carbon dioxide from propane. Usually 
99), (100), and (102) are neglected because of inter- 
ference by the doubly charged isotopes of mereury 
that are alw ays present. Occasionally one or more 
of these peaks is useful. It can usually be ascer- 
tained which peaks are significant in the study of 
pollution because, in the absence of pollution, they 


(36 


use 


carbon 
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should all appear in the constant ratios of the mercury 
isotopes. 

Although several attempts have been made to 
simplify the study of pollutants by fractionating 
the condensate [2], they have not been very useful 
because the separations are not sharp and the spectra 
of the fractions are almost as complex as those of 
the total condensates. Although further study of 
fractionation is planned, all data in this paper refer 
to total condensates only. 


5. Graphic Arrangement of Data for Study 


A great deal of use has been made of the suggestion 
of Sibyl Rock of the Consolidated Electrodynamics 
Corporation that comparisons made between 
spectra by plotting them with the peaks arranged 
in the order of descending height in one of the spectra, 
called the “Reference.” The peak heights are 
plotted on a logarithmic scale 

This method of plotting has several desirable 
features. The reference spectrum approaches a 
fairly smooth curve without maxima or minima. <A 
spectrum compared with it shows differences of com- 
position by deviations from the curve which, if they 
are at all large, take the form of a series of maxima 
and minima. A third spectrum shows similarities 
to the second by the location of the high and low 
points. At any peak the linear distance between 
the two plots represents a ratio between the abun- 
dances of the constituents that produce the ion repre- 
sented. If several peaks, the (117), (119), and (121 
peaks of carbon tetrachloride for example, represent 
only a single compound in both the reference and the 
sample compared with it, they all stand out equally 
from the reference line although they may be of 
greatly different magnitude. The logarithmic scale 
makes relative quantities of rare ions in two spectra 
as prominent as those of more abundant ions. This 
Is especially important because the peaks most useful 
in identification are often small. These features 
help greatly in the identification of compounds in 
what seems otherwise a hopeless tangle of data 
Plots are made on tracing paper so that two or more 
can be read when superposed, and punchings 
facilitate accurate registration. 

In order to make convenient intercomparisons 
between them, all the spectra from the California 
samples were plotted in the order of a single spectrum 
designated “Ref.” Ref was chosen as representa- 
tive of the vers similar spectra of condensates ob- 
tained in Los Angeles and Pasadena during the 
worst smogs that occurred in a period of 3 months 


be 


in the autumn of 1951. The order of peaks in the 
original Ref was as follows: 43, 27. 41. 26. 42. 39. 
57. 56. 55. 70. 25. 58. 69. 71. 38. 91. 40. 7&8. 83. 37. 
S4. 31. 50. 92. Oa. 42. 40, Of, 85. 31, 98, 52, 54. 86. 
67, 68, 77, 111, 82, 65, 112, 49, 63, 95, 46, 106, 79, 
59, 81, 62, 130, 132, 60, 105, 74, 76, 73, 61, 47. 90, 


‘ 
96, 66, 93, 64, 126, 87, 89, 114, 117, 119, 48, 75, 113, 
134, 94, 110, 120 109, 125, 80, 103, 107, 104, 121, 88, 
115, 140. 
Figure 3 represents, by an envelope, the spectra of 
three samples taken in Los Angeles and three in 


\ 


Pasadena. The envelope consists of an Upper | 


joining points representing the highest percenty, 0 
. . ay 


of each ion in any of the six condensates and a Joy 
line representing the lowest percentage. The brok 
line between the two is Ref The croup of hie 
peaks (117), (119), and (121), near the lower o, 
of the plot represent carbon tetrachloride jn o, 
sample, and the peaks (64) and 48) represent sulf 
dioxide in another. Both may be the result qi 
contamination of the samples in the laboratory. 
Samples taken later in various places produg 
other peaks below the level of in t 
samples used in the original selection of Ref. — It wa 
convenient to adopt a standard order for plotti; 
these and other peaks, such as (29) and (30), thar 
had been disregarded initially. This order js gp. 
proximately that in which the peaks happened {| 


detection 


become of interest, and is as follows: 128, 194 
123, 127, 138, 139, 131, 129, 142, 136, 137, 118. 12 
141, 108, 116, 147, 154, 151, 152, 156, 148, 149 
143. 158, 146, 155, 145, 29, 30, 31, 35, 36, 99. 109 
164, 165, 166, 168, 170. In the figures that follow 


PERCENT 


this order w ill be considered a part of the “order of 
Ref,” not all of which is alwavs used. 

The plotting of peaks in other orders than that of 
Ref is frequently desirable. It is now the usual 
practice to plot a spectrum first in the order of Ref 
for comparison with plots of other mixtures in order 
to make as good a preliminary estimate as possible of 
what is present. The sample under study is then 
plotted in the order of magnitude of its own peaks 
and the spectrum of the components thought to ly 
present is built up as they are tentatively identified 
This procedure gives most of the advantages of 
comparision with Ref, particularly the recognitior 
of groups of peaks that indicate a single compound 
and permits the modification of the tentative analysis 
as often as necessary with a minimum of replotting 


ION CURRENT, 


TOTAL 


6. Regularities Found in Samples From 
Southern California 


The most remarkable thing about the condensates 
obtained from the atmosphere of Southern Cali- 
fornia during smogs, was the close resemblance be- 
tween their mass spectra. This is strikingly shown 
by the narrowness of the envelope of the stx samples 
represented in figure 3 taken in different parts of the 
area during a period of 3 months. Although these 
samples were admittedly selected because of theu 
close agreement, many other samples showed strong 
resemblances to the group and tended to show a 
progressive series of differences. By this is meant 
that the deviations from Ref were random 
Instead, differences usually appeared 1 much th 
same series of peaks, and if a certam peak height 
differed from that of Ref more in one sample than in 
another, several other peak heights showed  cor- 
respondingly greater deviations in the first sample 
than in This would. of course, neces- 
sarily be the case if the peaks considered represented 
a single compound; but the relationship appeared to 
be too general to be explained by different  con- 

2 Two explanations 


not 


the second. 


centrations of 1 or 2 substances 
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eement between spectra of several condensates 


appear probable, 1) that the atmosphere contained 
varying proportions of two principal groups of pol- 
lutants from different 2) that a pre- 
dominant pollutant from a single source was under- 
going change by reactions in the atmosphere, the 
progress of which was indicated by the changing 
spectrum. It is not surprising that local or tem- 
porary variations confused the picture considerably. 

Accompanying the samples were some rather 
vague statements about weather conditions, v isibility, 
irritation, ete., from which it was possible to 
arrange the numerous samples roughly in the order 
Hach of numerous 


sources, or 


| 
eve 
of apparent severity of smog 


£4157 ‘7 


take n d iring pe riods of he avy smog in Los Ange le 8 and Pasade na, 


, Ref 


peaks and functions of peak heights, especially ratios 
between selected peaks, were plotted with respect to 
the numbers of the sample arranged in this order of 
apparent severity. Several of these plots indicated 
fairly definite trends; for example, the most abundant 
peak, (43), appears to increase during smog. How- 
ever, the most positive trend appeared among the 
three ratios, (58)/(53), (72)/(67), and (86)/(81), all 
of which increased strongly and fairly consistently 
with increasing indications of smog. A _ further 
discussion of regularities within groups of spectra 
is given after the next section on “‘normal traffic gas.” 
The trends mentioned in this paragraph all relate 
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to the composition of the mixture of pollutants in the 
air and not to its over-all concentration, with which 
no certain correlations are apparent. 

When studies of pollution in other cities were begun 
it was immediately evident that the spectra of their 
atmospheric condensates usually resembled strongly 
those from Los Angeles. The obvious explanation is 
that a large part of the pollution came from a similar 
and when weather conditions and 
places of sampling are also considered, there appears 
to be nothing common to the probable origin of the 
pollutants except automobile traffic 

Both exhaust gas and gasoline vapor are chemically 
so complex that their determination as individual 
compounds is impossible by any means now available 
But because they are released from so many vehicles, 
they are probably much the same everywhere; and it 
should be possible to treat “traffic gases’? (a term 
intended to include an average amount of gasoline 
vapor from filling stations, leaks, etc.) as a unit 
In order to recognize pollutants of a different 
character and to locate their sources, it is important 
to differentiate them from traffic gases as a whole 
Unless we can do this, by mass spectra or otherwise, 
we lack the information needed to deal with the 
control of pollution at its various sources 


source or sources; 


7. Mass-Spectra Pattern of Normal Traffic 
Gas 


The first step was to determine with a fair degree 
of certainty, though not necessarily with great 
accuracy, the mass spectrum we may expect from 
traffic gas. For this purpose we had available the 
spectrum (A) of a sample taken from a tunnel in 
Pasadena in 1951 during heavy traffic, and another 
(B) taken in a tunnel in Washington, D. C. in 1954 
under similar conditions. The Pasadena sample 
represented both exhaust gases and vapors from the 
engines, and the generally polluted atmosphere of 
the city drawn into the tunnel as the cars passed 
through, but it was reasonable to suppose that the 
high concentration of traffic would result in 
the outside pollution becoming negligible or nearly so 
Actually, the one or two recognizable constituents 
of outside smog not attributable to motors indicated 
that about 20 percent of the total condensate 
collected came from outside The Washington 
sample (B) was taken in a region free of industrial 
contamination and was believed to be almost free 
of pollution from any source except automobiles 
Both these samples represented continuously moving 
traffic. To supplement them three samples were 
taken; (C) from a busy intersection where traffic is 
interrupted by lights; and (D) and (E) from points 
of heavy but uninterrupted traffic in residential 
districts. <A link of a sort was established with the 
Southern California samples by plotting the ratios 
between certain peaks previously mentioned as best 
indicating the appearance of smog The several 
samples, including two from Detroit and four from 
Washington. were arranged chronologically the 
order of sampling. The result figure 4. The 
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FIGURE 4. Variation in the ratio of the peak height of thre 
pairs of mass numbers as affected by the nature of the sample 
S, Southern California; D, Detroit; W, Washington 
Samples used in deriving patterns: R, Ref; T, NTG 


ordinates are logarithmic, and the curves are sepa- 
rated by shifting the seale 1 decade for each. Hori- 
zontal lines represent the ordinate 1.0 for the cury 
that intersects it The fact that the three ratios ar 
closely related without regard to locality is up- 
mistakable. The six samples used to establish Ref 
fig. 3) are indicated in figure 4 by the letter R, and 
the five sanrples chosen to represent normal traffi 
gas, by T 

The five spectra of traffic gas were plotted and 
critically examined, point by point. The agree- 
ment was not as good as had been expected, and the 
reasons for disagreement are obscure. Hlowever, it 
could be assumed that if such a thing as a normal 
traffic gas exists, there is nothing in any of the sam- 
pling locations to remove anything from the mixture 
but it is possible, even probable, that significant addi- 
tions may On the other 
hand, there is the possibility that something was lost 
from any of the samples during handling, by re 
action, solution in stopcock lubricants, or the like 
In view of these possibilities it was not believed that 
statistical rules of probable error would have much 
significance, and it decided consider each 
mass number individually and to select what seemed 
to be a “best value” rather than to take any kind 
of mathematical mean The result ts plotted with 
peaks in the order of the customary Ref, in figure 5 
The solid line is the selected spectrum, which will be 


come from other sources. 


wis to 


referred to as normal traffic gas or NTG. Ref 3 
shown as a dotted line. It must be admitted that 
initially no great confidence was placed in the 


accuracy of this picture of normal traffic gas, but 
subsequent use has caused greater significance to be 
attached to it 
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Figure 5 Spectrum of normal traffic gas plotted in the same order as Ref 


NTG 


‘8. Comparison of Normal Traffic Gas with 


Ref 


A comparison of the two curves of figure 5 at once 
shows that compounds of certain chemical series are 
more abundant in one than in the other 

The members of a homologous series of compounds, 
for example, the paraffin hydrocarbons, differ in 
molecular weight by multiples of 14. Thus methane 
has a molecular weight of 16, ethane 30, propane 44, 
butane 58, ete A positive ion formed by the loss 
of an electron from the molecule, called the parent 
ion, has the mass number of the original molecule 
However, many more 1ons formed from most 
compounds by the breaking up of molecules than by 
the loss of single electrons. An examination of the 
structural formula of any paraffin hydrocarbon will 
show that if it is broken in any manner into two 
fragments and no more, each fragment will have a 
mass of | or of | plus a multiple of 14; and that mass 
will be 1 less than one of the parent masses of the 
Thus. when butane loses 
is left, if 
methyl group) ts 


same series of compounds 

hydrogen atom, mass (57 1 carbon atom 
with its attached hvdrogen a 
detached, the residue has mass (43), 1 unit less than 
the parent mass of propane detaching an ethyl 
group leaves 29, 1 unit less than the parent mass of 


ethane. ete Although ions are also formed with the 





. Ref 


breaking of more than 1 bond of the molecule, the 
ions of the parent mass and of 1 less than the parent 
mass of some member of the homologous series, are 
predominant enough to be significant of the chemical 
nature of the parent compounds. 

With this in mind, table 2 was prepared to afford 
a comparison of the chemical nature of the com- 
pounds represented by NTG (normal traffic gas) 
and Ref (the condensate from smog, assumed to be 
formed in considerable part, at least, by chemical 
reactions of traffic in the air). It should be 
noted that the same series of numbers may include 
the parent peaks of homologous series of compounds 


gas 


of quite different character. 


In the following discussion it will be assumed that 
Ref is the product of reactions of NTG in air rather 
than a mixture with things from other sources, 
although this can be only roughly true, and the 
mixture represented by Ref will for brevity be refer- 
red to as “oxidized”’ as distinguished from ‘‘normal’’ 
traffic gas. Table 2a indicates that compounds of 
the same parent masses as the paraffins increase in 
relative amount as the atmospheric reactions pro- 
ceed. Because all numbers in the table are percent- 
ages of a total sample, the simplest explanation for 
the over-all effect is that compounds other than the 
paraffins are disappearing, probably with the for- 
' nonvolatile particles or water, carbon 


mation of 
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dioxide, carbon monoxide, etc., not included in the 
part of the spectrum considered, leaving the inert 
paraffins as a greater percentage of the total. How- 
ever, we find that in the aggregate it is only the 
parent peaks that have increased. If the first ex- 
planation were the only one, the odd- and even- 
numbered peaks should have changed in proportion. 
A more interesting explanation than the first is that 


aldehydes and ketones are among the products of 


the smog-forming reactions. Aldehydes and ketones 
produce a much larger proportion of parent ions 
than do the easily broken paraffins. 

Table 2b shows the surprising fact that in general 
the peaks thought to be most characteristic of the 
olefins remain practically constant in relattve abun- 
dance. In the literature bearing on the Los Angeles 
smogs the unsaturated hydrocarbons, and the olefins 
in particular, are given a major part of the blame 


for smog formation, and irradiation experiments of 


olefins in air have shown that the blame is not mis- 
placed. At present no adequate explanation of the 
relative constancy of the “olefin peaks’ has been 
thought of by the authors. One possibility is that 
most of them are derived from cycloparaffin rather 
than olefin sources and that the cycloparaffins remain 
substantially inert in the atmosphere. 

Table 2c shows a striking decrease during smog 
formation of nearly all the parent and parent-minus- 
one peaks of the group that includes acetylene and 
the dienes. The only exception is acetylene itself, 
which remains relatively constant as indicated by 
its parent peak, (26). The relatively smaller num- 
ber of (25) ions in traffic gas may well be the effect 
of a reduced number of fragments of heavier mole- 
cules rather than a decrease in acetylene itself. It 
seems reasonably certain that the decided decrease 
in the other peaks in this group during conversion 
to smog is the result of reactions in which their 
parent molecules are used up. 

Table 2d shows a decrease in the 
hydrocarbons only less striking than in the preceding 
group, especially in the higher members of the 


aromatic 


series. 

Table 2e represents the 
of a rather wide variety of oxygen-containing com- 
pounds, of which the alcohols, and particularly 
ethyl alcohol are likely to be the most abundant. 
Generally the peaks are not in the ratio that would 
indicate ethyl aleohol alone, however. This series 
shows such striking increases that it most 
probable that it is derived in large part from the 
products of traffic gases reacting in the air. 

Many of the peaks in figure 5 have been labeled 
Those of the series beginning with (25) and (26 
are marked with a line above the number, those in 
the series beginning with (32) and (31) are under- 
lined. With the exception of (25), all the 16 peaks 
of the first series are above Ref by surprisingly 
consistent amounts. A secondary regularity is to 
be seen in the fact that, on one side of peaks (53) 
and (54), which are about equally above Ref, the 
even-numbered peak is larger than its adjacent 
odd-numbered peak, whereas above (54) the re- 
is true. Peaks (95) and (96) seemed out of 


peaks characteristic 


seems 


verse 18 
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place until it was discovered that those numbers 
in Ref were strongly affected by trichloroethylene. 
When they were corrected for the chlorine compound, 
they fell in the positions on the descending Ref 
yrve marked by (), which makes them consistent 
ith other members of the group 

Again with one exception, (74), the underlined 
numbers are all below the line. Peak (74) is prob- 
ably affected by an interference, as (95) and (96) 
were, though the source of the interference has not 


heen determined. 


9. Estimation of the Kind and Amount of 
Traffic Gas in a Local Situation 


It would be unreasonable to suppose that all the 
regularities occurring in the relation between NTG 
and Ref could have appeared had not the selection 
of the two patterns been substantially representa- 
tive of two general phases of air pollution. Con- 
tinuing the assumption that smog is derived prim- 
arily from traffic gas, a pollutant having the spec- 
trum of NTG should be converted, in time, to one 
with the spectrum of Ref. These two spectra are 
so different that it is at least desirable to distinguish 
between them whenever a pattern for a traffic gas 
is used in the study of an unknown sample. 

It has become a fairly well established practice 
in the petroleum industry to determine the various 
types of compounds in petroleum products, paraffins, 
aromatics, cycloparaffins, olefins, ete., by consider- 
ing ratios between selected groups of peaks; and 
reported comparisons with synthetic mixtures show 
that high accuracy is frequently obtainable. An 
attempt was made to apply a similar method, 
using the ratios represented in figure 4, to determine 
the state of reaction of the pollutants from motor 
traffic. 

The recorded values of the several peaks (with 
one addition) and the computed ratios are shown 


in table 3. 


TABLE 3 Peal princi pa ly used in sdging the extent of 
oxidation of traffic ga 
Peak NT¢ Re R NTG 
NTG Re 
1 0. 66 2 4! 2.4 
s) ‘4 1. 40 2. ¢ 0. 38 
1. 20 { 0 2.2 
2 . 33 ‘0 , 0. 66 
67 AN 2 0. 28 3 
st Is 1.4 0.71 
81 10 121 0. 30 3.3 
_~ : 
72 4 ts 9 
~“ s! 15 
Aver { 2 


Two approximate formulas were developed for 
estimating the extent of oxidation, but the rigid ap- 
plication of neither of them was very successful for 
two apparent reasons, (1) pollutants unconnected 
with traffic gas interfere with some of the peaks used, 
and (2) the reactions of oxidation do not all oecur 
at the same rate. A rough estimate of the extent of 
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oxidation can usually be made by plotting the spec- 
trum of the sample in the order of Ref and comparing, 
preferably by superposing, with the spectra of NTG 
and Ref as shown in figure 5. 

The total percentage of traffic gas is estimated at 
the same time. It will be noted that the portions of 
Ref and NTG between ordinates 10 and 1.5 (peaks 
27 and 25 in the usual order) coincide very closely. 
The plot of the sample is adjusted vertically to obtain 
as good a fit of this portion of its curve to those of 
figure 5 as is practicable but with one limitation, 
sample peaks in other parts of the spectrum must 
not fall so low that their explanation would require 
the elimination of an unreasonable fraction of the 
traffic gas. Peaks much higher than either Ref or 
NTG should not be disturbing, in fact they are what 
we are looking for. When the two sets of curves have 
been fitted together as well as practicable, the vertical 
distance between the indices (ordinate 1) of the two 
scales shows the percentage of the total sample to be 
attributed to traffic gas, normal or oxidized. The 
fraction ““X,”’ the percentage of conversion of normal 
to oxidized traffic gas, is estimated by noting the re- 
lation of peak heights of the sample to those of Ref 
and NTG that are widely separated from one another. 
Peaks (53), (54), (58), (67), (72), (81), and (86) ap- 
pear to be the most useful for this purpose. If 
maxima are not noticeable at (53), (54), (67), and 
(81), oxidation of the sample has proceeded far. 
Minima at (57), (72), and (86) are less significant 
because the presence of interfering substances accen- 
tuates common peaks but raises low points of the 
spectral curve. 

Some of the peaks found in the spectra of the traffic 
gases are indicative of known compounds with very 
little probable interference. This is especially true 
of the parent peaks of the aromatic hydrocarbons. 
The aromatics are important constituents of the 
traffic gases and should remain a part of the total 
spectrum of a sample until a decision has been made 
as to its probable content of traffic gas. However, 
the aromatics can be independently determined and 
it would be unfortunate to have the analysis affected 
by their amounts in Ref or NTG. Hence the pat- 
terns of the aromatics are “‘taken out” of the initial 
patterns of Ref and NTG, and are applied separately 
in the interpretation of the spectrum of the sample. 

The chlorine compounds trichloroethylene and 
carbon tetrachloride were present in all samples used 
in the selection of Ref and in 2 or 3 of those used for 
NTG; but they are not considered normal constitu- 
ents of traffic gas and they have also been taken out 
of the patterns of Ref and NTG. With somewhat 
less assurance, the patterns of ethyl alcohol, acetylene 
and sulfur dioxide were also taken out of the traffic 
gas patterns, but the spectra of these compounds are 
subject to so many interferences that they are deter- 
mined only after most other identifiable constituents 
have been taken out. 

While the spectrum of the sample is superposed on 
figure 5, a notation is made of any peaks that seem 
abnormally high, and an effort is made to identify 
these peaks, as far as possible, with the chemical 

they most likely to represent. 


compounds seem 



























Sometimes they indicate a whole group of comp rubnidds 
such as gasoline vapor, in excess of that 
present in traffic gas. 

The procedure is now to attempt to interpret the 
spectrum of the sample in the usual way except that 
the expurgated patterns of Ref, NTG, and other 
approximately known mixtures are 
though they represented single substances. It is 
helpful to plot the spectrum of the sample in the 
order of abundance of its own peaks, and to compare 
with it from time to time the composite spectrum 
of whatever seems to have been identified. Rela- 
tions among the residual peaks are more r “aclily seen 
in graphic than in numerical form, and when the 
analyst has done the best he can, the remaining plot 
indicates at a glance how his effort 
been 


10. Application of the Method to a Specific 
Sample, Louisville D 6 


1 
normally 


treated us 


has 


successful 


Even with the general method of approach that 
has been outlined, each sample presents a unique 
and complicated problem that is solved largely by 
trial and error. The interpretation of the spectrum 
of one sample, designated D-6, from Louisville will 
be followed in considerable detail to illustrate the 
method employed 

The spectrum of the sample was first plotted in 
the order Ref and the graph superposed on the 
patterns of both Ref and NTG, shown in figure 5 
The result is figure 6, in which the spectrum of the 
sample is represented by a solid line and the two 
types of traffic g broken ones From the 
agreement of the half dozen larger peaks, it is evident 
that the sample contains about the percentage of the 
more abundant hydrocarbons that is usually found 
in traffic fas However, the plot also indicates that 
trafic gas contains much greater quantities of a large 
number of other things than are found in the sample 
and these cannot be accounted for except by assum- 
ing that (1) there are great differences in the com- 
position of the traffic-gas in Louisville and elsewhere 
or (2) that much of the major hydrocarbon fraction 
of the sample was derived from something other than 
traffic gas. 

Rejecting the first assumption, the relative posi- 
tions of the two plots were shifted to those of figure 7. 
in which there is as good agreement among minor 
constituents represented by the numerous smallet 
peaks as seemed readily attainable. The relation 
between the indexes of the two plots (ordinate 1 of 
NTG is in parentheses) shows that the dotted curve 
now represents the contribution of 30 percent of 
NTG to the total spectrum. Ref has been omitted 
from the figure because the sample was taken early 
in the morning when much conversion to the condi- 
tion of Ref seems unlikely and because the oceur- 
rence of high and low points in the two spectra at 
the same mass numbers is frequent enough to indicate 
normal rather than oxidized traffic gas. 

If only a third of the principal group of hydro- 
carbons is from traffic gas, it is necessary to account 
otherwise for the other two-thirds 
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probable that it was raw gasoline vapor from a source 
not at first identified. There was available for com- 
parison the spectrum of a composite sample of 
“average’’ gasoline from Los Angeles, made in 1951 
This has a strong resemblance in the proportions 
among its higher peaks to the spectrum of both the 
sample and NTG, but the percentage of these peaks 
is somewhat higher than in the sample. A matching 
of the curves as well as practicable by superposing 
them in the same way that the sample and Ref were 
previously matched, indicated an excess in the gaso 
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line of about 25 percent. From these data it 


appeared that in addition to 30 percent of NTG the 
sample should contain about 56 percent of gasoline 
vapor. A first calculation gave generally high values 
for the principal peaks that were high in the gasoline 
pattern, and the fraction assumed to be present was 
then rounded off to 50 percent 

In the remainder of this discussion the interpreta- 
tion of the spectrum is represented by figure 8. The 
sample spectrum was replotted in the order of magni- 
tude of its own peaks, The open circles indicate the 
contribution to the spectrum of 30 percent of NTG, 
the crosses represent the sum of 30 percent NTG 
and 50 percent of Los Angeles gasoline vapor, and 
the solid dots represent the sum of everything even- 
tually added together in the effort to account for the 
entire sample. 

‘the figure a glance. The 
distance between a cross and the solid line indicates 
material to be accounted for from sources other than 
traffic gas and gasoline. Most of the masses at 
which there were large gaps to be accounted for were 
found to be important peaks of some of the com- 
pounds listed in advance as probable emissions from 
the industrial district called Rubbertown. Thus 
peaks (26) and (25) surely represent acetylene, (88) 
represents chloroprene, ete The | 


shows a great deal 


substances 


show that the 


gasoline only 


of those 


identification of 
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that cannot be attributed to common 
traffic gas and gasoline, is the most important reason 
for assuming and applying a pattern for traffic gas. 


The substances identified are summarized in table 4 


It will be observed that most of the totals that are 
much above the line coincide with the crosses and 
sample did not produce particular 
peaks as large as was expected from traffic gas and 
They are to be interpreted as mean- 
ing that the traffie gases and motor fuels whose pat- 
terns were used did not have exactly the composition 
1 Louisville. This is not surprising. It is 
gratifving that such local differences affected as little 
as they did the identification and approximate deter- 
mination of pollutants from other sources. 

The peak heights of (78), (91), and (92), obtained 
by adding the patterns of NTG and gasoline and 
indicated by crosses, are above the dots representing 
“secounted for.’ The reason for this is that (78), 
the parent peak of benzene, and (92) the parent peak 
of toluene have no probable significant interferences, 
amount recorded by the spectrometer must, 

be considered to be the amount in the 
The dot was therefore placed on the line 
representing the sample. None of the aromatic 
peaks higher than (92) was reported; accordingly 
(91) accounted for was calculated from (92) and the 
pattern of toluene. It falls somewhat below the line, 
indicating that there was a certain amount of aro- 
matic gage of higher molecular weight, but not 
enough to be reported. It is much below the (91) 
to be die ‘ted from NTG alone. Peaks (67) and 
(81) are high, possibly because the butadiene which 
was one of the gases from Rubbertown may have 
contained some higher homologs. With two excep- 
tions, (71) and (85), the other peaks notably above 
the line are, in the order in which they occur in the 
figure, (55). (56). (70). (69), (83). (68), (98), (82), 
(97). (79). (112). (96), and (111). All of these could 
be olefin peaks, but are more likely to be peaks of 
the cyeloparaffins. Their scarcity in Louisville prob- 
ably results from the fact that the California petro- 
leums and those of Venezuela from which Washington 
is probably largely supplied, have high cycloparaffin 
contents that are lacking in the midcontinent petro- 
leums. It was found that samples from Detroit and 
Cincinnati showed a similar dearth of cycloparaffin 
peaks, though not . the same extent as this sample 
which was so largely raw gasoline. It was noted 
that the ale of accounted for peaks to those of the 
sample tends to increase with increasing mass. A 
possible explanation for this is the fact that most of 
the gasoline in Louisville probably came from a tank 
farm in which the vapor escaped above the liquid 
and was generally lighter than the liquid left behind, 
whereas the gasoline used as a pattern in the ac- 
counted for spectrum was completely vaporized. 

Three chlorine compounds, chloroprene, vinyl chlo- 
ride, and dichloroethylene, were easily identified. 
The first two were on the list, with hydrogen chloride 
and chlorine itself, of compounds to be expected 
from Rubbertown. Because acetylene and chlorine 
readily form dichloroethylene, this compound may 
be the product of reaction in the atmosphere. 


and the 
therefore, 
sample. 



























The dots significantly below the line represent 
portions of the sample spectrum, (37), (49), (64 
(35), (36), (47), (76), (87), and (89) still not com- 
pletely accounted for, but all nearly all of them 
are peaks that occur more abundantly in the spectra 


of chlorine compounds than elsewhere, and they seem 
to represent a background of numerous chlorine de- 
rivatives present only in very small traces. 

Conclusions regarding sample D—6 are summarized 
in tables 4 and 5 

The relative sensitivities of compounds are usually 
given in the tables of their patterns, hence relative 
concentrations of the identified chemical compounds 
need not be much in doubt. They are shown for 
sample D-6 in table 5. It is a majer weakness of 
the method of analysis that no definite sensitivity 
can yet be assigned to the patterns of traffic gas and 


TABLE 4 Inte rpretation of Louisville sam ple D-6 


General fractions \fole percent 


Traffic gas Uw 


Gasoline vapor LAL 
Principal a. - 
Identified compounds identifying neignt, 
of total 
peak 
spectrum 
One-ring aromatic hydrocarbons 
Above C3H 10 Nome None 
X viene, CsH 106) race race 
Toluene, C;He v2 0. 109 0.17 
Benzene, CyH 78 223 2 
Terpenes us None No 
Styrene, CaHs 14 Nome N one 
Carbon tetrachloride, CCl 117 None Nome 
Trichloroethylene, C,HC] 130 Nome Nome 
Chloroprene, CyHsC1 al 0 665 3.71 
1,3-Butadiene, C,H, 4 208 0.74 
Acrylonitrile, CyHyN i O51 lf 
Vinyl acetylene, C,H 52 123 tt 
Vinylebloride, CyH yC! 62 Os 9 
Dichloroethylene, CoH )C! il 151 24 
Acetylene, C,H 2h 6.79 10. 31 
Ethanol, C,H,O 31 1. 32 2.4 
Sulfur dioxide, 8O 1s 0. 106 { 
Hydrogen chloride, HC! a 141 ‘ 
Chiorine, Cl 4 1h 
* Includes ethyl! benzene 
Not quantitatively significant of concentration in the original atmosphe 
because a large fraction is lost with the carbon diovxids Chiorine and ilft 
dioxide are entirely incompatible in aqueous solution and their coexistance 
the sample to the point of introduction into the spectrometer n some doubt 
Although chlorine and hydrogen chloride were on the list of pollutants discharged 
into the air by industrial plants, it ilso uncertain that (35) and (36) represent 
these two gases rather than unidentified organic compounds contair 
TABI } 5 Appro rmate concentrations of ide nlitied compound 
in Louisville sam ple D ‘) 
Column A— Height of identifyin Ji ik of total ion current 
Column B—Ratio of sensitivity of identifying peak to base peak 
Column C Ratio of sensitivity or ise peak to base peak of butane 
Column D—A number proportional to mole fraction of identified com, 
A/BC 
Column E Mole fraction of pollutant if identified compound are 
ggregate 2 
Compound 4 B ( I) } 
Toluens 0. 109 0.77 1.17 0.12 0.17 
Benzene 223 1.00 1. St lf 
Chloroprene Phy 1.00 (0). 2 2 55 1 
1,3- Butadiene 208 0. 88 67 0.51 i 
Acrylonitrile O51 OO ) 11 ! 
Vinyl acetylene $23 1.00 ) 15 ti" 
Vinyl chloride O95 0. 76 74 17 
Dichloroethylene 151 1.00 74 " Pal) 
Acetylene, total 6.79 1.00 " 7.08 
(from industrial sources 
Ethanol, total 1. 32 1.00 79 1. 67 2.43 
from industrial sources j 
Sulfur dioxide 0. 106 0.49 { tt 7 
Hydrogen chloride 141 1. OM ‘ a) 8 
Chlorine 117 
rotal 13.74 20). OM 
* Not quantitative gnificant for the reason stated in footnote (b) of table 4 


gasoline, and this makes it impossible to determi, 
with accuracy even the relative mole fractions 
these mixtures and of the individually identifie 
compounds. In pre paring tab le 5, it assumed 


was 


that the 13.74 percent of the ion current carried by 
‘ide ntified” compounds represented a 20 mole per. 
cent of the s ample, and that the remaining 869 


percent 
percent of traffic gas and gasoline. It follows tha 
for the traffic gas and rasoline BC= 1.08; that js 
unit volumes of their mixture contribute as much to 
the ion current in the spectrometer as 1.08 volumes 
of butane. 


11. Louisville Sample ST-1 


Other samples from Louisville can be reported jp 
much less detail. Several of them present points of 
special interest, however. 

Sample ST-1 was taken from 9 to 10 p. m. in the 
West End of Louisville at a time of ve ry little alr 
movement, the direction of which was not de ‘finitely 
determined. The ratios plotted for other samples ip 
figure 4 and often used in judging the percentage of 
change from NTG to Ref, were all lower than for 
NTG itself. This, and the fact that the sample was 
taken after sunset combine to make any previous 
photochemical reaction in the atmosphere improb- 
able. From a comparison of a plot of the sample 
with one of NTG, it was estimated that 85 percent 
of the sample was traffic gas, but as the interpreta- 
tion of the spectrum progressed this estimate was 
changed to 75 percent. 

The procedure illustrated in figure 4 was followed 
The spectrum of the sample was plotted in the order 
of its own peak heights and initially compared with 
the sum of NTG and the aromatic 
indicated by ‘their parent peaks. 


hy crocarbons 


Among the peaks 


prominent as “unaccounted for’ at this stage were 
(93), (94), (107), and (108) which were thought to 
represent phenolic compounds; (47), (49), and several 
other pairs of peaks two mass units apart, inter- 


preted to represent chlorine compounds; and (55), 
(69), (70), and others that seemed to represent 
a concentration of olefin or naphthene compounds 
much too great to be explained by traffic gas. 

It appeared probable that the prominent 
and (93) peaks should be attributed 
CH,C,H,OH. No pattern was available 
sample of commercial ‘‘Cresol, U.S. P.’’ was run to 
provide one. This gave a spectrum containing 
some ions too heavy for pure cresol and much more 
(94) than seemed likely to be formed in quantity 
from a compound with a parent mass of (108). It 
was, therefore, decided to take out the pattern of 
phenol from that of the Cresol, U P. and to 
consider the remainder be the pattern of every- 
thing in cresol heavier than phenol. Phenol and 
based on (94) and (108), were then taken out 
of the spectrum of the sample. There still remained 
a large (93) be accounted for. Terpenes, in 
particular the pinenes, seemed to offer the only 
patterns that would account for (93) without large 
deficits in other peaks, and the (93) was accordingly 
attributed to pinene 
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Two other easily recognizable but complex groups 
of compounds were obviously present: (1) Chlorine 
derivatives and (2) unsaturated hydrocarbons, appar- 
ently mostly of the olefin ’ but with some | 
acetylene or diene members. — These were clearly 
not from the same sources in Rubbertown that 
supplied the chief pollutants of sample D-—6; for | 
example, chloroprene was present only as a trace, if 
at all, and there was no evidence of acrylonitrile, | 
vinv! compounds, or of gasoline in excess of that to 
be attributed to traffic gas. 

The possible number of compounds that might 
have been constituents of either the chlorinated or 
the unsaturated groups of pollutants is many times 
greater than the number of peaks available for their 
identification; a complete analysis in terms of indi- 
vidual chemical compounds is accordingly out of the 
question, just as is that of the hydrocarbons of 
However, some of the constituents are 


series 


traffic gas 
almost certainly identifiable, and the general nature 
of the two groups seems fairly certain to be that 
indicated in table 6. both in chemical character, and 
approximate molecular weight 

The actual order of taking out 
stances was as follows: Cresol, phenol, CoHy, 


the various sub- 
and 
other aromatic hydrocarbons heavier than xylene, 
C.Hy» representing the xylenes and ethylbenzene, 
styrene, toluene, benzene, N TG, a mixture of alpha- 
and beta-pinene in equal proportions, ethanol, 
dichloroethvlene, dichloromethane, methyl chloride, 
trichloromethane, dichloroethane, dichloropropane, 
3-methyl-l-pentene, 3-hexyne, cis-2-pentene, and 
finally isohexane (2-methyl pentane After the last 
of the unsaturated hydrocarbons had been taken 
out, there remained prominent peaks at (86) and 
71) in a ratio that corresponds to isohexane, which 
was then taken out. To have attributed (86) and 
71) to gasoline vapor of the composition assumed 
in the case of sample D-6 would have produced 
serious deficits in other peaks, particularly (72) 
Several other possibilities were considered but led to 


lle sample ST-1 


PaBLe 6 Interpretation of spectrum of Louis 
P ’ Peak Mole 
S ‘ he ht, ry 
total 
peak pectrum 
Normal traffic ga 
Phenolic compound 
Cresol, C;HsO 108 0. 552 1 84 
Phenol, CsH»O r 151 0. 66 
One-ring aromatic hydrocarbon 
Above, CcH Ww O73 97 
Xylenes, CsHy» and ethyl benzen 10 360 67 
Toluene, C;H« 2 1. 6s 2 OF 
Benzene, CoH 7s 1.79 1. 32 
Styrene, Css 104 0. 030 0. 02 
rerpenes 
Pinenes, CyoH ; Ty 78 
Unsaturated hydrocarbons not a 
counted for by NTG 
3-Methyl-l-pentene, C,H “4 i) 1. 80 
}-Hexyne, CyH 2 336 0. 61 
cis-1-2-Pentene, C;H 1) FM) 167 
Isohexane, CsH 71 "4M) 1. 20 
Ethanol, C,_H;OH 7 1.92 2 O4 
Aliphatic chlorine compounds 
Methyl! chloride, CH.C1 0 0. 40 0. 7 
Dichloromethane, C H.C} 410 oA) 42 
Chioroform, CHC] a3 p24 | 1. 23 
Dichloroethane, C.H«C} 62 303 0.51 
Dichloropropane, ¢ HeC] ill 126 13 
Dichloroethylene, C H.C} } Tt 14 
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similar difficulties. It seemed most likely, therefore, 
that there existed some source of fairly isolated 
isohexane, perhaps used as an intermediate in chemi- 
cal manufacture or as a solvent. 

Other residuals are either too small to be of probable 
significance, are of uncertain origin, or are probable 
fragments of parent molecules too heavy to have 
been included in the recorded spectrum. Thus (58) 
may represent propionaldehyde, acetone, or un- 
accounted for butane; (81) is probably a fragment of 
an unidentified diene or acetylene derivative with a 
molecular weight of (96) or more; and (95) and (97) 
are probably derived largely from unidentified 
chlorine compounds with at least 2 carbon and 3 or 
more chlorine atoms. 

Figure 9 shows the degree of agreement reached in 
accounting for the sample, in this case by the relation 
of solid dots to the unbroken line. As a matter of 
interest, a broken line was added to show the spec- 
trum of D-6 plotted in the same order. The large 
differences between these two samples is in striking 
contrast to the generally close agreement among 
southern California samples, illustrated in figure 3. 
The numbers above the solid line identify peaks 
much higher in D-—6 than in ST-1, those below the 
line those that were much lower. Numbers below 
breaks in the dotted curve represent mass numbers 
not reported in sample D-6. It is probable that the 
small peaks attributed to fragments of unidentified 
chlorine compounds in D-6 came from the same 
general source as the aliphatic halogen derivatives 
in ST—1 during a brief shift of wind direction. 


12. Analysis of Louisville Sample X~—1 


Sample X-1 (fig. 10) shows several points of 
interest. Whereas the major hydrocarbons’ pattern 
of ST-1 is well accounted for by normal traffic gas, 








or 
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Fiaure 10 Total contribution of identified substances to the 
spectrum of Louisville sam pe Y—/ 

, Sample swccounted for 
and that of D-6 indicates the presence of a major 
amount of raw gasoline vapor, the pattern of X—1 is 


best accounted for by the use of Ref and the assump- 


tion that oxidation of traffic gas had proceeded 
nearly as far as in the Los Angeles smogs. The 


analysis was eventually based on the assumption 
that the sample was two-thirds traffic gas which had 
undergone a 75 percent conversion to Ref 

Industrial impurities included the group of pol- 
lutants, in sample D-—6, attributed to Rubbertown, 
ethyl alcohol, and a definite though not large portion 


of phenol and cresol, presumably from the same 
source as the same compounds in ST-1. The amount 
of (93) unaccounted for by toluene, phenol, and 


cresol was negligibly small and much less in propor- 
tion to these compounds than in ST-1. This may 
indicate that the terpenes in ST-1 were from a source 
other than that of the cresol, although it was at first 
assumed that they represented effluents from the 
same plant. Evidence of either the numerous group 
of chlorinated aliphatic compounds or the olefins in 
ST-—1 is lacking in X-1. 

Table 7 lists the compounds identified in X-—1 
Both normal butane and normal pentane are 
undoubtedly constituents of Ref. Their inclusion 
in the accounting as additional to the usual pattern 
of Ref is an accounting device and does not neces- 
sarily indicate that they were isolated pollutants 
They may indicate some admixture of gasoline vapor 
from the same source as that in D—6 but not dis- 


tinguishable as such, and the evidence of butane 
could be almost as well interpreted as acetone 


However, a strong indication of butane occurs in 
. several Louisville samples, and it may be an indus- 
trial pollutant distinct from gasoline. The identifica- 
tions of carbon disulfide, fluorotrichloromethane 
(Freon 11), tertiary-butyl alcohol, and penty! 
benzene are based each on a single peak and may be 


TABLE 7 Inter pretation of Louisville sam ple Y-] 


Principal Peak Mole 
Substance identi height fraction 
fying i total ( 
peak pectrum 
NTQG 16 
Ref 50) 
Pentyl benzene, }{* oe a . 4 “2 
Xylene, CsH 106 291 73 
Toluene, C;Hg 74 4 55 
Benzene, CyH 7s 2.19 3.@ 
Cresol, CH yC,sH,OH 108 0. 053 0.08 
Phenol, CsH;sOH “4 ORO ll 
Styrene, CsH<¢ 104 O15 0 
Chioroprene, CyHs5C1 AS 031 06 
Vinyl chloride, C,H C1 2 +. 164 & 51 
1,3-Butadiene, C,H a] 0. 408 1.74 
Acrylonitrile, C3;H N 53 267 0. 91 
Normal pentane, C;H 72 O78 12 
Normal butane, C,H 44 2. S02 3.80 
Ethyl alcohol, C,H;sOH $1 1. 309 2.24 
Acetylene, C,H zt 5. 142 10. 70 
Carbon disulfide, CS 76 0. 125 0.13 
Fluorotrichloromethane, C FC} 103 O13 03 
Tertiary butyl! alcohol, C,HsOH 59 073 4 


in error. The first three 
used in chemical industry, and their occurrence ip 
detectable amounts seemed more probable than that 
of other compounds found to have patterns that 
could be fitted readily into the picture. No pattern 
was available for pentylbenzene, and its assignment 
was based on the following facts: (1) There was a 
small but significant peak at (148) but none were 
observed at (133) or (119); (2) mass (148) is one of 
the series of parent peaks of benzene rings with 
paraffin side chains, with which there is usually very 
little interference from other compounds; (3) from 
the parent molecular weight, there must be five side- 
chain carbons; (4) an inspection of the patterns of 


other benzene-paraffin hydrocarbons shows | that 
they tend to fragment by breaking away alkyl 
groups. The ‘vy are prone espec tially to lose the alkyl 


in the case of toluene, hydroge ‘n from the 
carbon attached to the benzene ring; (5) if a methyl 
group were attached to the ring we should expect a 
major peak at (133), if ethyl, at (119) 
of both these peaks indicates that the side chain 
must be pentyl; (6) the (91) peak which we should 
expect to be the most prominent in the case of the 
pentyl benzene is not otherwise fully accounted for 
Sensitivity was estimated from sensitivities of the 
parent peaks of other aromatic hydrocarbons. 

Peaks that are significantly over-accounted for 
again contain a disproportionate number that are 
attributable to cycloparaffins in the traffic gas 
patterns. The excess of peak (40), which includes 
argon based on oxygen in the sample, may have been 
affected by the initial liquefaction of a small amount 
of oxygen when the sample was taken. Among 
the ions not completely accounted for, (49) and (50 
are prominent. This is of interest because the same 
two peaks appear prominently in the uns accounted- 
for portions of several Louisville —. They 
suggest the presence of biacetylene, ¢ ",Hy, but this 1s 
a relatively rare compound that would not be ex- 
pected to appear in significant amount. 
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13. Other Louisville Samples 


Several other samples taken in different parts of 
fouisville during a period of clear weather with 
little wind are summarized in table 8, and the result 
of the final accounting ts presented in figures 11 to 15. 

Sample G-1, figure 11, was very largely traffic gas, 
apparently a little less than half oxidized, with the 
addition of a rather large amount of acetylene. 
Although the aromatics and alcohol were reported 
separately because their determination is inde- 
pendent of traffic gas, they were less than are usually 
found in traffic gas. The only evidence of pollution 
from Rubbertown consisted of traces of vinyl chloride 
and dichloroethylene such as might have been carried 
by an atmospheric swirl of short duration. Of the 
peaks not fully accounted for, (49) and (50) were 
again the most prominent, and cycloparaffins as 
indicated by peaks (69), (70), (83), (84), (97), (98), 
111), and (112) were again overaccounted for. 
Difluorodichloromethane was introduced into the 
sample during laboratory operations and will not be 
mentioned in the discussion of other samples in which 
it occurred; but the small amount of carbon tetra- 
chloride is believed to have been from a Louisville 
source. The sample, on the whole, was a typical 
traffic gas almost free of industrial contamination. 

Sample G—2, figure 12 was almost like G-1. The 
sample was traffic gas apparently about one-eighth 
oxidized and entirely without evidence of industrial 
pollution, if we except a small amount of tetra- 
chloroethylene and the unidentified constituent 
represented by peaks (49) and (50). To permit a 
comparison of some constituents computed directly 
from the spectrum of G—2 and independently from the 
patterns used for traffic gas, the latter have been 
marked with X at masses (78), (92), (106), and (31), 
masses which are usually significant of benzene, 
toluene, xylene, and alcohol, respectively. 

Sample G-—3, figure 13, was again traffic gas only 
slightly oxidized and free of the usual industrial 
contamination of the city, but strongly polluted with 
carbon tetrachloride and a little chloroform. Peaks 
49) and (50) were again high, though not much so, 
and (59), attributed to tertiary butyl aleohol in 
sample X—1, was again prominent 

The traffic gas of G—5, figure 14, was about 40 per- 
cent oxidized. There was evidence of industrial 
pollution and on the basis of previous samples, it was 
thought to come from at least four separate sources. 
The industrial pollutants identified were cresol, phenol, 
chloroprene, carbon tetrachloride, and dichloro- 
methane. Acetylene and alcohol were not in excess 
of the usual content amount in traffic gas. 

Sample G—6, figure 15, appeared to have a base of 
about 40 percent of completely oxidized traffic gas 
Ref), and 15 percent of gasoline vapor. It con- 
tained small amounts of the Rubbertown pollutants, 
chloroprene, butadiene, acrylonitrile, and acetylene, 
less of cresol, phenol, and pinene. There were 
methyl chloride, carbon tetrachloride, tetrachloro- 
ethylene, and some evidence of carbon disulfide. 
There was less benzene, toluene, and alcohol than 


expected from traffic vas The listed paraffins, 
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Figure 11. Total contribution of identified substances to the 
spectrum of Louisville sample G1 L. 
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TABLE 8 Interpretation o ndensed po fs f Lou 


' (3-4 


rraffie ga 
NTG ‘ 1s F 4) 
Ref 2 : 
Gasoline vapor 7 S¢ 
Paraffin hydrocarbons: t 
Normal butane ~“ be 
Isobutane 
Normal pentane 0.12 12 to 
I hexane 


: 1.3 0 

Olefins: ¢ ° x 

cia-2-Pentene f 11 re 

1-Hexene 

+Methyl-l pentene wT) 07 | i) 
Other aliphatic hydrocarbons | Ww 

Acetylene } 7 | j 1 2 1 ‘ 3.1 | 

3-Hexyne f 0. 64 h 

1,3- Butadiene ¢ { 1 7% 17 

Vinyl acetylene tit 
Benzene homologs: ¢ Q) 

Benzene 2 2 $ O2 1.2 1.8 0. 1. 23 3] 

Toluene 2 9 ( 0. 03 1¢ 6 1 22 02 N 

X viene and ethyl benzene ‘ 7-2 0 a 7 0. 43 

Above CaH j 12 13 2 p 

CyHw i" 9 { 

Above CoH 

Cul ’ | 
Other hydrocarbons 

Styrene 02 { 

Pinene s < 
Phenols - 

Phet ol w ! i xi 2? 

Cresol s4 Os 7 (2 / 
Chlorine derivatives | ( 

Methyl chloride ‘ 1.2 | 

Dichloromethane ‘ : Of 

Chloroform 2 1. 92 

Carbon tetrachloride 2 111 0.19 17 

Dichloroethane 9 51 

Dichloropropane 


13 
Vinyl chloride y. ] " 


s ; 

Dichloroethylene 29 i 2 

Tetrachloroethylene 1] 07 

Chloroprene 0. Of 114 
Difluorodichloromethane 1 a 22 22 4 10.9 
Fluorotrichloromethane 


Acrylonitrile ( a] (). 21 
Ethy!] alcohol 24 2 04 2 24 0. 24 0. 44 0.51 0.4 45 
Tertiary butyl! alcohol 


0.14 
Carbon disulfide 1 a 
Approximate total pollutant in atmosphere, pp l 17 1.0 18 
* In addition to gasoline represented by patterns of NTG and Ref 1 Laboratory contamination, not included in computed percentage 
> In addition to any of compound represented in patterns of NTG, Ref, and Acetylene from this sample lost in pretreatment of sample 
average gasoline { These values calculated from measured pressure of the total condensate 
Total of compound n sample 
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normal butane and 2-methyl-pentane, were more 
likely a portion of the gasoline vapor than relatively 
pure materials from other sources. Their appearance 
in the accounting probably results from the fact that 
the pattern of Los Angeles gasoline was derived from 
completely vaporized fuels of somewhat different 
general composition. Vapors from stored gasoline 
must contain a higher fraction of the lower-boiling 
constituents, and pentane and hexane are under- 
estimated when the regular gasoline pattern is used. 
There is some evidence of additional olefins, ac- 
counted for as isobutene and 1l-hexene. They 
may have been a part of the gasoline vapor, but they 
may also have come from the same sources that 
provided the olefins in sample ST-1 that could not be 
attributed to gasoline. Of the over-accounted-for 
peaks in this sample, (83), (98), (97), (111), (112), 
and (126) are evidence of the cycloparaffins in the 
Los Angeles gasoline pattern; and (81), (94), (95), 
and (96) are also accounted for entirely by the 
patterns of gasoline and Ref used. The spectrum 
actually obtained showed a large quantity of di- 
chlorodifluoromethane which was probably intro- 
duced in the laboratory and was, therefore, taken out 
of the spectrum in advance of the accounting 
represented by table 8 and figure 15. The Freon 
was based on the total (87), part of which must have 
come from something else, for there was not enough 
(31) and (S88) to complete the pattern, 


14. Pasadena Sample S23 


Among the many old spectra of samples from 
Southern California, one was of special interest 
because it was taken during a period of what appeared 
to be, from the vague description available, one of 
severe smog; but from a formula based on ratios of 
certain peaks (see fig. 4), the percent oxidized was 
only 49.6. Even if the concentration of pollutant 
was unusually high, a relatively slight effect would 
have been anticipated from the extent of its oxidation. 

The problem was attacked by taking out from the 
spectrum the pattern of an equal mixture of Ref and 
NTG, followed by the separation of individual com- 
pounds. It was soon evident that the peaks rela- 
tively high in NTG were being “over-accounted for.”’ 
In effect, the procedure from then on was one of 
trial and error that resulted in the interpretation 
shown in table 9 and represented by figure 16. 
Table 9 indicates relative concentration only in terms 
of the percent of total ion current occurring at the 
principal identifying peak. From this it appears 
that the sample consisted of traffic gas about 90 
percent oxidized to Ref, some aromatic hydrocarbons 
heavier than benzene, a complex group of ecyclo- 
paraffins in much greater proportion than is usual 
even in Southern California, and a number of sub- 
stances too varied and in too small quantity for 
identification. Although the evcloparaffins are rep- 
resented in the table by a half dozen compounds that 
have a combined spectrum substantially like that of 
the otherwise unaccounted for residue, it is not to be 
supposed that they the only cycloparaffins 
present Residues of (125) and (126) probably repre- 
sent cycloparaffins that could not be fitted by known 


are 
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TABLE 9. Interpretation of Pasadena smog No. S—23 


Peak | 
height, % 


Principal 


Substance identifying 


peak of total 
spectrum 
Aromatic hydrocarbons 
Cy Hy 133 0.005 
CyoHy 134 045 
CoH - 120 061 
CsH 106 331 
Cr;H¢s 92 . 912 
Cele 78 1.094 
Tetrachloroethylene, CoC], 129 0.012 
Trichloroethylene, C;HC] 130 124 
Carbo: tetrachloride, CCh. 117 . 058 
Acetylene, CoH, 25 906 
Ethanol, CeHsOH.. 31 244 
Cyclonaraffins 
+ 2 140 028 
Celis lil 326 
CeHy 112 14 
CrHis Ys 129 
Cry 97 O68 
C;H 70 839 
NTG 10 
Ref 7 


patterns without an excess of lower masses. The 
catalogs of spectra contain patterns of more than 50 
cycloparaffins in the range of molecular weight in- 
volved, and a much larger number is possible. The 
chlorine compounds, alcohol, acetylene, and benzene 
are present in about their usual proportions in the 
atmosphere of the region. 

It seems probable that the heavy aromatics and 
the cycloparaffins that distinguish this sample came 
from a common source that may have been crude oil 
or a heavy petroleum residue. Because this was a 
unique sample among many from the same location, 
it is suggested that a job of paving or road oiling 
may have been in progress in the neighborhood 
when the sample was taken. The numerous small 
unaccounted for residues probably represent slight 
pollution from varied industrial or domestic sources, 
including waste disposal. 


15. Detroit Sample No. 5 


A sample taken in Detroit during a period of 
smog is interpreted so far as seemed practicable in 
table 10 and igure 17. The plot, arranged in the 
order of descending peak heights of the sample, 
shows excellent agreement with a mixture of 60 
percent of NTG and 40 percent of Ref down to a 
peak height of about 0.1 percent of the total ion 
current. Beyond that, there are nearly 50 small 
peaks that are in major part unaccounted for. They 
indicate a very complicated mixture of pollutants, 
none of which can be present to the extent of more 
than a few tenths of 1 percent of the total condensate. 
In this, they resemble the less numerous unaccounted 
for residues of Pasadena sample S-23 previously 
discussed. In the figure the points representing the 
accounted for portion of these peaks have been 
labeled, and ions not accounted for at all are indicated 
by numbers near the bottom of the figure. 

The amount of benzene present is a little less than 
that normally found in traffic gas; toluene is nearly 
normal, but all the other aromatics much above 
normal. Alcohol and acetylene are each a little 
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Figure 16 Total contribution of ide ntified tances to the spectrum o Pasadena sample S—? 
TaBLe 10. Composition of Detroit smog No identified among the traffie gases It is hence im- 
possible to say with certainty what they represent 
Princiy \ reasonable guess can be made at some of them 
; however, on the cround that some compounds are 
more commonly used or more volatile than others 
7 lor the first reason, (58) is more likely to represent 
Aromatic I 4 . 
CoH { acetone or butane than acetaldehyde, and for both 
Cui + i reasons, (116) is more probably from an ester than 
( H 2 from an acid or alcohol Butvl acetate, of parent 
( 1M { ; . : 
CoH, ‘9 mass (116), is a widely used solvent in plastic and 
= . lacquer industries. As in Louisville, peak (111), 
Carbon tetrachloride, C4 which seems always to be the most prominent peak 
tr mh { rT . . 
sents Coil of the ecycloparaffins when heights are plotted 
— logarithmically, is more than accounted for by 
Ref 38 traffic gas alone 


less than would be expected for traffic gas suta 
diene appears from (53) and 


for traffic but compounds with six o1 


54) to be about normal 


vis, 


-bond appear in unusual quantity and are indicated 
by underscored (S1), (S82). (95). (96 LO9), 23 
(124), (137), and (138). Identifying connections 
have not been established between the many small 


peaks of high mass and their fragments have not been 


more 
carbon atoms and with two double or one triple 


16. Cincinnati Samples 


Two samples taken in Cincinnati during a period 
of high barometric pressure, were found to be unique 
among urban air samples in that they 
basically of halogen compounds instead of hydro- 
carbons. One sample is represented by table 11 and 
figure 18. Hardly more than trace 
cresol and phenol were recognized, and no more than 
the normal amounts of the individually determinable 
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constituents of traffic gas, aromatics, acetylene, and 
aleohol. As in Detroit, there numerous but 
very small unidentified peaks that indicated slight 
and extremely varied industrial pollution not closely 
enough related to suggest a common source 

In addition to the difluorodichloromethane believed 
to have been introduced at the time of analysis, the 
samples contained carbon tetrachloride, chloroform, 
tetrachloroethylene, and fluorotrichloromethane. The 
source of the halogen compounds and its relation 
to the sampling point is unknown, but it is incredible 
that halogen derivatives in such quantity can be 
widespread in the city \ point of general interest 
is that the portions ol the sample attributed to 
traffic approximate composition of 
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spectrum of Detroit sample, Smog No. 5. 


.., accounted for, 


Total contribution of identified substances to the 
S$ per trum of Cincinnati Smog Vo. 3 


Sampk 


Figure 18 


sccounted for x, accounted for by tr iffic gas 


Ref in one sample and at least half Ref in the other. 
As in Louisville and Detroit, it was evident that the 
traffic gases had a smaller cycloparaffin content than 
the traffic gases of Los Angeles and Washington. 


17. Summary and Conclusions 
The composition of the condensate, at liquid 
oxygen temperature, from any urban atmosphere is 
too complex for resolution into its individual com- 
ponents; but by considering pollutants from each of 
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various sources as substantially constant in composi- 
tion, it is possible to identify many of the sources 
and to determine approximately how much each has 
contributed to the total. This method of study has 


been applied to condensates from seven cities: 
Los Angeles, Pasadena, and Wilmington, Calif.; 


Louisville, and Cincinnati. 
California cities were numer- 
ous and enough alike to be treated, for the purpose 
of this paper, together. When minor variations 
have been allowed for, substantially all the samples 
fall readily into a mixtures of regularly 
changing composition that could be regarded either 
as blends different proportions of two basic and 
nearly uniform mixtures, or as various stages in the 
transformation of one basic mixture by chemical 
reactions into another. At one end of the range of 
composition found in Southern California is a sample 
taken from a highway tunnel in which it is probable 
that the pollutant was nearly all from moving auto- 
mobiles, and at the other end are samples, very much 
alike, taken at different places during periods of 
severe If we assume that mixtures of the 
second type are derived from the first by chemical 
reactions, it is evident that they are not end products 
of those reactions, but are the result of an approxi- 
mately steady state that is approached while the 
initial substances are regularly poured into a steadily 
reacting mixture. As a matter of convenience the 
first stage is referred to as “normal” and the second 
as “oxidized” traffic gas without implying that these 
terms tell the whole story. The term “traffic 
is intended to include engine exhaust gases and the 
amount of gasoline vaporized in the ordinary han- 
dling and use of the fuel. 
Condensates from the 
cities reveal basic pollution, 


Detroit, Washington, 
The samples from the 


series of 


smog. 


gas 


atmospheres of eastern 
very much like that of 


southern California, in which variations in the pro- 
portions between normal and oxidized traffic gas 
are clearly discernible. With these traffic gases, 


which are of substantially constant but virtually un- 
known chemical composition, a variety of other pollut- 
ants occur in relatively small amount. These minor 
pollutants usually come from local industrial sources 
Typical mass-spectral patterns have been selected 
for the initial and final of traffic and 
when these are taken out of the spectra of conden- 
sates under study, the principal pollutants not attrib- 
utable to automobile traffic are much more easily 
recognized and determined. . 

From the studies so far made by this method the 
following conclusions may be drawn: 

The pollutants in the air of Washington, D. C.. 
which has almost no chemical industry and in the 
summer no obvious source of appreciable pollution 
except traffic resemble the pollutants of 
Angeles very closely. However, the traffic 
are seldom or never in an advanced state of oxidation. 

2. As ‘orollary, we may conclude that industrial 
pollution in the Los Angeles area is minor or that it 
is of a character, a mixture of fuel vapors and prod- 
ucts of combustion, not readily distinguishable from 
traffic gas. 


Stages fas, 


ASCs, Los 


FAases 
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3. A group of mass peaks, believed to represent 
cycloparaffin compounds, is much more prominey; 
in California and Washington traffic gases than jy 
those of Louisville, Detroit, and Cincinnati. This is 
attributed to the derivation of the motor fuels on the 
coasts and in the interior from petroleums of diffep. 
ent types. With this exception, no general disting. 
tion has appeared between traffic gases in varioys 
places. 


| Inf 


4. Some compounds that occur in all traffic gases 


have identifying peaks that are relatively free of 
interference by other compounds ordinarily encoup. 
tered. They include benzene and toluene as indi- 
vidual compounds, and heavier aromatic hydroegr. 
bons as groups of the same molecular weight, acety- 
lene, and ethyl alcohol. They are determined and 
taken out of the spectrum of a sample as individual 
compounds, not as part of traffic gas. Concentrations 
much greater than normal for traffic are 
noticeable and are attributed to other sources. 

5. In contrast to the uniformity of samples from 
southern California, samples from Louisville have 
shown large variations that cannot be attributed to 
traffic gas. Most of the chemical materials respons- 
ible for the larger variations could be identified, 
Some could be identified only partly or not at all: 
but essentially the same pattern, repeated in several 
samples, indicated pollution from a common source 
which it should be possible to locate by further study, 

6. Repeated joint occurrences of the same polu- 
tants in samples from Louisville indicate that they 
originate in the same plant or nearly adjacent plants. 

From the analytical data alone, the following 
grouping appears to be probable: 

a. A source of gasoline vapor in large quantity. 

b. A source or neighboring sources of acetylene 
and vinyl chloride. 
Cc. A source or 
acrylonitrile, chloroprene, 

chloreoethylene 

d. A source of benzene. 

e. A source of cresol, phenol, and styrene. 

f. A source of pinene. 

¢. A source or sources of ethyl alcohol. 

h. A source of aliphatic chlorine derivatives in 
considerable variety 

i. A source of olefins in variety. 

j. A source of carbon tetrachloride. 

k. A source of tentatively 
tertiary butyl alcohol 


vases 


neighboring sources of butadiene, 


vinyl acetylene, and di- 


ion (59 attributed to 


1. A source of ions 149) and (50) not certainly 
identified with any compound 
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Infrared Emission Spectrum of Silicon Carbide Heating 
Elements * 


James E. Stewart! and Joseph C. Richmond 


\ method for determining the spectral emittance of a silicon carbide heating element 


Globar 
four different 
816°, and 982° ¢ 


temperatures 


900°, 1,200°, 1,500 


was devised and data were obtained over the wavelength range 1.25 to 15.25 yp for 
The spectral emittance at all four temperatures, 482°, 649°, 
, and 1,800° F 


, was found to be approximately con- 


stant except for a gradual rise in emittance from 1.25 to about 4 ~ and two minimums in the 


curves at longer wavelengths 


The first of these minimums, which occurred at about 9 gz, 


was ascribed to minor amounts of SiQ2 present on the surface of the silicon carbide heating 
element, whereas the second at about 12.5 » corresponds to one of the strong Raman lines 


for Sif 


1. Introduction 


There much current interest in heat-transfer 
processes at high temperatures. The extreme tem- 
peratures encountered in jet aircralt engines have 
long presented problems; but in modern aircraft 
performing at supersonic velocities even skin tem- 
peratures are well into the range where metals show 
Because 


Is 


undesirable creep and corrosion properties. 
radiative processes play a primary role in the trans- 
fer of heat to and from such parts it is desirable to 


obtain data on the thermal emission of various 
metals and coatings of engineering interest in a 
temperature range from 482° to 982° C (900° to 


1.800° F) 

Some work has been done in measuring the total 
emittances of various materials of interest. Total 
emittance, the ratio of the total radiant energy per 
unit area emitted by a specimen to that emitted by 
a black body at the same temperature, not : 
satisfactory measure of the amount of radiant energy 
emitted or absorbed by a specimen for use in heat- 
transfer computations because (1) the spectral dis- 
tribution of the radiant energy emitted by all sur- 
faces, including black bodies, varies with tempera- 
ture, (2) the spectral distribution of the radiant 
energy emitted by a surface may be markedly 
different from that of a black body at the same 
temperature, and (3) air and other absorb 
radiant energy at certain wavelengths that 
characteristic of the vas. 

Spectral emittance is the ratio of the radiant 
energy per unit area emitted by a surface at a given 
wavelength to that emitted by a black body at the 
same temperature and wavelength, for all wave- 
lengths at which the specimen or black body emits a 
significant amount of energy. The radiant energy 
transmitted from a hot to a cold surface can be 
computed if the temperature of both surfaces, their 
spectral emittances at their respective temperatures, 
and the spectral absorption of any material through 
which the radiant energy passes between the surfaces, 


IS 


fuses 


are 


are known 


* This project was under the spor 
National Ad ory Committee for Aeronautic 
Present addr Beckman Instruments, Inc., Fullerton, C 
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Absolute measurements of infrared radiation are 
difficult to make. The spectral distribution of 
energy emitted by a black body as a function of 
temperature has been known for many years, hence 
most emission measurements today are referred to 
a black body standard. By definition, black 
bodies absorb all radiation incident upon them, and 
emit according to the Planck radiation law. <A 
cavity with its walls at uniform temperature, having 
an appropriately small opening, furnishes a close 
approximation to the required conditions, and is not 
only the best standard for laboratory use, but also, 
under the proper conditions, is very nearly ideal in 
its physical characteristics. However, a small open- 
ing focused on a spectrometer requires that the 
spectrometer slit also be kept small, and this in turn 
requires a high amplifier gain and slow scanning 
speeds. It also makes focusing more difficult. To 
circumvent these difficulties it was decided to mea- 
sure the spectral emission of the various test surface 
with respect to a secondary standard and then to 
adjust the resulting data to emittance after measur- 
ing the secondary standard against a black body. 

A recrystallized silicon carbide rod (Globar) was 
believed suitable as the secondary standard because 
of its ruggedness, resistance to oxidation, and the 
ease with which the temperature can be controlled 
when the material is resistance heated by passing a 
current through it. 

The present paper describes the method used for 
measuring the spectral emittance of test bodies and 
for calibrating the Globar, used as a secondary 
standard, to a sufficient degree of accuracy to permit 
the comparison of emission properties of test bodies. 
Values given below for the spectral emittance of 
Globar are not to be taken as standard emittance 
data for silicon carbide for reasons to be discussed. 


2. Review of Earlier Work 


The total emittance of silicon carbide powder 
bonded with 25 percent of clay was reported by 
Pirani [1] * to be nearly 1.0 at 200° C, dropping to 
about 0.75 at 900° C. The spectral emittance of a 


as 


in bracket "indicate the literature references at the end of this paper. 











silicon carbide heating element (Globar) was first 
reported by Silverman [2] for the range 2 to 15 wu and 
for a single temperature of 1,102° C (approximately 
2,015° F). Silverman used a specially prepared 
hollow cylinder of silicon carbide with a large outside 
diameter. An axial slot cut through the wall per- 
mitted radiation from the cavity to be compared 
with radiation from the surface of the source. His 
emittance curve shows a gradual rise from about 
0.73 at 2 uw to about 0.84 at 5 wu, a fairly constant level 


to about 9.4 uw where a shallow minimum occurs 
followed by a stronger minimum at about 12.5 4g. 
He identified the minimum and maximum of the 
emittance curve at about 12.5 and 10 uw with the 


maximum and minimum, respectively, of the re- 
flectivity curve at 12 and 10 ug. 

Briigel [3] measured the spectral emittance of a 
standard spectrometer of silicon carbide 
(Globar) in the 0.6 to 15 uw region for a series of 
temperatures from 727° to 1,527° C (approximatel) 
1,340° to 2,780° F), with a comparison cavity black 
body at 727° C and instrument deflections calculated 
from those obtained at 727° C for measurements at 
other temperatures. He shows an emittance curve 
for 970° C (1,778° F) which rises from 0.75 at 2 u to 
a maximum of about 0.87 at 8 uw and falls to about 
0.77 at 15 4. His curve shows a minimum at about 
12.2 w but none near 9.4 uw, probably because his 
measurements were made at 0.5 uw intervals. Briigel 
reports a linear temperature dependence of spectral 
emittance ranging from about 3 10 per degree 
Kelvin at 2 w and 14 yu maximum of 
12 10~-° per degree Kelvin at 10 x. 

The reflectance spectrum of silicon carbide was 


source 


to a about 


first measured by Coblentz [4] and repeated by 
Schaefer and Thomas [5], who also measured the 
infrared absorption spectrum of a plate. The re 


flectance spectrum showed only the strong peak ut 
12 w and the minimum at 10 uw mentioned above 

Ramdas [6 recently redetermined the absorption 
spectrum between 1 and 20 u for a number of modifi 
cations of silicon carbide, using plates ranging be 
tween 0.14 and 1.0 mm in He found a 
very strong band probably centered near 12 u and 
numerous weaker combination Ramdas as 
signs the 12 uw band to a vibration mode involving an 
asymmetric stretching of C—Si bonds 

The Raman spectrum of silicon carbide has been 
observed by Naravanan ri who reported a strong 
line at SIS em 


thickness 


tones 


12.2 


3. Experimental Method 


A Perkin-Elmer model 21 double-beam infrared 
spectrophotometer [8] was used with a sodium 
chloride prism for all measurements. Ordinarily 


the Perkin-Elmer instrument is used for transmission 
measurements and is arranged as shown schemati 
cally in figure 1. A Nernst glower source, N, is 
viewed by two pairs of mirrors which produce 
“Images, magnified two times on the monochromator 
slits. The instrument the ratio of the 
radiance * of the two images and this ratio depends 


records 


Radiant flux per unit surface 


406 








; 
. - 
a f 
= fk | 
> u ~ | 
¢ 2 N 
oO Hf | 
z 
4 
2) ~ 
a Su, | 
| zt I 
FIGURE | Schematic drawing of 


double-beam infrared 8s pectro pho- 














tometer when arranged for trans- 
mission measurements 
N, Nernst glower source; M, mirror; 8, specimet 
~ 
| Nn; " 
2 _ 
< f 
= 8) 
2) ™~ 
x ON 
x 
= 
5 fl 
> 1 
5 us 
= IGore QM, 
| ray 
¢ 
~ 
ox 
Figure 2. Same instrument as in 
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on the absorption of a transmission specimen placed 
ats 

For measurements the equipment is 
modified as shown schematically In) figure 2 The 
mirror, MI,, is rotated as indicated so as to produce 
an unmagnified image of either the Globar souree. 
G,. or the black body, B, on the slit Because the 
path lengths of the two beams are approximately 
the same, the absorption by atmospheric water 
vapor and carbon dioxide will be nearly balanced. 


CIHISSION 


The design of the black body is shown in figure 3. 
The core is made of graphite and the two electrodes 
not shown In the figure), one at each end of the core, 
are connected to a welding transformer of variable 
power input During a test, nitrogen 1s permitted 
flow into the Vycor eylinder surrounding § the 
cavity in the graphite This prevents rapid oxida- 
tion of the graphite An asbestos baffle serves to 


Lo 


prevent the nitrogen from impinging directly on 
the graphite, which might otherwise produce a cool 
spot al this area 
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n flowing through Vycor cylinder prevents rapid oxidation of graphite 
juring determinations at elevated temperatures 


The temperature at the inside of the cavity of the 
graphite core was measured with a Chromel-Alumel 
thermocouple During operation neither — the 
thermocouple bead nor the ceramic insulator could 
be observed through the slit, thus indicating good 
temperature uniformity within the cavity. In each 
test the temperature was maintained within +1.7° C 
of that desired by manually adjusting the power 
nput 

The temperature of the silicon carbide heating 
element was measured by a Chromel-Alumel couple 
inserted into a hole that was drilled into the side of 
the element with a Cavitron [9]. The temperature 
of this element was also maintained to within 

1.7° C of that desired 

In making a determination, the silicon carbide 
element and graphite core were brought to the 
same temperature. Mirror M, in figure 2 was ad- 
isted so as to focus the image of the element on 
the slit of the monochromator and a trace repre- 
senting emission with respect to the Nernst glower 
source was obtained over the wavelength region 
25 to 15.25 u Next, mirror M, was rotated 
htly so as to fill the slit with radiation from the 
black body and a second trace Was obtained The 
ratio of the heights of the two curves at any point, 
corrected for the difference introduced by the slight 
shift in angle of mirror Ml,, was then taken as the 
emittance of the silicon carbide surface at that 
particular wavelength. In making these caleula- 
tions, the angular variation in reflectivity of the 
aluminum coating on the mirror as well as the effect 
of changing optical aberrations with angle were 


slig 


assumed negligible for the small changes in angles 
nvolved 


4. Results and Discussion 


Figure 4 isa plot of the emittance of silicon carbide 
as a function of the wavelength of emitted energy 
for 482°, 649°, 816°, and 982° C (900°, 1.200°, 1.500 
and 1.S800° F The results are in qualitative agree- 
ment with those of Silverman 2| and of Briigel [3], 
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FIGURE 4 Spectral emittance of Globar heating element at 
four temperatures. 


but the values tend to be somewhat lower than those 
reported by these two investigators. The change 
in spectral emittance with temperature is small; 
therefore, the lower values obtained in the present 
series of measurements are probably not due to 
differences in the range of temperatures in the three 
laboratories. 

The measurements were reproducible within a 
range of about +2 percent at the shorter wave- 
lengths (A<.5 w) and about +1 percent at the longer 
wavelengths. Fluctuations in the temperature of 
the Nernst glower are probably responsible for much 
of the uncertainty. Within these experimental 
limits no temperature dependence of emittance was 
found for regions sufficiently remote from the mini- 
mums near 9 and 12 uw. The minimums however, 
were definitely temperature sensitive. As seen in 
figure 4, they shift to longer wavelengths as the 
temperature is increased, and in fact appear to do 
so as linear functions of temperature: \,=8.97 
(1+4.1K10°-°7T) and Ag—12.17 (14+-5.0810~-°7). 
(7--temperature in °K.) This’ temperature de- 
pendence was investigated further by determining 
the positions of the peaks at additional temperature 
of 193°, 260°, 316°, 398°, and 1,115° C (380°, 500°, 
600°, 749°, and 2,039° F). The broadness of the 
peaks, particularly at the lower temperatures, makes 
accurate determinations of the minimums difficult. 
It is clear, however, that the 12 u« peak should occur 
at about 12.2 » at room temperature, in excellent 
agreement with the Raman observation. The other 
peak should be found at about 9.0 » at room tem- 
perature. McMahon [8] has observed a_ similar 
temperature shift of the peak observed near 9.5 py 
in the emission spectra of glasses. He attributes this 
to an increase in the Si—O bond distance as a result 
of thermal expanston. 

Beeause no band had been observed at 9 u in the 
infrared reflection or absorption spectra, or in the 
Raman spectra of silicon carbide, the appearance of 
a peak there in the emission spectra 1s puzzling 
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Both Ramdas [6] and Schaefer and Thomas [5] made | from different parts of the Globar. In this y | 
their observations on plates oriented with the basal | the existence of a steep thermal gradient along th 





plane perpendicular to the incident radiation. Hence | Globar was demonstrated, with the ends of he} Th 
vibrations perpendicular to the basal plane might | Globar perhaps as much as 10 percent cooler ad 
not be observed [10]. An attempt was made to | the middle. An averaged intensity is observed ar ay 
observe the infrared absorption spectrum of a ran- | a spectrometer having longer slits. — eal 
domly oriented sample of silicon carbide by suspend- The peak in Briigel’s curve showing the temper _ 
ing some very fine powder in a potassium iodide | ture rate of change of emittance as a fenctiol le 
disk. The differences in refractive indices (about wavelength is produced, or at least augmented. nam 
1.65 for KI and 2.6 for silicon carbide at 750 my the shift in the 9 and 12 yu peaks. The observation -~ 7 
were too great to permit observation of a good spec- | by Pirani [1] that total emittance of a mixture 7 a 
trum. A broad band was observed at about 12 u | silicon carbide and clay decreases with increasing oak 


and a Christiansen window [11] occurred at 9.5 u, | temperature is consistent with a spectral emittanes whic 
: . ; : “A “alae ‘ ey 
no doubt where the refractive indices ol silicon that is independent of or even shiehtly increasing TI 
. age ~~, ig 
carbide and potassium iodide coincide. No band | with temperature, for as the temperature is increased | + vic 
SU res 


was found at 9 x. the peak of the black body emission curve moves to’ that 
A plausible explanation of the 9 vi band is that it shorter wavelengths where the emittance of silicoy a 
arises from an impurity. The various forms of | carbide is decreasing. The clav in Pirani’s sample i 
silica are known to have strong reflection bands near | also undoubtedly affects the emission properties . ri 
9 yu, due to an asymmetric stretching of the Si () SI As the purpose in measuring the spectral emittance “ 
bonds. Furthermore, these bands seem to shift to | of the silicon carbide heating element is to be able 
longer wavelengths as temperature is increased [12] to determine the emittances of coated metal speci- T 
Silica is a likely impurity in silicon carbide because | mens, a typical example of this application is shown! top 


the material is prepared by reduction of silica with | in figure 5. The upper curve is the emission, normal etry 
carbon at high temperature and also because X-ray to the surface, from a sand-blasted stainless-steel ri 
[13,14] and electron diffraction [13] studies have | strip coated with A-418, a material consisting of ‘ ‘ 
shown that cristobalite is formed by oxidation of | chromic oxide, barium borosilicate glass. and » Be - 
silicon carbide at temperatures encountered by heat- | relative to silicon carbide at 982° C (1,800° F). The 

ing elements in ordinary use. §-Cristobalite, the | lower curve is the emittance adjusted to be relative 6. 
high temperature form, has its reflection maximum | to a black body. : 

at about 8.8 uv at 300° C [12]. Because a peak in , 

the reflection curve is influenced by both the refrac- 








tive index and the absorption coefficient, it does not a ; nee a em 
necessarily coincide with a vibrational frequency | of 
of the material. The true frequency of the vibration - , J\ | 
producing the 8.8 u reflection peak in the 6-cristo- i. Q } | en 
balite has not been measured, but by analogy with Yi y | | ter 
a-cristobalite and quartz it probably lies at a wave- - r X 1 | rat 
length longer than 8.8 u by a few tenths of a micron Pa / } } an 
In addition, the bonding material in Globar is likely O A \ thi 
to produce its own emission peaks. / R / 

No change in the cristobalite peak, or in the entire 8 00 £ of 
Globar spectrum, was detected during the period of a / Is 
time (about 50 hr) these and other measurements | 3 ,, y. ry PB | ol 
required. Nevertheless, because of the known re — 4 / ay 1) 
changes in emitting properties of Globar sources | Z [ of ‘YP en 
with time it is necessary to compare the Globar | z °° { pt 
with the black body at intervals during a series of | & ] Fal lt 
measurements 5 70 / - | gi 

The discrepancy between the results reported here | — Vi 
and the earlier observations could be due to differ- 60 f [ | t] 
ences in the composition or aging of the specimens i ; 
studied or to experimental difficulties associated : 
with temperature equilibrium and measurement sa [ 

The spectral emittance of the Globar used in this / 
study was checked in the region between 1.0 and 40 
2.5 » by R. G. Johnston and E. D. Tidwell, using a 
spectrometer and black body furnace of the Radiom- ; - r - - : - A 
etry Section of the Bureau. Agreement was obtained WAVELENGTH , MICRONS 
throughout the region only by adjusting the tempera a ae a a a eo 
ture of the Globar to obtain agreement at 2.5 Ul tance (lower curve) of a barium bores licate coat ng o7 tain! 
The instrument used had very short slits (3 mm tel Ciiatiaees Diteen ainnals , | 
which permitted the observation of energy emitted or 0 Ae F | 
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5. Conclusion 


The spectral emittance of silicon carbide seems to 
be independent of temperature in the range 482° to 
gg2° C (900° to 1,800° F) within the experimental 
oror of the measurements, except for two minimums 
sear 9 and 12.2 uw. The wavelengths of these peaks 
appear to vary nearly linearly with temperature and 
to fall at longer wavelengths at higher temperatures. 
The peak near 12.2 yp undoubtedly arises from a 
vibration of the silicon carbide itself, whereas the 
peak near 9 yu probably is contributed by cristobalite 
which is known to form when silicon carbide oxidizes. 

The observed variation in emission from different 
regions along the Globar rod leads to the conclusion 
that a Globar source is not a precise standard for 
emissivity measurements unless the emittance from 
a fixed portion of its surface is determined by com- 
parison With a black-body source 


The authors express their gratitude to R. G, Johns- 
ton. E. D. Tidwell, and E. K. Plyler of the Radiom- 
etrv Section, and to D, G. Moore and W. N. Harrison 
of the Enameled Metals Section for helpful disecus- 


sions and assistance 


6. Appendix. Emission, Emittance, and 
Emissivity 
The use of the words emission, emittance, and 


emissivity in this paper follows the recommendations 
of Worthing | 15] 

Emission is the act by which radiant 
energ\ is emitted by a bod, as a consequence of its 
temperature only. This term for the 
rate of such emission in units of energy per unit area 
and unit time It is influenced by the composition, 
thickness, and surface texture of the specimen. 

Emittance is a property of a specimen, regardless 
of its composition, thickness, or surface texture, and 
is the ratio of the emission of the specimen to that 
of a black body at the same temperature, 

Emissivity is a property of a material, and is the 
emittance of a specimen of the material having a 
polished surface, of sufficient thickness to be opaque. 
lt represents a minimum value of emittance as re- 
gards the effect of surface texture, and a maximum 
value as regards the effect of thickness (opacity) 
the specimen 


or pre CESS 


iS also used 


of 





[10] 


[11] 


[12] 


[14] 


[15] 
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Robinson Universal Buffers 
Robert G. Pike and Donald Hubbard 


\ direct comparison of the chemical reactivity at the surface with that in the bulk 
volume of buffer solutions over the range pH 7 to pH 11.8 was made, using Corning 015 


glass as the indicator 


By observing interferometrically the surface alterations of optically 


flat specimens brought about by exposures to aqueous buffer solutions under controlled con- 
ditions of time, temperature, and pH it was possible to evaluate, at least semiquantitatively, 
the difference between the chemical reactivity of the surface and the bulk volume of the 


solutions 
ion activity of pH 9.4 


For the buffer at pH 7 the surface reactivity appeared to correspond to a hydrogen 
The apparent concentration differential between the surface and 


bulk volume falls off sharply for buffers of increased pH, amounting to only 0.2 pH unit 


for the buffer at pH 11.0, and becomes nondetectable at pH 11.8. 


Unfortunately the chem- 


ical durability characteristics of the glass were such that no differentiation could be made 
between the hydrogen ion activity of the surface and the bulk volume in solutions more 


acid than pH 7 


1. Introduction 


Although direct experimental verification is often 
difficult, there is abundant theoretical evidence that 
the concentration at the surface of a solution may 
differ considerably from the bulk-volume concentra- 
tion. J. Willard Gibbs emphasized this point when 
he derived his mathematical relation between surface 
concentration and surface tension [1].' This rela- 
tion, as simplified and expanded by Wolfgang 
Ostwald is u c/RT-dy/de, where u is excess solute 


in unit surface, ¢ is concentration or activity, and 4 
is surface tension [2]. 
Inspection of this equation reveals that if the 


surface tension Y increases with the concentration 
then uw is negative and the surface concentration is 
less than the concentration of the bulk of the solu- 
tion. However, if Y decreases as the concentration 
increases, u is positive and the surface concentration 
is greater than the concentration of the bulk volume. 
Finally, if the surface tension is independent of the 
concentration, then the concentration of the solute 
in both the surface layer and the bulk of the solution 
will be the same.” 

In a previous investigation into the relation be- 
tween chemical durability of glass and the voltage 
anomalies of the glass electrode, a pronounced in- 
creased reactivity on optical glass flats by the sur- 
face of Britton-Robinson universal buffer solutions 
had been noted [5]. The present investigation was 
undertaken to study these effects systematically over 
an extended pH range. 


2. Experimental Procedure 


Instead of using surface tension as an uncertain 
and uninterpretable indicator of concentration dif- 
ference in the complex buffers used, the increased 


_—_—_—— 


Figures in brackets indicate the literature references at the end of this paper 


Gibbs’ equation has been verified experimentally [3] and the difference in 
concentration between the bulk of the liquid and the surface determined for 
different substances by sé veral methods /4 
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chemical activity at the surface compared with that 
in the interior of the solution (bulk volume) was 
studied on a glass of known durability characteris- 
ties in order to obtain direct comparison of the hy- 
drogen ion activity [H*] at the surface and in the 
bulk volume. A glass sensitive to attack in alkaline 
solutions and having a composition near the lowest 
melting eutectic of the soda-lime-silica ternary 
phase equilibrium diagram, viz., 23% of Na,O, 5% 
of CaO, and 72% of SiO, (Corning 015 glass) was 
selected as the indicator in this study [6, 7, 8]. 
Glass was chosen primarily because its use eliminated 
oxygen as a possible contributing factor in the inter- 
pretation of the results. 

Small specimens of this glass having surfaces suf- 
ficiently flat and well polished to give satisfactory 
interference bands when placed beneath a fused 
silica optical flat, were exposed by partial immersion 
in a series of buffer solutions under controlled .condi- 
tions of time, temperature, and pH. This procedure 
leaves a portion of the specimen unexposed to either 
the surface or the bulk volume of the test solution. 
The excess chemical activity of the surface (attack 
on the glass) compared with that of the bulk volume 
of the solution was obtained from the alteration of 
the interference pattern caused by exposures of 6 
hours at 80° C. 

The solutions employed were Britton-Robinson 
universal buffers, because with these a controlled 
pH range from pH 2 to pH 11.8 could be obtained 
without introducing any new ionic species [9]. Bas- 
ically, these buffers consist of a mixture of phosphoric, 
acetic, and boric acids in which the desired pH ‘is 
obtained by adding the proper amount of sodium 
hydroxide. 

The experiments were performed with a thin layer 
of liquid petrolatum on the surface of the buffers to 
minimize evaporation during the period of the test. 
This simplified procedure was accepted and adopted 
only after ascertaining that exposure in a closed 











system without an oil laver did not give appreciably 
different results. 


3. Results and Discussion 


Figure 1 shows typical examples of interference 
patterns obtained for specimens of Corning 015 after 
partial immersion in buffers of different pH values 
for 6 hours at 80° C. In each of the photographs 
the lower part of the field is that of the immersed 
portion of the specimen, which when compared with 
the unexposed portion of the field gives the surface 
alterations of the specimen brought about in the 
main bulk and at the surface of the buffer. For the 
portion of the specimen exposed to the buffer at 
pH 4 (fig. 1,A) the displacement of the fringes to 


the left indicate a swelling of about 2/10 fringes 
Such swelling is characteristic of many glasses in 
acid buffers [10, 11]. At pH 7 (fig. 1,B) and at 


pu g (fig. 1c) a perceptible swelling for the portion 
of the specimen exposed is evident in the bulk solu- 
tion, but a distinct attack or surface cut (1. e., a 
displacement of the fringes to the right) appears at 
the surface of the solution. Such surface cuts were 
noted and briefly discussed in an earlier publica- 
tion [5]. At pH 10 (fig. 1,D) the specimen exhibits 
attack for the immersed portion with additional 
attack at the buffer surface. 

Figure 2 shows the composite curve obtained when 
these and additional data from table 1 were plotted 
with surface alterations as the ordinates and pH as 
the abscissa The closed circles represent the values 
found for alteration of the specimen occuring in the 
bulk volume of the solution. The open circles give 
the accompanying attack (surface cut) attributable 
to the chemical activity of the surface of the solu- 
tion. Inspection of figure 2 reveals that these sur- 


face cuts may be interpreted directly in terms of the 
effective pH for the surface of the buffer solutions 
For example, for the buffer at pu 7 the alkalinity of 
its surface has caused an attack on the glass equal 
by a 


to the bulk-volume attack caused buffer at 
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FIGURE 2 Chemical durability characteristics of Corning 015 
glass o the ranqe pu 2 to pH 11.8 
Exposure 6 hours at 80° ¢ 4, Surface alteration of the specimens in the bulk 
me of the buffers; B, surface alteration of the specimens brought 
the increased surface reactivity of the buffet 


pH 9.4. In other words the surface of the buffer 
at pH 7 appears to be more alkaline than the bulk- 
volume value by at least 2.4 pH units 
for the buffers at pH 8 to pH 11.8, the observed 
differences are given in table 2 

At pH 11.8 the effective reactivity of the buffer 
and that of its surface either approach the same 
value, or the use of this glass as an indicator loses its 
This latter can be dis 


effectiveness suggestion 


missed, because from the ascending shape of the 
chemical durability curve it is evident that the glass 





FIGURE 1 Interferometric patterns illustrating the chemical durability characteristics of Cor ning O15 glass for exposures of 6 hours 
at SO°-C' 

A, Pattern for pH 4, illustrating the swelling generally initiated wid buffers; B, patter or pH 7 showing slight attack at the buffer surface and swelling below; 

C, pattern for pH 9 illustrating liberal attack at the buffer surface and swelling below DD, pattern for pH 10 exhibiting attack for the immersed portion wit! xiditional 


attack at the buffer surface 
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TABLI 2 Indicated surface pH of various buffers in the 
alkaline range 
Compiled from figure 2 
Bulk-volume Surface Difference 
ilkalinity ilkalinity 
pu pu pu 
ri u4 2.4 
» UH 1.6 
, 10.4 1.4 
10 10.9 oY 
11 11.2 2 
11.8 11.8 0 
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§ pet imens 
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should be a very sensitive indicator at this high 
alkalinity. The attack by the bulk volume for 6 
hours at pH 11.8 was so large, it was feared that 
any accompanying surface cut might have been 
obscured. To ascertain for sure if such a cut did or 
did not exist another specimen was given a 1-hour 
exposure at pH 11.8. Inspection of figure 3 obtained 
with the Fabry-Perot interferometric technique [12] 
reveals no indication of an excess surface attack. 
The reason for the surface and the bulk volume of the 


buffers exhibiting equal activities at high alkalinities 
may be (a) that the surface approaches saturation 
or (b) that the same molar change for the surface at 
pH 11.8 that appeared at pH 7 would be a very 
small pH change indeed at the higher concentration. 

Unfortunately no data were obtainable in the acid 
range because this glass gives a small uniform 
swelling over the range pu 2 to pu 6.5, making any 
differentiation of activity between the surface and 
bulk volume impossible. 


4. Concluding Remarks 


The difference in chemical reactivity at the surface 
and in the bulk volume of buffer solutions has been 
examined over a broad pH range using Corning 015 
glass as the indicator. This glass served as an 
acceptable indicator over the range pH 7 to pH 11.8, 
but was ineffective for buffers more acid than pH 6.5. 
The increased chemical reactivity of the surface of 
buffers may be interpreted in terms of an increase in 
the effective pH at the surface of the buffer solution. 
This difference in effective pH between the surface 
and bulk volume was greatest near pH 7 where it 
was to be expected, and decreased with increasing 
pH. 

Liquid-line corrosion has been observed on many 
substances other than glass. Examples are calcite in 
hydrochloric acid, sucrose and sodium carbonate in 
water [13], and metals in aqueous solutions [14]. 
Several theories have been advanced to account for 
the surface cuts in the cases cited, but none of these 
theories alone are compatible with the data obtained 
in the present study. For example: (1) Increased 
oxygen content at the surface of the liquid can be 
eliminated from consideration because a material was 
used that was already oxidized; (2) increased 
diffusion rate of active ions, increased reaction at 
the surface, or protection of the main body of the 
specimen by the downward streaming of reaction 
products might explain the increased surface attack 
in some systems, but for the Corning 015 glass a 
change from swelling in the body of the solution to 
attack at the surface (as shown at pH 7) can be 
accounted for only by a change in alkalinity at the 
surface of the liquid. 

In analyzing these results, it is important to point 
out that knowledge of the value of the pH of a solu- 
tion is not sufficient to predict its chemical reactivity. 
The attacking ions must be of a type that will form 
soluble compounds with the constituents of the glass 
network. For example, most glasses are inert to an 
aqueous solution of NH,OH at pH 13.3 [15]; yet 
these same glasses will be attacked rapidly by NaOH 
prepared at an equivalent pH. Thus, the results 
obtained in the present study do not apply to all 
alkaline solutions of varying pH but ay to the 
particular series of buffers used. 

From inspection of the photographs it can be seen 
that some gain in accuracy and increase in sensitivity 
of the method can be obtained by tmproving the 
quality of the polish and optical flatness of the sur- 
faces of the specimens. It was shown that further 
refinement could be brought about by increasing the 
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sensitivity of the interferometer by employing the | [7] 
Fabry-Perot [12] or the Saunders’ multiple reflection 
techniques [16]. [8] 


There is at least one other particularly Interesting 
observation from these data that not primarily 
concerned with the surface versus bulk-volume con- | 
centration of buffer solutions, but should not be 


is 


overlooked. For example, inspection of the inter- | [10] 
ference patterns of figure 1 shows at a glance why 
glass specimens partially immersed in aqueous solu- 
tions should be mechanically weakened at the liquid 
air interface. ia) 
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Photometric Determination of Tungsten in Steel and 
Titanium Alloys With Dithiol 


Lawrence A. Machlan and John L. Hague 


\ method is described for the direet photometric determination of tungsten in steels. 


The sample, 0.1 to 0.2 2 in weight 
perchloric acid mixture is added, 
diluted to volume 
treated with dithiol 
molybdenum and 
chloroform Copper 
filtration 
acid, 
of the 


certain other 


stannous chloride 


l 


in diluted sulfurie acid solution (143 


is dissolved in aqua regia, and a sulfurie-phosphoric- 
Che solution is evaporated to fumes of sulfurie acid and 
with diluted sulfurie acid 


3 An aliquot portion of the solution is 
containing sulfur dioxide, and 


potential interferences are removed by extraction with 
and other elements forming insoluble dithiolates are removed by 
The tungsten-dithiol complex is then formed, after the addition of hydrochloric 
and dithiol, and extracted with butyl acetate. 
tungsten complex is determined at approximately 6385 mg. An 


The absorbancy 
accuracy of 0.005 


pereent of tungsten, or better, is indicated in the range 0.05 to 0.50 percent of tungsten, 


and of about 0.001 percent for amounts of less than 0.05 percent of tungsten 


The applica- 


tion of the method to titanium allovs is described in an appendix section. 


1. Introduction 


Tungsten ts finding an increased usage, in the 
range of 2 percent or less, as an alloving element 
in ferrous metallurgy. Typical of some of these 
alloys are the “nondeforming” die and chisel steels 


usually carbon-chromium-tungsten alloys), age- 
hardening chromium-nickel-tungsten alloys, low- 
tungsten tool steels such as the AISI types M 1 


and M 7, and alloys finding application in the high- 
temperature field such as the 19-9DL types of 
stainless steels. Aside from alloy addition, many 
steels, through the use of scrap, contain residual 
amounts of tungsten, which may vary from a few 
thousandths to several tenths of a percent. 
Gravimetric procedures are not well adapted to 
the determination of small percentages of tungsten. 
lron and molybdenum retard the precipitation of 
small quantities of tungsten by the usual cinchonine 
procedure. Titanium, niobium, tin, and tantalum 
are often encountered, and any procedure devised 
must provide for their presence. Of the reactions 
available for the photometric determination of 
tungsten, the most commonly employed are those 
that use thiocvanate, hydroquinone, and dithiol 
as reagents. The thiocvanate reaction [5]! has 
found acceptance [11, 12] primarily for the analysis 
of tool steels, but some interference is encountered 
with molybdenum and vanadium. The hydro- 
quinone reaction has been used in the Bureau’s 
Analytical Chemistry Section for some years in 
combination with an alpha-benzoinoxime precipita- 
tion [8], but difficulties are encountered when 
present In facet, the color has been 
recommended |9| for the determination of niobium 
and tungsten in stainless steel, using photometric 
measurements at two different wavelengths. — In 
any case, titanium [10] gives a much more intense 
color and provides a serious interference. 
The dithiol (1,2-dithiol-4-methvlbenzene) reagent 
7, 13] appeared to offer advantages for the deter- 
mination of low percentages of tungsten in spite of 


—_—_——— 


' Figures in brackets 


niobium is 


it the end of this paper 


indicate the literature references 
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its obnoxious odor. Potential interferences are few, 
and most of these can be removed, with the same 
reagent, under conditions such that the tungsten 
color does not develop. This reagent has been 
applied to the determination of tungsten in titanium 
metal [18], in tantalum, titanium, and zirconium 
metals [6], and in steel [8, 14, 16, 19). 

The method presented here is the result of a study 
of the conditions necessary to make the reaction 
specific for tungsten in the presence of the large 
number of elements encountered in ferrous metal- 
lurgical analysis. A clear solution of the sample in 
a sulfuric-phosphoric acid mixture can be obtained 
by a simple fuming operation. The resulting solu- 
tion is well adapted [1, 4, 17] for the removal of 
molybdenum, the principal interfering element, with 
dithiol because considerable variations in acidity 
and temperature can be tolerated without tungsten 
reacting. Dithiol solution is added; the molyb- 
denum and other elements that form soluble com- 
plexes are extracted with chloroform. The insoluble 
complexes, primarily copper, are removed by filtra- 
tion. Hydrochloric acid, stannous chloride solution 
containing a trace of copper, and dithiol are added, 
and the tungsten-dithiol color is developed by heating 
the solution for 10 minutes on a steam bath. The 
solution is cooled, and the tungsten complex extracted 
with butyl acetate. The absorbancy of the extract 
is measured at 635 my and the tungsten determined 
by reference to a calibration curve or table. 

~The present method will not compete timewise 
with spectrographic methods. It was designed to 
provide values for tungsten, in the range 0,001 to 
0.5 percent, for the standardization of a group of 
spectrographic steel standards containing additions 
of 22 elements and the rare earths. A single deter- 
mination can be carried through in several hours; a 
group of eight determinations can be completed In 
a day 


2. Apparatus 


Colorimeter. An Evelyn type of colorimeter, 
including filters, voltage stabilizer, and galvanometer 








was used for the quantitative measurements reported 
Matched test tubes (22 by 175 mm) were used as 
absorption cells. 


Filtering Arrangement. See figure 1. 





Ficure 1 


Arrangement for filtering 


3. 


Sulfurie-Phosphoric-P¢ re hloric Acid Misture for 
Fuming. Add 100 ml of sulfuric acid ? to 250 ml of 
water, cool, add 100 ml of phosphoric acid (85% 
and 50 ml of perchloric acid (70%). Mix thoroughly 
and cool. 

Dithiol Solution. Dissolve 5 ha of sodium hydroxide 


Reagents 


it, and store in a refrigerator, Just before use, add 
2 ml of thioglycolic acid to 125 ml of the dithig 
solution, allow it to stand a few minutes, and filte 
through a close-textured paper. If a  precipitat 
forms, which does nof redissolve on a few minute 
standing, too much thioglycolie acid has been added 
Prepare a fresh portion, using slightly less thioely. 
colic acid. Use the turbid solution for the develop. 
ment of the tungsten color. 

Stannous Chloride Solution. Dissolve 100 g of 
stannous chloride dihydrate and 5 mg of copper 
chloride dihydrate in 400 ml of diluted hvdro. 
chloric acid (1 1) by heating on a steam bath 
Cool and dilute to 500 ml with diluted hydrochloric 
acid (1 1). 

Sulfuric-Sulfurous Acid Mizture. Saturate diluted 
sulfurie acid (1+-3) with sulfur dioxide gas. 

Standard Tungsten Solution (1 ml<1.00 mq of 
tungsten). Transfer 0.1794 g of sodium tungstate 
dihydrate to a 100-ml volumetric flask, dilute to 
the mark with water, and mix. The sodium tung. 
state was purified and tested by the procedure 
suggested by Yagoda and Fales [20], though bigh- 
purity material is not required. 


Standard Tungsten Solution (1 ml<>100 gq of 
tungsten). Transfer exactly 10 ml of standard 
tungsten solution (1 ml<&1.00 mg of tungsten 


to a 100-ml volumetric flask, dilute to the mark with 
water and mix well. This solution should be 
prepared as needed. 


4. Procedure 


Transfer 0.200 ‘ of the sample to a 
Erlenmeyer flask, add 10 ml of aqua regia (3 parts 
of hydrochloric acid and 1 part of nitrie acid), and 
heat gently until the sample dissolves. Add 20 ml 
of the “fuming acid’? mixture and | or 2 drops of 
hydrofluoric acid, boil until brown fumes have been 
expelled, evaporate to fumes, and fume the solution 
vigorously for a few minutes to remove the perchloric 
acid. Cool, add 25 ml of diluted sulfuric acid (1+3), 
and heat on a steam bath to dissolve insoluble salts 
Cool, transfer the solution to a 100-m1° volumetri 
flask, and dilute to the mark with diluted sulfuric 
acid (14-3), 

Transfer 7.5 ml of diluted sulfuric acid (1+1 
5 ml. of sulfuric-sulfurous acid solution, and a suit- 
able aliquot ° of the sample to a 60-ml separatory 
funnel. Place the funnel in a water bath maintained 
at 30 10° C. add 5 ml of dithiol solution, shake 
the solution for 5 to 10 seconds. and allow it to 
stand for 5 minutes. Cool, add 10 ml of chloroform, 
shake it for 40 seconds, and allow the layers to 
separate. Drain and discard the chloroform layer. 
Repeat the extraction with 5 ml of chloroform, 
finally discarding the chloroform layer. Filter the 


co 
= 
zie 


to 


125-ml | 


in 50 ml of water and dilute to 500 ml. Add 1¢of | acid laver through glass wool (see fig. 1) into a 250-ml 


3,4-dimercaptotoluene (dithiol),’ shake to dissolve | separatory funnel. Wash the 60-ml_ separatory 
?If no diluuwn is specified, the concentrated analytical reagent is " ‘ For materials containing more than 0.5 percent of tungsten, use 0.100 g of the 

Diluted hydrochloric acid (141) denotes 1 volume of concentrated hydrochlori sam ple 

acid, sp gr 1.18, diluted with 1 volume of water Diluted sulfuric acid (14-3 For materials containing from 1 to 2 percent of tungsten, transfer to a 200-m 

denotes 1 volume of concentrated sulfuric acid, sp gr 1.84, diluted wit! volume volumetric flask 

of water A suitable aliquot will contain up te 20 mg of sample and less than 50 pg { 
+A satisfactory reagent of domestic manufacture is commercia iVailable tungsten 
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funnel twice with 2- to 3-ml portions of diluted 


euifuric acid (14-3). 

Place the 250-ml separatory funnel in a steam bath 
for 10 minutes, cool in a water bath maintained at 
ai) >to 40° C, add 2.5 ml of dithiol solution, and shake 
for 5 to 10 seo Return the funnel to the water 
bath and allow it to stand for 5 minutes. Cool, add 
5 ml of chloroform, shake for 40 seconds, allow the 
layers to separate, drain and discard the chloroform 
laver. Repeat the extraction with another 5 ml of 
chloroform, finally discarding the chloroform layer. 
Place the funnel in a steam bath and heat for 15 
minutes.’ 

Add 50 ml of hydrochloric acid and 5 ml of stannous 
chloride solution to the funnel, return it to the steam 
bath. and heat for 3 minutes. Add 15 ml of dithiol 
solution, mix well, and again heat on the steam bath 
for 10 minutes. Cool the solution in a cold-water 
bath, add exactly 20 ml of butyl acetate from a pipet, 
shake for 20 seconds, and allow the layers to separate. 
Drain and discard the acid layer. Transfer the butyl 
acetate solution to a dry test tube and stopper the 
tube. 

Transfer a suitable portion of the blank, or butyl 
acetate, to a 2-cm absorption tube or cell and adjust 
the photometer to the initial setting, using a light 
band centered at approximately mu. While 
maintaining this photometer adjustment, take the 
photometric readings of the butyl acetate extract.° 
Determine the percentage of tungsten from a cali- 
bration curve or table prepared as follows: Transfer 
0.200 ¢g of iron cont: aining negligible quantities of 
tungsten and molybdenum (NBS Ingot Iron Stand- 
ard Sample 55d was used) to a 125-ml Erlenmeyer 
flask, and continue as directed in the procedure for 
dissolving a sample. Add 10 ml of sodium tungstate 
solution (1 ml<100 ug of tungsten) to the flask just 
before adding the 20 ml “fuming acid’? mixture, 
and continue with the procedure until the solution is 
diluted to volume in a 100-ml volumetric flask (1 
ml10 we of tungsten). Prepare a blank iron solu- 
tion by the same procedure, omitting the 10-ml ad- 
dition of sodium tungstate solution. Add 25 ml of 
diluted sulfuric acid (1+-3) to five 250-ml separatory 
funnels. Transfer 1.0, 2.0, 3.0, and 5.0 ml of 
iron-tungsten solution (1 ml<10 wg of W), and 5.0 
$.0, 3.0, 2.0, and none ml, respectively, of blank iron 
solution ° the funnels. Heat the solution on a 
steam bath, and continue with the development of 
the tungsten color, beginning with: “add 50 ml of 
hydrochloric acid and 5 ml of stannous chloride solu- 
tion.” Plot the photometric readings of the calibra- 
tion solutions with respect to micrograms of tung- 
sten. Beer’s Law is obeyed under the conditions 
specified, and blanks are usually negligible. 


635 


none, 


to 


? The 15-minute heating period is requ to remove the sulfur dioxi de, which 
inhibits the development of the tungste« ind to remove the chloroform, 
which would increase the volum«e nt used to extract the tungsten- 
dithiol complex, 

* The absorbancy of this extract varies with temperature, about 0.2 percent per 
degree; measurements should be made in a range of 5° C if an accuracy of 
1 percent is required, 

* The addition of blank solution is made to provide an adequate amount of iron 
for the complete development tur ten color 
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5. Discussion and Results 


The results obtained by the recommended pro- 
cedure on a number of synthetic mixtures are given 
in table 1, well as a few values obtained on 
mixtures containing copper, molybdenum, and tung- 
sten, which were not filtered. Copper in 200-yg 
amounts does not interfere, as is shown in exper- 
ment 1. Molybdenum in the absence of copper is 
removed by the chloroform extraction of the dithiol- 
treated solution and does not interfere, as shown in 
experiments 2 and In mixtures containing both 
elements in appreciable quantities, the black pre- 
cipitate formed by copper on the addition of dithiol 
appears to occlude molybdenum, and the extraction 
step does not completely remove molybdenum. The 
result is that on solution of the black copper precip- 
itate during the development of the tungsten color, 
molybdenum is also present and causes a positive 
error, as shown in experiments 4 and 5. A filtration 
step was therefore added to the procedure. 


as 


TABLE 1. Results of determinations for tungsten by the recom- 
mended procedure in various synthetic mixtures 


Molyb- Tungsten 
Experi Iron denum Copper _ 
ment vdded added added 
Added Found Difference 
mg ud ud ud ud wd 

1" 10 200 0.1 +0. 1 
2a 10 30 20.0 2.0 0 
3a a) - 1) iO. 0 49.5 5 
ja 10 MO) 200 6.4 +i. 4 
5a 10 30 200 20.0 23. 5 +3.5 
‘ ~) 2.0 14.4 0.6 
7 l Mw 2.0 19.8 2 
x 5 MM) 2.0 19.8 2 
) v1) “) 2.0 19.9 l 
10 ou 0) 20.0 2.2 +2 
11 i. 1) ”) wD. 0 274.5 +4.5 
12 1OO ry 20.0 $3. 1 +13. 1 
13 100 W.U 1u.8 0.2 
14 10 4) 200 0.5 +.5 
1 10 4) 1.0 1.0 0 
1¢ Ww “) 2.0 1.9 | 
17 ‘AY r aK 5.0 5.1 1 
18 x) ”) 200) 20.0 2.0 0 
10 a») Ww) 4) 2.0 19.7 $ 
m) 24 100 10 Ww. 0 20.0 0 
a4 | v1) “) x) 0.0 19.8 2 
») x) “) Ww). 0 10.5 -.5 


» Not filtered to remove copper precipitate, 


Iron also inhibits the removal of molybdenum by 
the dithiol extraction step, and the extent of this 
interference is shown by experiments 7 to 12. About 
30 mg of iron can be tolerated satisfactorily in the 
presence of molybdenum, or combinations of mo- 
lybdenum and copper. The high sensitivity of the 
tungsten-dithiol reaction allows small amounts 
(0.005 percent or more) of tungsten to be determined 
quite satisfactorily, as shown by experiments 14 to 
22. Experiment 6 shows the necessity for the 
presence of a small quantity of iron; a low recovery 
is obtained in the absence of iron [2]. As is shown 


by experiment 13, the accuracy of the tungsten 
determination is satisfactory in the presence of 
considerable iron if molybdenum is absent. A few 


experiments on 100-mg samples of NBS Ingot Iron 





Standard Sample 55d indicated that an additional 
quantity (7 to8 ml) of sulfuric-sulfurous acid solution 
and of thioglycolic acid (0.2 ml) helped in the 
removal of small amounts of molybdenum, and that 
this standard, used in preparing calibration curves 
contains tungsten in amounts of the order of 1 or 2 
ppm. A few experiments, in which copper and 
molybdenum were largely removed by treatment 
with hydrogen sulfide and filtration, indicated values 
of the same order. 

The reaction of molybdenum with dithiol also 
appears to be sensitive to temperature. Although 
this particular point was not studied extensively, 
the reaction is too slow in the presence of iron at 
10° to 15° C to be useful. The reaction temperature 
of 30° to 40° C was selected on the basis of conven- 
ience and a reasonable time for the removal of 
molybdenum. 

Because the only difficulty encountered was the 
occasional incomplete extraction of molybdenun 
at the proper point (a complaint common to many 
of the published procedures using dithiol), photo- 
metric readings were also taken at approximately 
720 my to detect a failure of this sort. A change in 
ratio of the absorbancy at 635 my to that at 720 mu 
will indicate the presence of molybdenum; in fact, 
a procedure has been suggested [15] for the simul- 
taneous determination of molybdenum and tungsten 
by photometric measurements at two different wave- 
lengths. Dithiol and thioglycolic acid were ordered 
in small quantities and stored in a refrigerator 
Difficulties were not encountered with prepared 
solutions of these reagents if they were stored in the 
refrigerator and used within a few days 

During a study of the thiocyanate photometric 
method for tungsten it was found that a relatively 
rapid and complete reduction of tungsten by stannous 
chloride in strong hydrochloric acid solution occurred 
in the presence of a trace of copper. Because this 
reduction is also required in the dithiol procedure, 
a similar technique was applied. Variations in th 
recovery of tungsten with several time intervals for 
color development, and with several amounts of 
stannous chloride and dithiol solutions are illustrated 
in figure 2. The conditions selected in the recom- 
mended procedure are about the minimum require¢ 
to give a satisfactory development of the tungsten- 
dithiol complex in the range 1 to 50 ug of tungsten 
The volume at the beginning of the color-develop- 
ment step is usually about 30 ml, but in handling 
small percentages of tungsten requiring larger sample 
aliquots, a larger volume is usually obtained 
Volumes of up to 75 ml, containing 20 ug of tungsten, 
have been satisfactorily handled, as shown in figure 2 

A few experiments in which a second extractior 
with butyl acetate was made, indicate the tungsten- 
dithiol complex is completely (more than 99 percent 
extracted in one treatment. The results of experi- 
- ments in which extra dithiol is added to an extracted 
solution, followed by reheating, cooling, and extrae- 
tion of any additional complex with butyl acetate, 
indicate that the reaction is essentially complete in 
the 10-minute period specified in the procedure \ 
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FIGURE 2 Effect of variables on the de elopme nt of the tungster 


dithiol complex 


small amount (less than 1 percent) of the tungsten 
present is lost during the chloroform extraction of 
molybdenum. This loss was checked by the evapo- 
ration of several combined chloroform extracts with 
sulfuric acid and hydrogen peroxide, followed by a 
determination of the tungsten present in the resulting 
solution, as described in the recommended procedure, 

The results obtained on a few of the National 
Bureau of Standards Standard Samples of steel are 
shown in table 2. The values for NBS Standard 
Samples 159, 123, 160, and 101d have not been 
prev iously determined, because a satisfactory method 
for the determination of such small quantities of 
tungsten in these types of alloys has not been avail- 
The values for NBS Standard Samples 123b 
155, 153, 134, 134a and the NBS Spectrographie 
Steel Standards 442, 443, and 444 show satisfactory 
agreement, considering the complexity of the ma- 
terials and of the procedures used in the original 
analyses. The higher values given for the NBS 
Spectrographic Steel Standard Samples 446 and 846 


able. 
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TABLE 2 Results of deter 





minations for tungsten by the recommended procedure in various National Bureau of 


Standards Standard Sam ple 8 


rungsten 


Number | Weight 
N BS Standard Sampk of deter of sam pk Iype of steel 
Certificate Found Difference Range minations } 
value 
mag 
159 0.005 0.005 2 1) Cr 1.0, Mo 0.4, Ag 0.1 
23 005 0.005 to 0.006 3 2 Cr 18, Nill, Nb 0.4 
160 oll 0.011 2 Ww Cr 19, Ni9, Mo3 
101d O16 0.015 to 0.016 3 a Cr 19, Ni9 
442 0.08 O78 0.078 2 Ww 
Stainless-steel spectrographic stand- 
443 09 ORS 0.084 to 0.087 3 20 ards, 
444 18 167 0.01 0.165 to 0.169 Y 10 to 20 
123b Is 174 Ol 0.173 to 0.175 4 20 Cr 19, Nill, Nb O.8&, Ta 0.2 
155 517 19 +002 0.516 to 0.521 4 ito 10 Cr 0.5, W 0.5 
153 1.58 1.55 03 1.54 to 1.55 6 2.5 Mo 8, Co8, Cr 4, V 2 
134 1.82 1.87 + 08 1.87 2 2 Mo 7, Cr4, V 2 
134a 2.01 2.05 O41 2.04 to 2.06 t\ 2.5 Mo8&, Cr4, V1 
446 and 846 0.043 0.043 to 0.044 6 10 to 20 
446 and S84 O48 0.044 to 0.050 8 10 Cr 18, Ni9, Mo 0.4, Cu 0.2 
146 and 846 O55 0.052 to 0.0459 s yt) 
» Not filtered to remove copper precipitat 


where the copper precipitate is not removed demon- 
strate the range of error that will be encountered, 
and the necessity for the steps in the procedure that 
allow for more than residual amounts of copper and 
molybdenum. The ranges given in table 2 show the 
excellent precision of the method 

Table 3 shows tests made with 15 elements men- 
tioned as reacting with dithiol or otherwise inter- 
fering in the detection of tungsten, As will be seen, 
none of these elements in the amounts tested inhibit 
the removal of molybdenum or hinder the develop- 
ment of the tungsten color, and quite satisfactory 
recoveries of tungsten are obtained. Satisfactory 
determinations were also obtained in the presence of 
500 mg of citric or tartaric acid 


TABLE 3 Results of determinations for t ingsten hy the recom- 
mended procedure in the presence oT; possible interfering 
elements 

Molyt rungsten 
Exper Other el Iron jenum 
ment ment added vided vided 
Added | Found Difference 
pe ng b uM mi) Mg 

l 200 Sb 10 4) 20.0 14.49 0.1 
2 20) As 10 4 w.0 1¥.8 2 
; 200 Bi Ww 4) 2.0 14.9 1 
4 5,000 Cr 1 A) 20.0 19.8 2 
f 2 000 Co 10 4) 21.0 Ha). 2 2 

2.000 Co 10 4 7.0 19.9 1 
7 200 Pb 1 4 20. ( 19.8 2 
8 20) Min 10 ‘ a0 19.9 l 
4 200 He lt 4 20.0 1v8 2 
10 1400 N 10 a 20.0 1044 l 
1] 600 Pt 10 4 ”.0 19.8 2 
12 200) Re lt 4) 20.0 1u.8 2 
13 200 Se 10 ».0 19 9 1 
14 200 Ag 12 ‘ 71.0 19.9 1 
15 2,000 St 10 ‘ HN) ( %. 1 l 
lt 200) NV Ww q 25.0 24.4% l 
rly . 
The reaction of some of these elements with 


dithiol may be of interest. The rhenium compound 
with dithiol is visually indistinguishable from that 
of molybdenum sismuth forms a brown precipi- 
chloroform; tin behaves 
Silver 


tate, which is extracted yy 
similarly except that the precipitate is red 
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forms a black precipitate and, like copper, is removed 


by filtration. 
a brown or pink precipitate. 


Selenium behaves similarly, forming 
Platinum added as 


chloroplatinate divides, a part being removed by 


filtering the 


brown precipitate formed. 


The 


re- 


mainder reacts on the addition of stannous chloride 
during the development of the tungsten color to 


form 


acetate, 
the extraneous color does 


not 


1 yellow complex, which is extracted by butyl 
An “off-colored”’ solution is obtained, 
interfere when 


but 
the 


absorbancy of the solution is measured at 635my as 


directed. 


Cobalt forms a blue complex on the addi- 


tion of hydrochloric acid, but in the amount used in 


the test does not appear to interfere. 


The remaining 


elements did not give visible reactions in the concen- 
trations tested, 
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7. Appendix 


7.1. Determination of Tungsten in Titanium 

At the request of the Watertown Arsenal Task 
Force for the determination of tungsten in titanium 
metal, tungsten was determined by the dithiol pro- 
cedure in their samples WA-22 and WA-8. Because 
titanium phosphate in quantity is not very soluble 
in sulfuric-phosphoric acid mixtures, the procedure 
was modified to provide a suitable solution for 
titanium-base alloys. 


7.2. Modified Procedure 


Transfer 0.2 ¢ of the sample to a 125-ml Erlen- 
mever flask, add 20 ml of diluted sulfuric acid (1-+-3) 
and 2 or 3 drops of hydrofluoric acid. Heat gently 
until the sample is dissolved and add, drop by drop 
(about 6 to 8 drops are required) hydrogen peroxide 
(30 percent) until the purple color of reduced ti- 
tanium is discharged and a permanent vellow color 
of peroxidized titanium is formed. Evaporate the 
solution to fumes, and cool. Add 50 ml of diluted 
sulfuric acid (1+-3) containing a gram of citric acid, 
and digest a few minutes to obtain a clear solution. 
Cool, transfer the solution to a 100-ml volumetric 
flask, and dilute to the mark with diluted sulfuric 
acid (1+-3). 

Transfer a suitable aliquot, as indicated in table 
4, of the prepared sample solution to a 250-ml sepa- 
ratory funnel. Add 0.2 ml of ferric sulfate solution 
(10 mg of iron per ml of diluted sulfuric acid (1+-9 
7.5 ml of diluted sulfuric acid (1+-1), and 5 ml of 
dilute sulfuric acid (1+ 3) saturated with sulfur diox- 
ide, and complete the determination as described for 
steel. Because most titanium allovs contain little 
copper or other similarly reacting elements in suffi- 


cient amount to require a filtration, this step wy | 
omitted on these samples. 


on . . 
TABLE 4 Results of determinations for fungsten in titan | 
. 


materials 


rungsten Molvyvt | 
Sampk denum Iron Titaniy 


hase 


idded vdded 


1 10. ¢ wm 7 Lo 
9 10.0 10 4 om 
; 10.0 10.0 0 om 
{ Wd. 19.7 7 
0 40 5 0) | 
40 ss. 2 s ] 
- ; 
{ 0. 372 ’ 
| Os 
WA-S - I 
} $74 « ] 
71 2 
f 7 o 
WA-22. rons a 
) 056 on 
| ve ” | 


The results obtained on a few synthetic mixtures 
and on two samples, are given in table 4. Titaniun 
was added to the six synthetic mixtures as portions 
of a titanium solution prepared from National Bureay 
of Standards Standard Sample 173, a 6-percent 
aluminum, 4-percent vanadium titanium-base alloy | 
containing approximately 0.001 percent of tungsten, | 
The values obtained by the Watertown Task Force 
by a double a-benzoinoxime precipitation and hydro- 
quinone photometric method, average about 0.40 
percent of tungsten for sample WA-S and approxi- 
mately 0.05 to 0.06 percent for sample WA-22 
Greenberg [6] has reported 0.38 and 0.055 percent 
of tungsten. for the same samples by a dithiol pro- 
cedure, 

12, 1957. 
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| On the Most General Form of the Compatibility 
Equations and the Conditions of Integrability of 


| Tt | Strain Rate and Strain 
_|f Edmund H. Brown 


The most general forms of the equations of compatibility are derived by a very simple 
analysis in terms of the material strain rate tensor and the rate of change of the Riemann- 
Christoffel curvature tensor during the deformations. It is shown that recent statements 
on the significance of the compatibility equations, essentially that their satisfaction is 
equivalent to the condition that the space be locally Euclidean, are misinterpretations based 

on restricted forms of the equations. In addition, a new form of the conditions of inte- 
grability of strain is derived in terms of the rate of change of the curvature tensor. 


tin lilaning | 


te 


~ 


0 1. Definitions and Notations 

2 Because of the great variety of usage of authors in the field, a few of the definitions and 

2 notations employed here will be discussed briefly. However, a general knowledge of tensor 
| analysis and modern elasticity is assumed. For the most general, and at the same time most 


- mixtu readable introduction to tensor analysis, including applications to elasticity, the reader is 
xtures : . 


Titaniun referred to Brillouin [1];' for an excellent and thorough review of modern elasticity through 
S portions | 1953, despite a frequently cumbersome notation, the paper of Truesdell [2] is highly recom- 
ial Bureay mended; for more recent general discussions, the treatise of Green and Zerna [3], or the article 
-percent of Doyle and Ericksen [4] should be consulted. 


base alloy { 
tungsten , 
isk Force 


To introduce the strain tensor, we imagine a complete n-dimensional material continuum 
G which, with an arbitrarily chosen coordinate system z'‘, is defined by the form of the metric 


nd hydro- tensor g,,=49:,(2*). To show the covariant character of some of the quantities involved it 
pout 0.40 will occasionally be necessary to transform from one coordinate system to another. Following 
) ry Craig [5] we shall indicate quantities in one coordinate system by use of indices from the 
- abet group 7, j, k, h, l, and quantities in a second coordinate system by indices from the group 
thiol pro- Pp, 4, 7, 8, t. Thus, such a transformation will be defined by an analytic relation of the form 
g’=az'(x*). Further, for convenience, the components of the transformation matrices Oz*/dz’ 

or 0x7/dr* will be denoted by V'=0z"/d2r’ or Xi=02"/dzr*. Similarly, 0°x*/dxr*dz’ will be denoted 

by Y’,. During such a coordinate transformation, then, in one physical state, the components 


of the metric tensor will transform according to the tensor law, g,,=9,X/X{, the Christoffel 
symbols, or connections, will transform according to the Christoffel law, P',.=TI',,X$AX{X;+ 
X3,Xi, and the components of the Riemann-Christoffel curvature tensor by the tensor law, 
R votre = RijnXiXLXEX". Note that we use a short vertical stroke between indices to emphasize 
various symmetry properties. In addition, following Green and Zerna, we use a full vertical 
stroke behind a quantity, as in the identity g;,;|,=0, to indicate the covariant derivative. 

A material coordinate system is defined by assigning a value z‘ to each particle of the 
material continuum, which it retains during any subsequent deformation. We shall assume 
that the deformation, itself, depends on a single scalar parameter yu, independent of the 
coordinates z*, which in dynamic problems can be identified with the time ¢. Since the mate- 
rial coordinates of a particle do not change during a deformation, we must have 2*(u.)=2*(u). 
However, during the progress of the deformation the intrinsic quantities associated with 
the space, such as the metric tensor g;;, must change, assuming a continuous set of values 


9ij(i1), Guy(me), . . . . When we wish to indicate several stages during a deformation, the 
appropriate values of 4, mw, . . . Will be used, but more often only two stages, an ‘‘initial’’ 
state and a “‘final’’ state, are of interest, so that, for instance, we can write for the “‘initial”’ 
metric tensor, g,;, and for the “‘final’”’ metric tensor g,;.__ In such cases, we shall make a general 


practice of indicating a final state by placing bars over quantities which in the initial state 
are unbarred. 


Figures in brackets indicate the literature references at the end of this paper. 
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At this time, it is well to note that, parallel to the ‘“‘material’”’ representation of quantities, 
a “spatial” treatment can be used: A complete n-dimensional geometric, or spatial, continuum 
H can be specified, using an arbitrarily chosen fixed spatial coordinate system y', by the form 
of the spatial metric tensor, h,,=h,,;(y*). We imagine that the deformation of G takes place 
in /7, that is, that there is a one-to-one correspondence between the “‘particles’”’ of G and the 
“points” of 7 which, however, varies in form with uw. Thus, there is a “mapping” of G on H 
by y'=y'(z* uw), whose variation with u represents the deformation of the material continuum. 
Since during the deformation, that is, during the variation of u, the particles of G move through 
the points of FT, y‘(2* yu.) #y'(2*u,), but at each given point y‘, spatial quantities such as A,, and 
the spatial Riemann-Christoffel curvature tensor S;,, remain constant. There is nothing 
of course to prevent one assuming that the space /7 may change with parameters other than ug, 
a device which could possibly be of some use in discussions of plasticity; but for simplicity 
we shall avoid this here. In any case, a necessary condition that G remain in // is 


P..—s__.. dy? dy dy” dy’ 
— ee” Oxt Ox? Ox* Ox" 
at each value of u. 

Quantities such as strain can be related to either the material or spatial systems, whereas 
the displacement wu‘ is essentially spatial. Since we intend to discuss strain only, for con- 
venience we will use chiefly a material representation and make little further reference to the 
spatial treatment. We should remark, however, that one of the principal applications of the 
results for non-Euclidean continua is in cases where (unlike our definitions above) the dimen- 
sion of G is m, the dimension of /7 is n, and m n, since they are useful in reduced problems 
involving shells and curved rods. These problems, however, require the introduction of 
n—m other mth-order quantities (‘changes in curvature’, ‘‘changes in torsion’, ete.) in 
addition to the strain, so that we shall not consider such possibilities here. 


2. Introduction of Strain and Derivation of Results 
We shall take as a measure of finite strain the quantity 


12) ¢ (a) = M[gij(x*, m2) —Giy(2* m1) ], (1) 


or, if we wish to think only of some particular initial and final states, 
é, (2 lg (x q (ar* |, (2) 


since—even though it is less used in the modern theory of finite elasticity—it reduces to the 
familiar strain tensor of ‘‘classical’’ theory, such as is found in Love [6], in the infinitesimal case. 
According to the above definition, strain depends merely on the initial and final states, and is 
independent of the way in which the deformation took place. Such a definition is too general for 
strain to be uniquely related to stress in more complex cases involving permanent set or plastic 
behavior. Thus, we are led to the introduction of the rate of strain tensor é,; defined by 

, O06; 5 I O4;, 

Ou 2 On si 
We note that, since differentiation with respect to the scalar parameter yu does not affect the 
character of q,,, é,, like gy is a tensor in the z* coordinate system for a given value of wu. Sim- 


2e.,is a tensor. However, though 


ilarly, since both g,, and g,,; are tensors, the finite strain 
z*—r* for all uw, the identities (based on the hypothesis of a Riemann space, that is, that a 
unique measure exists) 9;,|;=0 and g;,;|.=0 do not allow us to conclude either that g;;|;=0, or 
Gy\x=0, since g,,; is not the metric tensor of G, nor Giz the metric tensor of G. In G, gy is merely 
a tensor field and not the metric tensor; similar statements can be made about the character of 


dy in G, 
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Now, considering é;; to be a tensor field in G, in accord with the usual requirements for 


raising and lowering indices, we must have 


$14 — gik ging l og” 
€ ( ( Crp —_ 
g's 29 
2 Ou 
and 
L(x OM oy” 
¢ je (g a ) 
4 On * Ou 


(4) 


(5) 


Just as the metric tensor g,; changes during the deformation, the Christoffel connection 


I based on Jiy Must change. We shall define the rate of change Ee ». of er by 


Or’ ;, 
Ou 


However, since [ 
by the Christoffel law, and none of the quantities 2*, 7, YX}, are functions of y, 


fir, Ob OV yixey: 


= ~ Bp X3XpA%; 


which shows that—unlike I 
first pointed out by Bergmann [7]. 


Just as I‘, can be expressed, for a Riemann space, in terms of the metric tensor, 


l Oo”) Odn. O 
I ,, git dn 4 =... In 
2° Or or Or" 
, can be written in terms of the strain rate é,,; for 
7 OM, 10g". Oda; , Odm Og L 4, O | Odny , One Oy 
k =s eT —— se © ——- — ae & 
Ou 2 Ou Lou Or Ou 2 Ou L Or Ox Or 
or, since g;,; is an analytic function of z' and yu, and z‘ and uw are independent, 
, ‘ O€, nk O€y) 
B= —2"Ty etl oe 
2 Or Or Ox" 
can be shown to be given by 


"a ‘ : O€ O€; O€ 
’ ( gg 2é al je J { 
Ek Jinks" k ~€ 1 +S T Or | 


or 
Ou 


so the E 


Similarly, # 


or 


Using the expression for the covariant derivative of the second order tensor é,,, the 


partial derivatives of € appearing above may be written 


O€» 


7 , 
FS) Fr €h é | , é I : 
I 
OF ni a 
2 Crk f I A €ni | 
Ir 
O€ ») 7 
or’ —~€ 4k |) €xl é,,1 4 
ys . 


(6) 


transforms, when the coordinates are changed at one physical state, 


(7) 


the quantities £*,, and thus Ey ,.=—gm"», are tensors, a fact 


(10) 


(11) 


three 








On adding these expressions, interchanging indices where symmetry permits, we obtain 


Oln, , O€n OF 


— = ¢ r r 2é,,T 13) 
Ox" t Ou Ou . 
Introducing this in the expressions for £',, and £,,, in terms of the strain rate, and making 
use of the fact that q” 0 results in the more elegant forms, 
Boy x= isle tena s—E vals 14) 
and 
Ee =e" é;' é ix\'. 15) 


We must, finally, consider how the Riemann-Christoffel curvature tensor R',,,, based on 


gi Varies during the deformation. Defining 4‘ ;,,,—Of',,, Ou, this rate of change of curva- 
ture is given by 


m) (20s). O (2t:#) Orn | pe OFS OF! x | pe OF'y 
ON Oa* \. Op Or" \ Ou Ou Ou Ou ; Ou 
OE’, oF 
Sat PE aD ia —El al a — Eig l ie ( So FE aD Ela EV ) 16 
or 
¢ li } E b l4 
and similarly, 
Ps kh ir€ - sen a —2¢,,R' kh KE, Aik E, kihs LS) 


both of which are especially elegant forms for the rate of change of the curvature tensor 

Just as in the case of E, », We Can express & », In terms of the strain rate, and, in doing so. 
we obtain the equations of compatibility in their most general form. By substitution and 
simple algebra these identities may be written 


Oifith — \Cijgi\ne— Cis en) T \C in| je — C ik gn € iniak — Caxlan) 19) 


It is well to examine the significance of these equations. There is no reason for us to 
assume that the strain rate is a tensor function (which would require that é,, satisfy certain 
conditions of integrability), for in the most general case, the é,, would be merely the com- 
ponents of a tensor field. However, the curvature tensor ,,,, is an intrinsic property of the 

*1j\kh 
of the spatial continuum H, in order that Gremain in 7. Thus, 4,, must be related to S, 
by an equation of the form 


material continuum G and must be related, as stated previously, to the curvature tensor S 


kh 


ne Gp O | Or‘ On Ov Oy’ 1) 
ee © Ou | Oy? dx! Or" Oz" ; 


We have assumed that 7 does not vary, that is, that the value of S,j,, at each point y* is 
constant. In particular, if S,.,=0, R,.,=0, and thus, of course 4,),, must also vanish. If. 
however, Sijx, 18 not zero, é jen Will be constant, or more strongly, zero, only if other special 
conditions hold 

The conditions of integrability of the rate of strain tensor must now be investigated. 
Using the relation for 4;),,, and the fact that R,,., is skew-symmetric in the indices of both 
pairs, we can write 


16) 


ind 


19) 


to 
Lin 
m- 
he 


kh 





However, the equations 
Ci hee jun t€ jh gen =Eijlen—islaa (22) 


are merely the conditions of integrability; the vanishing of either side of the above equation 
implies that if é,; is subjected to parallel transport around a closed curve, it will return to its 
initial value at the starting point, and thus, may be considered a tensor function of the coordi- 
nates x‘ rather than merely a tensor field. 

Thus, the rate of strain is integrable if any one of the following equations holds identically: 


€:j\en— €1j\me= O 
CR jn te Rosin =0 (23) 
( r+ Fite =O 


the last of the equations being equivalent to the requirement that jj), be skew-symmetric 
in its first two indices (it is already skew-symmetric in the second pair as a consequence of its 
definition). It should be noticed that these integrability conditions are distinct from the 
compatibility equations, and in particular, it is not necessary for 4;;., to vanish for é;; to be 
integrable. 

If, however, the strain rate is integrable, and we substitute the first of the above conditions 
in the compatibility equations, they take the simpler form 
( kh (Esp I €ixim) — (€ hi ik -€sx\ an) (24) 
In a great many practical cases, either the space is Euclidean, or some special relation 


may hold, such that ¢;;).,—0, which reduces the compatibility equations to the form 
(25) 


It is of some interest to see the forms which the last equations take when n<3. If n=1, 
all indices must be identical and we have, merely, 0=0. If n=2, only two indices can be 


distinct, resulting in only one independent equation 
€ xi iax t+ Cx i 2¢;) ky ixk, (26) 


indices not summed 
If n=3, three indices may be distinct, and two independent sets of, in all, six equations 


may be derived 
es; j ex %¢,, j ixk, (24) 


indices not summed (3 equations), and 


” ~ Fk, (28) 


indices not summed (3 equations 

The last two sets of equations, for n=3, are equivalent to those given by A. E. H. Love 
in Cartesian coordinates. Thus, our completely general form of the compatibility equations 
reduces to the more well-known form for the appropriate conditions. 

We note that, since « is a scalar parameter independent of z‘, derivatives of the infini- 
tesimal strain tensor de é,du may be written de é,du\,=€;\edu. Similarly, derivatives 
of the finite strain « Sé,du may be written ¢;|.= Sédule= S éyj|du. Thus, all the rela- 


tions and conditions derived above hold equally well for both the infinitesimal and finite strain 
tensors 
In conclusion, we can state the equations of integrability and the equations of compati- 


bility are two independent sets of conditions. Satisfaction’ of the:formerZimplies that the 
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space to which the deforming medium is assigned is Euclidean. Satisfaction of the latter 
implies the medium remains in that assigned space. The fact that both must hold for the 
solution of actual problems in elasticity should not lead to considering them as identical, for 
each has its own physical significance. 
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An Apparatus for Measuring the Piezoresistivity 
of Semiconductors’ 


a a 


R. F. Potter? and W. J. McKean 


\ detailed description is given of an apparatus and procedure designed to measure the 
piezoresistive effect in semiconductors over an extended temperature range. A tensile force 
up to 1 kilogram can be applied to the sample by means of a calibrated beam balance. The 


apparatus has been used for measurements on indium antimonide over the range 78° K to 
300° K, and tensile stresses of the order of 5X10? dynes per square centimeter can be applied 
) to samples that are cut in a special manner. 


In modern solid-state physics, phenomena such as | a form quite analogous to that for the elastic con- 
electrical conduction, Hall effect, and optical absorp- | stants. The II-constants are defined as follows: 
tion have been studied extensively because of their 


direct connection with a well-developed theory of 14.=Uy. 1 
| semiconductors. More recently several other effects A 

have been receiving an increasing amount of attention | a Mi2= Th, 29 (2) 
| by both the theorist and experimentalist; some of poX 

these are cyclotron resonance, photoelectromagnetic Ha THis us. 

effect, magnetoresistivity, and piezoresistivity. 
| The latter parameter has been studied for some | X,, is the applied stress, and p is the resistivity, 

time in connection with metals, but it is only in | where the cubic axes are taken as the reference axes. 

recent years that large anisotropic The change of resistance in a material is measured 


| comparatively 

changes in the resistivity with applied stress were | when a stress is applied in a given direction to a 
| measured for single crystals with cubic symmetry. | crystal while a constant current flows at a given 
Smith [1]° measured large effects for both n- and | angle to the stress direction. As in the case of 
p-type crystals of germanium and silicon at room | elastic constants for a cubic material, three inde- 
and liquid-nitrogen temperatures, thus adding fur- | pendent measurements must be made to determine 
ther to the evidence that the energy surfaces for | the elements of the tensor of formula 1. Smith 
electrons in these crystals are not spheres but | measured the longitudinal and the transverse effects 
ellipsoids. More recently, Morin, Geballe, and | in two samples of different orientations. The 
Herring [2] reported on measurements in which the | fourth determination served as a check on the other 
piezoresistive effect was produced with compressive | three. The longitudinal effect is defined as the 
stress. Herring [3] has discussed the above phe- | relative change observed in the resistance per 
nomena in terms of a theory of “simple many- | applied stress when both current and stress have the 
valley”’ semiconductors, which seems to describe the | same direction. The transverse effect is measured 
conduction bands for silicon and germanium, and | when current and stress are perpendicular to each 
possibly several III-V compound semiconductors. other. 

The piezoresistivity coefficients are the elements In converting the resistance to resistivity it is 
of a fourth-rank tensor; the form has been given by | necessary to consider the contraction due to strain. 
Smith [1] for cubic crystals of the O, and 7% classes, | The piezoresistivity coefficients are given by the 
which apply to crystals of the diamond and zine | formulas 
blende lattices, respectively, and is viven by 





» Apr AR po , , 
IL Ls —-— (§ Sie—Si2) (long.) 
II I] I] 0 0 0) 1 PS tee & 11 12 13 § 
ie Hy Tn @ @ @ » Apa AR». pg gies 
Il jo - Seer t (S12 +43) (trans.), 
Podn Rox 
Wy Wye Oy, 0 0) 0 (3) 
I] (1) 
CG 0 0 hy, 0 0 where Sj, are the elastic compliances for the orienta- 
tion used. Values for S,;, have been determined for 
0 0 0 O Iy 0 silicon [4], germanium [4], and indium antimonide [5). 
The II, are related to the elements of the tensor (1) 
0 0 0 0 0) IIgs =, bv 
This rese rch was supported by the United States Air Force through the Air Y nr 
Force Office of Scientific Research of the Air Research and Development Com- Ih, IT); 2ril 
mand under Contract CSO 670-53-12 
? Present Address—U. 8. Naval Ordnance Laboratory, Corona, Calif , , ‘ 
* Figures in brackets indicate the literature references at the end of this paper. | ITi2 IT}. + YII ? (4) 
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where Il’ IT, IT; _ IT 4s, and r and 7 are orientation 
factors related to the direction cosines of axes along 
the length and width of the sample as follows: 

Tr nil, 


[?mi +-nimi 


¥=GE + mims + njni. 

A program was recently undertaken at the Bureau 
to study the piezoresistivity of the cubic semi- 
conductor indium antimonide. (Concurrent studies 
were made of the magnetoresistivity of this material. 

Because of the nature and properties of the 
material, several considerations had to be taken in 
the design of an apparatus to measure piezoresist- 
ance. First was the matter of temperature. Reason- 
ably pure indium antimonide is intrinsic at room 
temperature, whereas it is certainly extrinsic at 
78° K; therefore, the temperature range down to that 
of liquid nitrogen had to be the minimum range 
studied. Indium antimonide is extremely brittle; a 
suitable method of clamping the sample had to be 
devised in order that known stresses could be applied 
along the axes of the sample. It is felt that the re- 
sultant design more than adequately solved the 
special problems posed by these properties. 

A tensile-stress system was selected as being the 
simplest to apply and to measure in a closed system. 
As shown in figure 1, a balance beam (S) graduated 


in 100-g intervals up to 1 kg was devised and alopo | 
which a weight (U) could be moved. Counterbalane } 
weights (J and K) could be adjusted so that the beam 
balanced with zero load. The load was transmitted to 
the sample by means of a thin piano wire (N), whieh 
was coaxial with an inner stainless-steel tube (Q), 
which in turn was coaxial with the outer stainless. 
steel tube (P). As the load was applied or removed. 
the wirc suffered a plastic strain, which could be com- 
pensated by moving the fulerum (R) up or down by 
means of a nut (D). A given load could be applied 
or removed quickly by turning the handle (C 
which activated the attached cam arrangement. 
Figure 2 shows a schematic cross section of the 
sample holder. The problem of applying pure tensile 
stress to the sample was solved as follows. <A mag. 
netostrictive machine was used for cutting flat 
specimens into a suitable form for resistivity measure- | 
ments by the standard four-probe method; a special 
die was made to produce samples of the shape shown 
in figure 3. Each end of the sample was placed ina | 
clevis with a pin through the eye (see fig. 2). The 
upper clevis (L) was attached to the piano wire; and 
the bottom clevis (S) was the upper part of a screw, | 





A spherical-shaped nut (V) on this screw seated ina 
mating socket allowed the sample to aline its axis in 
the direction of the applied tensile stress. This 
arrangement also permitted further adjustment in | 
the vertical position of the sample. It was eventually 
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Figure 1. Schematic drawing of the balance-beam mechanism. | 
The entire system is enclosed with the sample. The balance-arm rider and the vertical position of the fulerum can be adjusted by external manipulation. For 
dentification, the components are lettered as follows: A, Rider pickup arm; B, sliding O-ring vacuum seals (2): C, quick tension release; D, fulcrum height-adjusting 
screw; E, keyway; F, key; G, knife edge (hardened steel); H, tension wire bushing; I, Teflon insulator: J, zero-adjusting weight (coarse): K, zero idjusting weight 
(fine); L, opening to vacuum line; M, static ring vacuum seals (2); N, tension wire (0.006-in.-diam. piano wire); O, upper guide (Teflon); P, outer stainless-stee 
tubing (44-in. o. d. by 0.020-in. wall); Q, inner stainless-steel tubing in. o. d. by 0.004-in. wall); R, fulerum (hardened steel); 8, beam; T, Pyrex-glass tubing; U, 
rider; V, detents (11, 100 g apart); W, front cover for balance mechanism; X, gla nspection port; Y, 14-lead hermetic seal. 
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— along determined that a sample of this type and size (ap- 





erbalance | HUH Y proximately 1072 em? cross section) can be safely re- 
the beat y q \ t~< A cycled many times with loads up to 500 g. 
vues Yuin R . An aluminum shield may be placed around the 
1, Whieh Yard d sample holder. On the outside of this shield is 
tube | Q), ya | hel C wound a Kanthal wire heater (U) for temperature 
Stainless. 1 | HW control. A brass can (T) was soldered in place with 
removed, a A i" Wood’s metal; the entire interior could be evacuated 
| be com. Aw t or filled with a gas of suitable pressure. For use at 
down by J Y) i 5 ' temperatures between 78° and 300° K a Pyrex en- 
© fpPlied p i F velope was placed about the sample can and the 
ldle Y 


P stainless-steel tube (see fig. 4). This envelope was 
attached to a roughing pump, and the pressure of the 
gas inside varied from 1 atm to a few microns of mer- 
cury; this allowed the heat leak from the sample to be 
controlled when a suitable refrigerant bath was placed 
a the envelope. (For lower te mperatures, e. g., 
Q 1° or 20° K, the envelope can be replaced by a Dewar 
. cle be sth ‘d for holding liquid helium or liquid hydrogen. 
Figure 4 shows the assembled apparatus with the 
Py rex envelope. 

Measurements on both n- and p-type indium anti- 
monide have been made at temperatures between 
78° and 300° K, and the results are reported else- 
where [6]. The resistances were measured by the 
potentiometric method, using a type K-2 potenti- 
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Figure 2. Schematic drawing of the sample and holder 





By means of the pins through the sample ends, the sample axis lines up with the 
pplied stress For identification the components are lettered as follows 
Outer stainless-steel tubing (44-in. o. d. by 0.020-in. wall); B, tension wire (0.006 

liam. piano wire); C, inner stainless-steel tubing e-in. 0. d. by 0.004-in 
wall); D, top; E, vacuum-tight removable solder joint Wood’s metal); F, | \ 
thermal insulator; G, stainless-steel screws; H, frame holder; I, lower guide 

Teflon);"J, heater leads; K, heater retainin crew: L, upper clevis; M, clevis 
pin; N, sample; O, contact leads; P, thermocouple; Q, frame; R, clevis pin; 5, 
wer clevis; T, outer can; U, heater; V, lower spherical adjusting nut W, Teflon 


tion; X, heater bindir post 





mn. For 
djusting 
y weight 
less-steel 
bing; U 





FicuRE 4. Assembled apparatus with the Pyrex envelope. 


° : he frame supports associated vacuum lines, Vv ilves, meters, and pressure gages. 
Figure 3 Typical sample used with the apparatus. It can also support a double Dewar for measurements in the liquid-helium range. 
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ometer a breaker-type d-c amplifier in place of the 
galvanometer. The output of the amplifier was fed 


to a driven chart recorder, so that as the load was 
changed the incremental change in resistance was 
recorded by the unbalanced potentiometer. Checks 
were made to insure that currents due to the un- 
balance has negligible effects. Effects due to 


changes in the ambient temperature were observed 
on the recorder and the necessary corrections made. 

Samples were prepared in the following manner. 
Disks 2 mm thick and about 1 in. in diameter were 
sawed from single-crystal material with (001) axes 
normal to their flat surfaces. These were lapped to 
a thickness of 1.0 to 1.1 mm, taking care that both 
sides were parallel. The disks were then mounted 
on plate glass with glycol phthalate. An X-ray 
Laue back-reflection photograph was taken normal 
to the exposed side of each disk and analyzed. 
From the analysis the axes of the specimens were 
selected to lie along a particular orientation. A 
rotary milling table was used on the magneto- 
strictive machine so that the sample could be ac- 
curately oriented with respect to the tool. Indium 
electrodes were soldered to the four side arms and 
along the sides of the dumbbell-shaped ends. For 
the longitudinal effect, the end electrodes were used 
for the current leads, and two side arms were used 
for the potential leads. For the transverse measure- 
ments, the areas between a pair of arms on each side 
of the specimen were copper plated by applying a 
copper-sulfate solution momentarily and then rins- 
ing, taking care that no short circuit across sn sample 
occurred. The measurements were made by having 
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one current lead and one potential lead opposite » 


another such pair while the stress was applied alo 
the sample axis. As pointed out by Smith [1], th 
current lines are not straight the sampk 
Hence, a correction must be made if the piezoregis. 
tive coefficients deduced from transverse and long). 
tudinal measurements show appreciable differences 

It is felt that the apparatus has solved the special 
problems that are inherent in this type of measure. 
ment and a material like indium antimonide. De. 
termination of the longitudinal effect by means of 
this apparatus if used in conjunction with hydro. 
static-pressure measurements (cf. Long [7], Keyes 
[8], or Bridgman [9]) would eliminate the need for 
the transverse measurements. With a slight modi. 
fication of the sample holder to insure a fixed orienta. 
tion of the sample with respect to a magnetic field, 
one could also determine changes of the Hall eo. 
efficient and magnetoresistivity caused by applied 


stresses. 


across 
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